AD-R169  928  DEVELOPMENT  OP  FINITE  ACTIVE  CONTROL  ELEMENTS  FOR  LARGE  1/ 
FLEXIBLE  SPACE  ST.  .  (U>  MASSACHUSETTS  INST  OF  TECH 
CAMBRIDGE  MR  SPACE  SVSTEHS  LAB  D  M  HILLER  ET  RL. 
UNCLASSIFIED  B2  JUN  83  HIT-SSL-6-8S  RF0SR-TR-86-8488  F/O  28/11  NL 


I  »!** <Li  L*  U  Li-*,yi  f»|  ,«|  .<4 


unclassified 


SECURITY  CLASSIFICATION  OF  this  page 


REPORT  DO 


U.  REPORT  SECURITY  CLASSIFICAI  ION 

unclassified 


2«  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DEC  LASS  IF  I  CATION  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERtS) 

MIT-SSL  #6-85 


II 


3.  DiSTfffBUTlON/AVAILABlLITY  OF  REPORT 

unclassified /unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERIS) 

AFOSR*TR-  b  0  -  o  4  «0 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6B.  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

Space  Systems  Lab.,  MIT  tifoppuca bie> 

Dept.  Aero. and  Astro.  AFOSR 


6c.  ADORESS  |C|».  Sloie  end  /.IP  Code,  j  7t>.  ADDRESS  /City.  Stott  and  ZtPCodei 

37—361,  MIT  ca  aa  j-.  Bo 1 1  i n g  AFB 

Cambridge,  MA  02139/\L/—  |  09  9 28Washin8ton .  DC  20332-6448 


contract  <'!F49620-83-K-0026 


3a.  NAME  OF  FUNDING  SPONSOR  NO 

8b.  OFF  ICE  SYMBOL 

ORGANIZATION 

(If  applicable / 

AFASF 

Bolling  AFE 

Washington,  DC  20332-6448 

11  TITLE  (Include  Security  Classification/ 

(see  reverse) 

PROGRAM 
ELEMENT  NO. 

61102  F 


PROJECT 

TASK 

WORK  UNIT 

NO. 

NO. 

NO 

2302 

B1 

12.  PERSONAL  AUTHORIS) 

Miller,  David  Veslev 


13a  type  OF  REPORT 
annual 


->  ,  T  AT  I  ON 


Crawley,  Edward  Francis 


'.jmr  T0  8t/5/i«r«»^  ”'o,t  ,v'  ■ 


FILE  COPY 


17 

COSAT l  CODES  | 

FIELD 

GROUP 

SUP  GR  1 

18  SUBJECT  TERMS  /Continue  on  reverse  if  necessary  and  identify  by  block  number i 

active  control;  distributed,  hierarchic  control; 
flexible  space  structures;  distributed  processing; 
active  damping. 


19  ABSTRACT  iContmue  on  reverse  if  necessars  and  <drnfi/>  by  block  number< 

P-f  » 

Four  main  objectives  were  achieved  through  the  research  presented  in  this  report . First , 
a  quasi  free-free  beam,  simulating  a  flexible  space  structure,  was  equipped  with  inertial 
proof-mass  actuators  and  sensors  capable,  in  principle,  of  functioning  in  the  space  environ¬ 
ment,  i.e.,  space-realizable.  -Secondly , -^tuning  rules  were  derived  which  determine  the 
optimal  actuator  passive  stiffness  and  damping  which  minimizes  the  control  effort  required 
while  increasing  the  modal  damping  in  the  structure.  ■Tbir41y^> active  control  using  two 
local  and  one  component  level  processor  was  demonstrated.  Lastly,  multi-input,  single  out¬ 
put  collocated  feedback  tests  were  performed. 

Optimal  passive  vibration  absorber  designs  were  derived  to  provide  maximum  structural 
damping.  Theoretically the  addition  of  an  absorber  mass  equaling  0.5%  of  structural  mass 
can  result  in  a  single  mode  structural  damping  ratio  of  5%.  Analysis  of  multimode  damping 
using  a  single  absorber  indicated  that  the  absorber  stiffness  should  be  tuned  to  the  lowest 
mode  in  order  to  maximize  the  achievable  damping  in  all  the  modes.  {see  reverse) 


21  ABSTRACT  SECURITY  CLASSIFICATION 


20  Distribution  availability  of  abstract 
UNCLASSIFIED  UNLIMITE  o  £  SAME  AS  RPT  _  OTIC  USERS  D 


22*  NAME  OF  RESPONSIBLE  iNDiViOuAl 

Dr.  Anthony  Amos 


UNCLASSIFIED 


22B  TELEPHONE  number 
l Include  Area  Codr, 

202-767-4937 


I  OF  1  JAN  73  IS  OBSOLETE 
l 

Jt  . 


22c  OFF  »CE  SYMBOL 

AFOSR/NA 


unclassified _ , _ _ 

SECURITY  CLASSIFICATION  Of  this  pace  _ 

11)  Development  of  Finite  Active  Control  Elements  for  Large  Flexible  Space  Structures. 

>>•  ’\v 

19)  ^Passive  actuator  components  and  active  feedback  gains  were  derived  simultaneously 
and  sequentially  yielding  identical  results  whose  passive/components  equal  those  of  the 
passive  vibration  absorber. 

Passive  absorber  and  active  feedback  damping  tests  were  performed.  A  single  mode 
damping  ratio  of  4.2%  was  achieved  by  adding  2.3%  of  structural  mass:  77%  of  which  corres¬ 
ponds  to  'dead'  absorber  mass.  Multi-input,  single-output  control  provided  damping  ratios 
ranging  from  2%  to  3%.  This  functioning  collocated  control  element,  composed  of  actuators, 
sensors  and  local  and  component  level  processors,  constitutes  the  first  stage  of  research 
and  experimentation  into  distributed,  hierarchic  active  control  of  elastic  structural 
behavior.  - 


w  -w ■»-  ■ 


AFOSR-TR-  86-  0  480 


development  of 

finite  active  control  elements  for 
large  flexible  space  structures 

David  V.  Miller 
Edward  F.  Crawley 

MIT-SSL  ? 6—85 


SPACE  SYSTEMS  LABORATORY 
DEPT  OF  AERONAUTICS  AND  ASTRONAUTICS 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
CAMSRIDOE.  MA  OCT3B 


i 

<*.  \  ’  !•; 


i’. '  L  ’  i'  Bast  ; 

•  1  iffilti 


'  V  O  v"  *  •  vr*  •  -  '  v"  ■  '  -  •  v-v'Ov’v*  ^  « 


TABLK  Of  CONTENTS 


Report  Documentation  Page .  p.  1 

Table  of  Content?.  .  p.  3 

List  of  Tables .  p.  5 

List  of  Figures .  p.  7 


1.0  Research  Overview 


p.  11 


2.0  Design  Optimization  of  Passive  Vibration 

Absorbers  and  Inertial  Actuators  . 

2.1  luiimg  a  Single  DOF  Passive  Vibration 
AbsorLci  to  a  Single  DOT  Structural  Model 

2.1.1  bu-udy  State  Design  Criterion  .  . 

2.1.2  Pole  Placement  Design  Criterion 

2.1.3  Weighted  Quadratic  Cost  Design 

Criterion  . 

2.2  Tuning  a  Single  DOF  Passive  Vibration 
Absorber  to  a  Two  DOF  Structural  Model 

2.3  Sing'll  Structural  DOF  with  a  Single  DOF 

Active  Control  Actuator  . 

2.4  Conclusions  . 


p.  18 


P-  2A 


p.  25 


P.  34 
1-  3P 

p.  40 


p.  5C 


p.  56 


3.0  Control  Processor  Architecture  .  p.  58 

3.1  Architecture  for  a  Distributed  Hierarchy 

of  Control .  p.  60 

3.2  Compatible  Hardware  Design  Options  .  p.  64 

3.3  Present  Processing  Capability  .  p.  67 

3.4  Second  Generation  Processing  Network  ....  p.  69 

4.0  Experiment al  Hardware  and  Software  Desigr.  ....  p.  71 

4.1  Test  Article .  p.  71 

4.2  Actuator  Implementation .  j  .  73 


ir 


ty  Codes 

jod/or 

spuria! 


4.3  Sensor  Implementation  . 

4.4  Signal  Processing  Implementation  .  .  .  . 

5.0  Control  Experiments  . 

5.1  Design  Considerations  for  Experimental 

Phase  . 

5.2  Single  Input,  Single  Output  (SISO) 

Collocated  Controller  . 

5.3  Multiple  Input,  Single  Output  (MISO) 

Controller  . 

5.4  Tests  Using  Passive  Vibration  Absorbers  . 

5.5  Damping  Using  Active  Actuators  and 

Vibration  Absorbers  . 

6.0  Conclusions  . 

References  . 

A. O  Appendix  A:  Optimizing  the  Passive  Parameters 

of  Passive  Vibration  Absorbers  and  Active 
Control  Actuators  . 

A. 1  System  Models  . 

A. 2  Quadratic  Cost  Optimization  of  Parameters 

A. 3  Results  and  Evaluation  . 

A.  4  Conclusion  . 

B. O  Appendix  B:  Information  Transport  in 

Structures  . 

B. l  Discretization  of  Local  Elastic  Behavior 
B.2  Wave  Propagation  Model  Formulation  .  .  . 


LIST  OF  TABLES 


2.1  p.  145  Eigenvectors  for  Two  DOF  Model  for  a  Mass  Ratio  of  fi  = 

0.02. 

2.2  p.  145  Parameter  Values  for  Tuning  to  the  Low  and  High  Mode 

of  a  Two  Mode  Structural  Model. 

2.3  p.  146  Eigenvectors  and  Frequencies  for  Different  Absorber 

Frequency  Ratios.  {Three  DOF  Model). 

2.4  p.  147  Actuator  Design  Cases  Studied  in  Chapter  2  and 

Appendix  A. 

4.1  p.  148  Modal  Characteristics  of  a  Free-Free  Beam. 

4.2  p.  149  Force  Actuator  Specifications  and  Components. 

4.3  p.  150  Sensor  Specifications. 

4.4  p.  150  Signal  Processing  Equipment  Specifications. 

5.1  p.  151  Open  and  Closed-Loop  Characteristics  of  Controlled 

Modes 

5.2  p.  151  Analytic  Results  for  Optimally  Tuned  Vibration 

Absorber 

5.3  p.  151  Experimental  Vibration  Absorber  Results  for  Eighth 

Mode  Using  a  Tip  Absorber  at  Various  Frequency  Ratios. 

5.4  r  .  152  Summary  of  Experimental  Vibration  Absorber  Results  for 

Modes  lesled. 

A . 1  p.  153  Parameter  Optimization:  Initial  Conditions. 

A. 2  p.  153  Quadratic  Cost  Versus  Frequency  Ratio  and 
Nondimensional  Damper. 

A. 3  p.  153  Optimal  Passive  Parameters. 

A. 4  p.  154  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  Ip 
(Passive  Vibration  Absorber). 

A. 5  p.  154  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  la 
(Control  Actuator). 

A. 6  p.  154  "COARSE”  Quadratic  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber). 

A. 7  p.  155  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator). 

A. 8  p.  155  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber). 

A. 9  p.  155  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator). 

A. 10  p.  156  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20p 
{Passive  Vibration  Absorber). 

A. 11  p.  156  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator). 

A. 12  p.  156  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  2Ip 
(Passive  Vibration  Absorber). 

A. 13  p.  157  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  21a 
(Control  Actuator). 

A. 14  p.  157  "FINE"  Quadratic  Cost  Table.  Cost  Map,  System  lp 
(Passive  Vibration  Absorber). 

A. 15  p.  157  "FINE"  Quadratic  Cost  Table.  Cost  Map,  System  la 
(Control  Actuator). 


A. 16  p. 
A.  17  p. 
A.  18  p. 
A. 19  p. 
A. 20  p. 
A. 21  p. 
A. 22  p. 
A. 23  p. 


158  "FINE"  Quadratic  Cost  Table. 

(Passive  Vibration  Absorber). 
158  "FINE"  Quadratic  Cost  Table. 
(Control  Actuator). 

158  "FINE"  Quadratic  Cost  Table. 
(Passive  Vibration  Absorber). 

159  "FINE"  Quadratic  Cost  Table. 
(Control  Actuator). 

159  "FINE"  Quadratic  Cost  Table. 
(Passive  Vibration  Absorber). 

159  "FINE"  Quadratic  Cost  Table. 
(Control  Actuator). 

160  '  flNt  '  uuadratic  Cost  Table. 
(Passive  Vibration  Absorber). 

160  "FINE"  Quadratic  Cost  Table. 
(Control  Actuator) . 


Cost  Map, 
Cost  Map, 
Cost  Map, 
Cost  Map, 
Cost  Map, 
Cost  Map, 
Cost  Map, 
Cost  Map, 


System  20p 
System  20a 
System  20p 
System  20a 
System  20p 
System  20a 
System  2Ip 
System  2 la 


B.l  p.  161  Fourth  Derivative  Finite  Difference  Operators. 


6 


LIST  OF  FIGURES 


2.1  p. 


2.2  p. 

2.3  p. 


2.4  p. 


2.5  p. 


2.7  p. 

2.8  p. 


2.9  p 

2.10  p 


2.11  p. 


2.12  p. 

2.13  p. 


162  Inertial  Dampers. 

A.  Passive  Inertial  Vibration  Absorber. 

B.  Electromagnetic,  Inertial-Reaction  Actuator 
with  Passive  Components. 

163  Single  Degree  of  Freedom  Actuator  Coupled  to  a  Single 
Degree  of  Freedom  Representation  of  a  Structural  Mode. 

164  Magnitude  of  Response  (X/dg,p)  Versus  Nondimensional 

Frequency  (u/pq)  for  a  Typical  Value  of  Mass  and 
Frequency  Ratio. 

165  Magnitude  of  Response  (X/dg^,)  Versus  Nondimensional 
Frequency  (w/po)  for  the  Optimal  Choice  of  Frequency 

Ratio  Given  a  Mass  Ratio  of  0.02. 

166  Nondimensional  Natural  Frequencies  (7)  of  the  Two 

System  Mooes  as  a  function  of  Ab;  orber  Frequency  Katii 
16; . 

'■  '<1  lude  o'  Per  poors  Versus  Nondiaensiunol 

Frequency  (u/p^)  for  the  Optical  Cnoice  of  Frequency 

Ratio  Given  a  Mass  Ratio  of  0.2  (Compare  with  Fig. 
2.4). 

168  S-Plane  Representation  of  the  Poles  of  the  Vibration 
Absorber  as  a  Function  of  Increasing  Damper  Value. 

169  Comparison  of  Optimal  Steady-State  and  Pole  Placement 
Design  Characteristics. 

A.  Peak  Magnitude  of  Response  (X/dg^)  Versus 

Nondimensional  Frequency  (w/po). 

B.  Time  Response  to  Initial  Velocity. 

170  Single  DOF  Actuator  Coupled  to  a  Two  DOF  Structural 
Representat ion . 

171  Magnitude  of  Response  (X/dg^,)  Versus  Nondimensional 

Frequency  (w/po)  for  an  Absorber  Tuned  to  the  Lower 

Structural  Mode  (Curve  1)  and  Tuned  to  Both  Modes 
(Curve  2 ) . 

172  Magnitude  of  Response  (X/dg^,)  Versus  Nondimensional 

Frequency  (w/pQ)  for  an  Absorber  Tuned  to  the  Higher 

Structural  Mode  (Curve  1)  and  Somewhat  Tuned  to  Both 
Modes  (Curve  2). 

A.  Response  of  All  System  Modes. 

B.  Detail  in  the  Vicinity  of  the  Low  Modal 
Frequency. 

C.  Detail  in  the  Vicinity  of  the  Higher  Modal 
Frequency. 

174  Relative  Modal  Motions  of  Mass  Elements  for  Different 
Structure/Absorber  Frequency  Ratios. 

175  Comparison  of  the  Frequency  Response  of  Several 
Options  of  Active,  Passive  and  Active/Passive, 


2.14  p. 


3.1  p. 

3.2  p. 

3.3  p. 


3.4  p. 

3.5  p. 

3.6  p 

3.7  p 

3.8  p 

3.9  p 

3.10  p 

4.1  p 

4.2  p 

4.3  p 


4.4  p. 


4.5  p. 


4.6  P. 

4.7  p. 

4.8  p. 
5. 1  p. 


5.2  p. 


Simultaneous  Design. 

A.  Passive  Components 

B.  Active  Components 

176  Root  Locus  Plot  of  the  System  Poles  for  Various  Mass 
Ratios  Obtained  From  the  Different  Optimal 
Formulations. 

177  Hierarchic  Processing  Network. 

178  Levels  of  Modeling  and  Information  Feedback  in  a 
Distributed,  Hierarchic  Control  Format. 

179  Options  for  the  Selection  of  Bandwidths  and 
Interactions  of  the  Various  Levels  of  Control. 

A.  Control  Option  Providing  Extended  Bandwidth  in 
cne  low  autnority  Controller, 
is.  Overlapping  Low  and  High  Authority  Control. 

180  Single  Bus  Hardware  Arrangement  for  Multiprocessor, 
Hierarchic  Architecture. 

180  Multibus  Hardware  Arrangement  for  Multiprocessor, 
Hierarchic  Architecture. 

181  Interrupt  Flow  Chart  for  Single  Bus  System. 

182  Interrupt  Flow  Chart  for  Multibus  System. 

183  First  Generation  Processing  Hardware. 

184  First  Generation  Control  Software  Flow  Chart. 

185  Second  Generation  Processing  System. 

186  Experimental  Hardware  Configuration. 

187  Swinging  Parallelogram  Beam  Suspension. 

188  Electromagnetic,  Inertial-Reaction  Actuator 

Conf  i  gnrs>*  ' 

A .  Tor q ue  C on  figuration. 

B.  Force  Conf  itfurat  ion 

C.  Side  View. 

189  Models  of  Actuator  Configurations  Coupled  to  a  Single 
DO!  Structural  Representation. 

A.  Model  of  a  Single  DOF  Torque  Actuator 
Controlling  a  Single  Rotational  Mode  of  a 
Structure. 

B.  Model  of  a  Single  DOF  Force  Actuator 
Controlling  a  Single  Translational  Mode  of  a 
Structure. 

189  Block  Diagram  of  a  Digital  Double  Integrator. 

A.  Open-Loop  Configuration. 

B.  Closed-Loop  Configuration. 

190  Frequency  Response  of  Stabilized  Double  Integrator. 

191  4.tock  A??  vgj of  To",  ll  integration  and  Filtering 
System. 

192  Frequency  Responses  of  Single  and  Double  Integrating 
and  Filtering  Systems. 

193  Actuator  Effectiveness  as  a  Function  of  Position  on 
the  Beam  as  Indicated  by: 

A.  The  Sum  of  the  Squares  of  the  Controlled  Mode 
Shapes . 

B.  The  Product  of  the  Squares  of  the  Controlled 
Mode  Shapes. 

194  Root  Locus  of  the  Two  DOF  Model  when  Rate  Feedback  is 


5.3  p. 

5.4  p. 


5.5  p. 


5.6  p. 


5.  '  p 


p. 


5.9  p. 


5.10  p . 


5.11  p. 


5.11:  p. 


5.14  p. 

5.15  p. 

5.16  p. 


5.17  p. 


A.  I’urely  Beam  Rale  Feedback  with  no  Passive 
Damper  (.p  =  0). 

B.  Beam  Rate  Feedback  with  the  Passive  Actuator 
Damper  Included  (p  >  0). 

C.  Stabilization  of  Case  B  through  Additional 
Actuator  Rate  Feedback. 

195  Experimental  Instability  Introduced  by  High  Gain  Beam 
Rate  Feedback  (See  Fig.  5.2b). 

196  Experimental  Transient  Decay  of  Seventh  Mode. 

A.  Open-Loop. 

B.  Closed-Loop. 

197  Experimental  Transient  Decay  of  Eighth  Mode. 

A.  Open-Loop. 

B.  Closed-Loop. 

198  Experimental  Transient  Decay  of  Ninth  Mode. 

A.  Open-Loop. 

B.  Closed-Loop. 

199  Experimental  Transient  Decay  of  Tenth  Mode. 

A.  Open- Loup, 
i. .  Closed -Loop. 

200  Experimental  Irarisient  Decay  of  Eleventh  Mode. 

A.  Open-Loop. 

B.  Cl  used -Loop. 

201  Experimental  Transient  Decay  of  Twelfth  Mode. 

A.  Open-Loop. 

B.  Closed -Loop. 

202  Experimental  Transient  Decay  of  Thirteenth  Mode. 

A.  Open-Loop. 

B.  Closed- Loop. 

203  Experimental  Transient  Decay  of  Fourteenth  Mode. 

A.  Open  Loop, 
b.  Closed-Loop. 

204  Experimental  Forced  Response  Versus  Sinusoidal 
Excitation  for  Open  and  Closed-Loop  Controller . 

A.  Open-Loop. 

B.  Closed-Loop. 

205  Experimental  Fassive  Tuning  of  Vibration  Absorber  to 
the  Eighth  Mode. 

A.  Actuator  Frequency  is  2.60  Hz. 

B.  Actuator  Frequency  is  3.31  Hz. 

C.  Actuator  Frequency  is  4.50  Hz. 

D.  Actuator  Frequency  is  5.68  Hz. 

E.  Actuator  Frequency  is  5.93  Hz. 

208  Experimental  Passive  Tuning  to  the  Seventh  Mode. 

208  Experimental  Passive  Tuning  to  the  Ninth  Mode. 

209  Experimental  Forced  Response  Versus  Sinusoidal 
Excitation  for  Untuned  and  Tuned  Absorbers. 

A.  Open-Loop. 

B.  Closed- Loop. 

210  Comparison  of  Experimental  Results  and  Theoretical 
Predict  ions. 

A.  Minimum  Achieved  Damper  Settings  Near  Optimal 
Frequency  Ratios. 

B.  Comparison  of  Experimental  and  Predicted 


Damping  Ratios. 

Optimal  Damper  Characteristics  Versus  Mass  Ratio. 
Comparison  of  Damping  Tests  Performed  on  the  Ninth 
Mode. 

A.  Free  Decay. 

B.  Passively  Tuned  Decay. 

C.  Actively  Controlled  Decay. 

D.  Actively  Controlled  and  Passively  Tuned  Decay. 
Damper  Models. 

A.  Passive  Absorber. 

B.  Actuator. 

Slriirtnra'  Ponrocont-Blinnc  flonnled  t.O  One  DOF 
Absorber , 

A.  System  1  :  One  Mode  Structure  with  Absorber. 

B.  System  20:  Two  Mode  Structure  with  Absorber  on 
Outer  Mass. 

C.  System  21:  Two  Mode  Structure  with  Absorber  on 
Inner  Mass. 

Closed  Loop  Poles  of  System  1  Optimized  for  Z  - 

o 

10  1  0  o.T. 

Closed  Loop  Poles  of  System  20  Optimized  for  Z  = 

o 

[0  1  0  1.62  0  0jT. 

Closed- Loop  Poles  of  System  20  Optimized  for  Z  - 


Closed  Loop  Poles  of  System  20  Optimized  for  Z  - 
; 0  -0.618  0  2.62  0  0:T. 

Closed  Loop  Poles  of  System  21  Optimized  for 

(0  -0.618  0  2.62  0  0T. 

Structural  Discretization. 

bandedness  in  Discrete  Representations  of  Elastic 
St  ruet ures . 

A.  handedness  in  Fourth  Derivative,  Finite 
Difference  Operators. 

B.  Comparison  of  Bandedness  in  Finite 
Difference  Operators  and  Flexibility 
Influence  Coefficients. 

A  Finite  Bea™ 

Stillness  Terms  Derived  Using  a  Flexibility 
Influence  Coefficient  Formulation. 

Bandedness  in  Finite  Element  Representation  of 
Structural  Stiffness  as  Illustrated  by: 

A.  Adjacent  Terms  Corresponding  to  Unit 
Displacement . 

B.  Adjacent  Terms  Corresponding  to  Unit  Rotation. 
Characteristic  Wave  Propagation  Length  Versus  Damping. 
Ratio  of  Characteristic  Length  of  Propagation  and 
Structural  Length  Versus  Damping  in  Structural  Modes. 


1.0  RESEARCH  OVERVIEW 


The  long  term  goal  of  this  research  is  to  develop  techniques 
by  which  the  dynamic  elastic  behavior  of  large  space  structures 
can  be  effectively  modified  to  meet  mission  specifications 
primarily  through  the  incorporation  of  machine  intelligence.  By 
providing  active  structural  control,  through  the  means  of  widely 
distributed  sensing  and  actuation,  several  beneficial 
characteristics  can  be  attained.  This  report  discusses  these 
benefits  and  presents  solutions  to  some  of  the  problems 
encountered  in  practical  implementation  of  structural  control. 

Numerous  space  structures  may  require  some  form  of  active 
control.  The  lack  of  structural  stiffness  and  foundation  promotes 
the  occurance  of  large  amplitude  structural  vibrations.  Because 
of  the  lack  of  energy  dissipation  mechanisms  prevalent  on  earth, 
these  motions  are  detrimental  to  the  performance  of  precision 
pointing  equipment  and  maintenance  of  a  micro-gravity  environment. 
Malfunctioning  of  critical  station-keeping  and  communication 
systems  might  also  occur  due  to  large  and  persistant  vibration 
levels. 

Active  structural  control  is  one  possible  solution  to  both 
reducing  the  magnitude  of  vibration  levels  and  dissipating  the 
associated  energy  (References  1,2  and  7).  The  sensors,  processors 
and  actuators  comprise  the  closed-loop  controller  for  the 
structure,  or  plant.  This  closed-loop  controller  would  introduce 
a  variable  set  of  state-dependent  forces  into  the  dynamics  of  the 
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structure.  The  designer  can  tailor  the  dynamics  of  the  plant 
through  the  proper  selection  and  configuration  of  these  control 
components. 

The  modifications  of  particular  interest  are  those  that  are 
unattainable  through  passive  techniques.  Not  only  Hookean,  but 
non-Hookean  behavior  can  be  induced  in  the  structure.  This  allows 
the  designer  versatility  in  modifying  the  dynamics  by  techniques 
such  as  time-varying  gains,  dynamic  compensations  and 
noncollocated  state  feedback.  In  addition,  structural  response 
can  be  easily  modified  during  construction,  and  in  the  event  of 
mission  redefinition  or  component  failure.  Futhennore,  many 
actuators  contain  devices  which  are  capable  of  dissipating  energy 
in  the  space  environment.  While  these  benefits  of  active  control 
warrant  its  further  investigation,  tne  implementation  of  a 
structural  controller,  capable  of  functioning  in  space,  presents 
several  unique  problems  that  must  be  addressed. 

One  problem  involves  designing  a  space-realizable, 
closed-loop  control  system  for  experimental  testing  and  later  for 
orbital  application.  To  be  space-realizable,  all  control 
components  must  be  capable  of  functioning  in  the  space 
environment.  In  such  a  system,  all  sensing  and  actuation  must  be 
performed  independent  of  the  laboratory  frame.  For  this  reason, 
only  inertial  or  relative  motion  actuators  and  sensors  can  be 
considered. 

Modeling  of  large  space  structures  presents  a  second 
difficulty.  Their  large  and  extremely  flexible  nature  presents  a 
couple  of  problems.  First,  generating  Linear  Quadratic  Regulator 
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feedback  controllers  or  dynastic  state  estimators  requires  an 
accurate  model  of  the  structure.  Because  the  number  of  structural 
modes  is  potentially  infinite,  a  very  large  order  model  would  be 
required.  This  model,  even  if  generated  by  one  of  the 
sophisticated  techniques  available  today,  will  inaccurately 
represent  the  higher  frequency  modes.  Since  implicitly  both 
output  feedback  and  state  estimation  schemes  require  accurate 
models  to  insure  predictable  closed-loop  response  and  explicitly 
state  estimation  schemes  require  accurate  models  to  generate  state 
information,  these  inaccurate  structural  representations  will 
likely  cause  these  feedback  schemes  to  become  unstable  when 
implemented  on  the  actual  structure.  For  this  reason  an 
alternative  approach  to  the  modeling  and  control  problem  should  be 
investigated. 

Secondly,  a  large  space  structure  may,  depending  on  the 
geometric  properties  of  the  structure,  have  modes  that  are  close 
in  frequency.  This  can  pose  a  problem  in  the  rolloff  region  of 
the  controller  where  the  phase  and  gain  margins  become  critical. 
In  addition,  such  close  frequency  spacing  can  make  accurate  model 
generation  and  identification  more  difficult. 

The  goal  of  this  research  is  to  develop  a  control  technique, 
which  can  be  readily  used  for  any  structural  design  application, 
that  relies  on  limited  knowledge  of  basic  structural  behavior  and 
incorporates  space-realizable  components.  This  would  reduce  the 
dependency  on  accurate  high  order  models  which  are  inherent  to  the 
application  in  question. 

Schemes  will  be  investigated  which  provide  control  robustness 


(model  insensitivity)  as  well  as  accurate  low  order  modeling. 
Possible  solutions  to  the  actuation,  sensing,  modeling  and 
processing  speed  problems  will  also  be  discussed.  Actuators  and 
sensors  that  operate  on  inertial  elements  are  chosen  to  eliminate 
interaction  with,  or  reference  to,  the  laboratory  (earth)  frame. 

The  control  methodology  under  consideration  uses  hierarchic 
control,  through  distributed  processing,  to  perform  low  and  high 
authority  control.  Simple  distributed  processing  allows  control 
algorithms  to  be  processed  simultaneously  by  several  processors. 
This  enables  a  higher  net  processing  rate  than  use  of  a  single 
equivalent  processor,  and  in  turn,  allows  increased  controller 
bandwidth  and  therefore  the  control  of  higher  frequency  modes.  By 
arrangement  of  the  processors  in  a  pyramidal,  hierarchic  command 
structure,  additional  benefits  are  derived.  3y  using  hierarchic 
control,  different  control  casks  are  logically  distributed  among 
the  processors. 

The  logic  behind  this  hierarchic  control  is  to  distribute  the 
control  function  in  a  way  so  as  to  reflect  the  fundamental  dynamic 
behavior  of  the  structure.  Three  levels  of  processors  are 
envisioned.  These  levels  are  refered  to  as  the  zero  or  local 
level ,  the  one  or  component,  level  and  the  two  or  global  level 
processors.  These  three  levels  of  processing  parallel  the  three 
levels  of  detail  in  the  modeling  of  the  passive  structure:  the 
fundamental  local  material  properties;  the  internal  element 
behavior;  and  the  global  behavior  of  a  Finite  Element  Structural 
model;  respectively.  The  zero  or  local  level  processors  carry  out 
high  speed,  localized  low  authority  (almost  collocated)  control 
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using  a  limited  number  of  local  sensors  and  a  single  actuator. 

The  one  or  component  level  controller  uses  the  information 
and  actuators  associated  with  several  adjacent,  local  level 
controllers  to  perform  higher  authority,  noncollocated  control  on 
the  disturbances  of  wavelength  up  to  the  length  of  the  element. 
The  increased  length  over  which  the  actuators  and  sensors,  under 
the  control  of  the  element  processor,  act  allows  improved 
observation  and  control  of  longer  wavelength  (lower  frequency) 
modes.  The  term  controlled  element  refers  to  a  portion  of  the 
structure  that  is  controlled  by  one  component  level  processor  and 
its  associated  local  level  processors. 

The  two  or  global  level  would  perform  three  main  tasks. 
First,  it  would  act  as  the  high  authority  controller  for  the 
lowest  modes  since  these  modes  are  of  the  longest  wavelength,  are 
the  most  likely  to  be  excited  and  are  usually  most  critical  to 
mission  objectives.  Secondly,  a  global  level  processor  for 
systems  identification  and  a  global  processor  for  failure 
detection  and  control  reconfiguration  could  be  added  at  this 
level.  A  detailed  discussion  of  a  portion  of  the  global  level  can 
be  found  in  Reference  4.  Note  that  these  are  the  only  processors 
that  have  the  capability  to  observe  and  influence  the  structure 
along  its  entire  length. 

The  size  of  the  control  task  allocated  to  each  level,  i.e. 
the  number  of  sensor  signals  it  processes,  is  inversely 
proportional  to  each  level’s  required  control  rate.  The  local  or 
zero  level  controller  processes  a  limited  amount  of  local  sensor 
information  to  generate  a  control  command  for  an  actuator 
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controlling  short  wavelength,  high  frequency  waves.  At  higher 
levels,  more  noncollocated  information  is  used  to  generate 
commands  to  more  than  one  actuator,  thus  requiring  longer 
processing  delays  and  slower  control  rates.  Since  these  higher 
levels  control  lower  frequency  waves,  this  represents  an  ideal 
marriage  of  processor  task  and  controller  level  bandwidth.  In 
this  framework  of  hierarchic  control,  the  development  of  the 
'Finite  Control  Element*  will  be  discussed. 

The  contents  of  this  report  detail  the  initial  research  into 
the  experimental  implementation  of  'Finite  Control  Elements'  and 
its  application  to  low  authority  control.  A  discussion  of  the 
analytic  formulation  of  the  hierarchic  approach  can  be  found  in 
Reference  14.  The  questions  involved  in  the  implementation  of  the 
Finite  Control  Element  include:  actuator  selection;  actuator 
design  optimization;  sensor  selection;  selection  and  integration 
of  the  zero  and  one  level  microprocessors;  and  determination  of 
the  level  of  achievable  control. 

Chapter  2  discusses  design  of  an  inertial  actuator  and  a 
passive  vibration  absorber.  Several  Electromagnetic, 
Inertial-Reaction  Actuator  configurations  are  presented.  Tuning 
laws  are  derived  to  determine  actuator  designs  best  suited  for 
both  passive  damping  and  active  control.  In  this  chapter, 
actuator  design  options  are  investigated  which  provide 
effort-efficient  active  control. 

Chapter  3  gives  a  more  detailed  explanation  of  the  processing 
hierarchy  and  its  implementation.  The  details  of  the  data  links 
between  the  first  two  levels  of  controllers  and  the  microprocessor 
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configuration  used  in  the  experimental  phase  are  discussed. 

Chapter  4  discusses  the  design  of  the  experimental  test 
apparatus.  The  actuators  and  sensors  are  chosen  so  as  to  provide 
a  space-realizable  experiment  (i.e.  capable,  in  principle,  of 
functioning  in  space).  Test  article  design,  actuator 
construction,  sensor  selection  and  signal  processing  development 
are  discussed. 

Chapter  5  encompasses  the  experiments  that  have  been 
performed  to  date.  First,  some  preliminary  research  is  presented 
to  illustrate  the  goals  and  direction  of  the  laboratory 
investigation.  Next,  the  results  of  several  tests  are  compared  to 
expected  results  to  demonstrate  the  capabilities  of  the  system. 
Finally,  the  limitations  of  the  controller  are  discussed  along 
with  planned  future  improvements. 

Appendix  A  presents  actuator  optimization  and  state  feedback 
control  law  generation  using  weighted  quadratic  cost  and  Linear 
Quadratic  Regulator  formulations.  This  appendix  approaches  the 
same  problem  as  is  presented  in  Chapter  2;  but  from  a  control 
theory  perspective. 

Appendix  B  presents  some  of  the  initial  research  into  the 
concept  of  Finite  Control  Elements.  This  appendix  studies 
techniques  through  which  fundamental  structural  behavior  can  be 
exploited  in  a  control  environment.  Model  generation  for 
transverse  motion  is  developed.  Ways  in  which  the  computational 
burden  and  dependency  on  model  accuracy  can  be  reduced  are  used  to 
set  guidelines  for  the  development  of  control  algorithms,  sensor 
selection,  actuator  selection  and  actuator  placement. 


2.0  DESIGN  OPTIMIZATION  OF  PASSIVE  VIBRATION  ABSORBERS  AND 

INERTIAL  ACTUATORS 

This  chapter  discusses  design  considerations  for  two  devices 
intended  to  Modify  structural  behavior.  The  first  type  uses 
passive  Mechanises  to  dissipate  vibrational  energy.  The  second 
type  eaploys  active  feedback  to  alter  the  characteristics  of  a 
particular  structure.  The  passive  device  will  be  hereafter 
refered  to  as  a  passive  vibration  absorber.  The  device  that 
operates  in  conjunction  with  active  feedback  is  termed  a  control 
actuator.  In  the  following  sections,  one  device  of  each  type  is 
selected  and  analyzed  with  respect  to  its  ability  to  Modify 
structural  response.  For  each  of  these  two  types,  a  device  is 
selected  which  is  capable  of  operating  in  the  space  environment: 
i.e.  it  is  space-realizable.  The  ultiMate  goal  of  analyzing  both 
a  passive  vibration  absorber  and  active  control  actuator  is  to 
determine  whether  the  corresponding  concepts  can  be  incorporated 
into  one  device,  which  would  combine  the  attributes  of  each,  and 
therefore  could  achieve  structural  damping  performance  exceeding 
that  of  either  device  taken  separately. 

Many  types  of  control  actuators  are  available  for 
demonstration  purposes  in  the  laboratory.  However,  to  demonstrate 
a  space-realizable  controller,  the  forces  on  the  structure  must  be 
created  without  reacting  against  the  laboratory  frame  of 
reference.  In  other  words,  the  controller  must,  at  least  in 
concept,  be  functional  in  space.  Several  categories  of 


space- realizable  actuators  are  available.  Expendable  gas  jet 
thrusters  can  be  used  to  control  elastic  Bodes  by  accelerating 
expendable  mass.  Intrastructural  actuators,  such  as  piezoelectric 
crystals  or  hydraulic  pistons,  can  generate  control  forces  by 
pushing  against  adjacent  parts  of  the  structure.  Captive  mass 
actuators  can  influence  a  structure  by  reacting  against 
nonexpendable  mass  or  stored  momentum.  Captive  mass  actuators 
include  inertial  actuators  and  Control  Moment  Gyros.  Inertial 
actuators  operate  on  a  principle  of  momentum  exchange  between 
masses  and  the  structure.  Examples  of  inertial  actuators  are 
pivoting,  proof-mass  actuators  and  momentum  wheels.  Control 
Moment  Gyros  (CMG’s)  react  against  stored  angular  momentum.  Each 
of  these  actuators  influences  the  plant  in  proportion  to  some 
commanded  control  signal. 

Several  types  of  space-realizable  vibration  absorbers  can  be 
used  to  passively  dissipate  structural  energy.  Viscoelastic 
coatings  and  distributed  friction  dampers  can  be  used  to  dissipate 
vibrational  energy  directly  as  heat.  Intrastructural  dampers  and 
lossy  joints  generate  damping  forces  proportional  to  the  relative 
velocity  between  adjacent  elements  of  a  structure.  Captive  mass, 
inertial  vibration  absorbers  dissipate  relative  motion  between  the 
structure  and  a  proof  mass.  Each  of  these  devices  can  be  used  to 
modify  the  passive  plant  dynamics. 

From  amongst  these  various  concepts  for  space-realizable 
control  actuators  and  passive  vibration  absorbers,  this  chapter 
discusses  the  design  of  passive,  inertial  vibration  absorbers  and 
inertial  control  actuators.  Because  of  the  functional 


similarities  in  their  models,  the  possibility  is  created  of 
incorporating  both  concepts  into  one  device.  An  Electromagnetic, 
Inertial-Reaction  Actuator  is  defined  as  a  device  that 
electromagnet ically  accelerates  a  captive  mass,  hereafter  refered 
to  as  a  proof  mass.  By  creating  an  action  on  the  proof-mass,  an 
equal  and  opposite  reaction  is  exerted  on  the  structure.  The 
Electromagnetic,  Inertial-Reaction  Actuators  may  also  have  a 
passive  spring  and  damper.  The  passive,  inertial  vibration 
absorber  also  uses  a  captive  mass  to  create  a  relative  velocity 
across  some  damping  mechanism  located  between  the  structure  and 
the  attached  mass.  Figure  2.1  illustrates  the  similarity  between 
the  models  used  for  the  analysis  of  the  Electromagnetic, 
Inertial-Reaction  actuator  and  passive,  inertial  vibration 
'rh°  difference  between  the  two  devices  is  that,  while 
onlv  collocated  state  information  affects  the  behavior  of  the 
passive  absorber,  noncollocated  state  feedback  can  be  used  to 
drive  the  actuator’s  motor. 

Both  the  actuator  and  absorber  can  be  modelled  as 
spring-mass-damper  assemblies  (Fig.  2.1).  The  mass  element  is  the 
proof  mass,  which  can  represent  either  a  translational  mass  or 
rotational  inertia.  The  spring  element  is  a  centering  spring 
which  prsvj:: to  th.i  proof -.'xaor  fro*  developing  e  bias  displacement. 
The  damper  corresponds  to  an  adjustable  viscous  damper  and,  in 
addition,  represents  the  armature  friction  and  other  passive 
damping  sources  in  the  mechanism. 

In  the  actuator  model,  the  values  of  these  passive  parameters 


can  be  selected  to  help  reduce  the  amount  of  control  effort  needed 


to  achieve  a  desired  system  performance.  The  presence  of  these 
passive  components  also  enables  the  actuator  to  possess  the 
characteristics  of  a  passive  vibration  absorber. 

In  the  case  of  the  vibration  absorber,  the  parameters  may  be 
chosen  such  that  the  structure  is  damped  most  effectively. 
Understanding  the  limiting  cases  helps  clarify  what  is  meant  by 
effective.  If  the  device  in  Figure  2.1a  is  attached  to  an 
undamped  structure  and  the  damper  is  removed,  the  assembly 
contains  no  damping.  Alternately,  if  the  damper  is  set  infinitely 
strong,  the  proof  mass  is  locked  to  the  structure  and  again  no 
damping  is  introduced  into  the  structural  mode.  Since  the 
structure  is  undamped  in  both  limiting  cases,  there  must  exist, 
between  the  limiting  cases,  an  optimal  set  of  vibration  absorber 
parameters  which  introduce  the  greatest  degree  of  damping  into  the 
structure. 

The  principle  question  of  interest  in  this  chapter  is  as 
follows:  is  the  design  of  the  optimal  passive,  inertial  vibration 
absorber  identical  to  the  design  of  the  optimal  passive  components 
of  an  optimized  Electromagnetic,  Inertial-Reaction  actuator 
design.  If  so,  the  classical  vibration  absorber  analysis, 
presented  in  the  next  section,  can  be  used  to  design  the 
optimized,  inertial  reaction  control  actuator. 

Section  2.1  discusses  three  design  procedures  in  which  an 
inertial  vibration  absorber,  modeled  in  Figure  2.1a,  is  used  to 
damp  a  single  structural  mode.  These  procedures  differ  in  their 
definition  of  desired  performance  for  the  absorber/structure 
assembly.  The  absorber  parameters  are  selected  to  achieve  some 


desired  optimal  performance  in  each  of  the  three  cases.  This 
process  of  selecting  passive  absorber  parameters  will  hereafter  be 
refered  to  as  tuning.  The  three  resulting  absorber  designs  are 
then  compared  and  evaluated. 

The  first  design  procedure  analyzes  the  use  of  a  vibration 
absorber  to  alter  the  steady-state  response  of  a  single  structural 
mode  (dection  c.i.ij.  me  optimal  steady-state  solution  results 
in  the  passive  vibration  absorber  design  which  minimizes  the 
steady-state  response  of  one  structural  mode  subjected  to  a  broad 
band  disturbance.  This  absorber  minimizes  the  maximum 
steady-state  response  of  the  structure,  using  a  minimax  criterion. 

In  the  second  design  procedure,  the  optimal  pole  placement 
solution  yields  a  set  of  absorber  parameters  which  place  the 
system  poles  as  far  left  fn  the  complex  S-plrvne  as  possible 
(Section  2.1.2).  Bv  attaching  a  single  degree  of  freedom  (DOF) 
absorber  to  a  structure  modeled  with  one  mode,  the  complete  model 
results  in  a  structural  motion  which  exhibits  two  modes,  each 
represented  by  a  pair  of  complex  conjugate  S-plane  poles.  The 
objective  is  to  choose  absorber  parameters  which  place  these  poles 
as  far  left  in  the  S-plane  as  possible. 

The  third,  or  final  design  procedure  for  the  passive  inertial 
vibration  absoroe-*  rey.,T~,»:>  absorber  parameters  be 
determined  such  that  a  structural  response  to  a  specified  initial 
condition  minimizes  the  value  of  a  particular  quadratic  cost 
function  (Section  2.1.3).  The  cost  function  of  interest  is 
composed  of  weighted  quadratic  penalties  on  the  response  of 
various  portions  of  the  structure.  For  example,  the  penalty  could 
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be  placed  primarily  on  a  particular  system  mode,  or  all  of  the 
modes  could  be  evenly  penalized.  Since  energy  is  typically  the 
quantity  penalized,  a  cost  function  composed  of  the  weighted  sum 
of  the  squares  of  the  state  variables  is  used.  This  penalty 
function  is  integrated  forward  in  time  from  some  initial  time  and 
state  vector  to  some  arbitrary  final  time  (usually  infinity)  to 
yield  a  total  cost.  The  vibration  absorber  is  then  designed  to 
minimize  the  value  of  this  total  cost. 

Section  2.2  extends  the  steady-state,  vibration  absorber 
problem  to  that  of  an  absorber  attached  to  a  structural  model 
containing  two  modes.  This  is  done  to  provide  a  means  for 
extrapolating  the  absorber  design  rules  of  Section  2.1  to  the 
design  rules  for  an  absorber  that  is  to  be  used  to  damp  a  multiple 
degree  of  freedom  structure. 

Section  2.3  details  the  design  procedure  for  the  active, 
inertial-reaction  control  actuator.  In  this  section,  quadratic 
penalties  on  control  effort  are  added  to  the  quadratic  cost 
function  of  state  developed  in  Section  2.1.3,  and  a  set  of  passive 
actuator  parameters  and  active  feedback  gains  are  derived 
simultaneously.  The  passive  portion  of  this  simultaneous  design 
analysis  and  the  passive  absorber  of  Section  2.1  are  compared.  If 
the  passive  vibration  absorber  design  is  similar  to  the  passive 
portion  of  the  active  control  actuator  design,  the  optimal  passive 
absorber  can  be  used  as  the  passive  component  of  an  optimal 


control  actuator. 
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2.1  TUNING  A  SINGLE  DOF  PASSIVE  VIBRATION  ABSORBER  TO  A  SINGLE 

DOF  STRUCTURAL  MODEL 

This  section  develops  design  rules  for  selecting  the 
parameters  of  a  passive  vibration  absorber.  Rules  are  developed 
for  designing  aosoroers  based  on  three  design  criteria:  best 
steady-state;  optimal  pole  placement;  and  minimum  weighted 
quadratic  cost  performance.  The  three  resulting  designs  are  then 
compared. 

Since  a  modal  representation  of  a  structure  yields  decoupled, 
ordinary  differential  equations  expressed  in  generalized 
coordinates,  a  single  mode  was  deamed  the  simplest  model  which  can 
represent  tynirnl  structural  behavior.  Tn  each  of  the  three 
analyses,  the  rules  are  developed  for  tuning  the  absorber  model 
shown  in  Fig.  2.1a  to  a  structural  model  containing  a  single 
structural  mode.  Figure  2.2  is  the  resulting  model  analyzed  with 
a  single  DOF  absorber  and  a  single  structural  mode:  where  Kj  and 
Mi  are  the  modal  stiffness  and  modal  mass;  respectively.  Note 
that  throughout  this  section,  the  possibility  of  introducing 
active  feedback  control  through  the  electromagnetic  force  labled 
F ,  is  ignored  The  pass  v/ a  absorber  spring  (X^).  damper  (C)  and 
mass  (M2)  are  the  only  parameters  which  may  be  varied.  The  rules 
for  selecting  these  parameters  are  developed  below. 
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2.1.1  STEADY-STATE  DESIGN  CRITERION 


The  first  of  the  three  design  criteria  for  designing  a 
passive  vibration  absorber  is  that  of  the  optimal  steady-state 
response.  Under  the  optimal  steady-state  criterion,  the  optimal 
absorber  parameters  are  those  parameters  which  minimize  the 
maximum  response  of  the  structure  when  it  is  subjected  to  a  broad 
band  disturbance.  The  choice  of  absorber  parameters  alters  the 
transfer  funct  io:.  of  the  system.  This  solution  is  often  refered 
to  as  the  ’Classical  Tuned  Spring-Mass-Damper  Vibration  Absorber.’ 
This  solution  yields  a  set  of  passive  absorber  parameters  which 
equate  and  minimize  the  structure’s  two  peak  resonant  amplitudes. 
This  procedure  is  detailed  in  Reference  12  and  reviewed  here  for 
completeness.  A  closed  form  solution  from  a  control  point  of  view 
is  presented  in  Reference  9. 

Tiie  approach  used  in  the  classical  ’Vibration  Absorber’ 
problem  calls  for  modifying  the  transfer  function  which  relates 
the  disturbance  on  the  structure  (P  in  Fig.  2.2)  to  the  response 
of  the  structure  (Xi  in  Fig.  2.2)  through  selection  of  the  passive 
absorber  parameters.  Assuming  a  sinusoidal  disturbance  and 
motion,  the  transfer  function  from  P  to  Xi  (Eg.  2.2)  is  obtained 
from  the  equations  of  motion 

m1x1  +  exj  +  (kj+k2)xi  -  cx2  -  k2x2  =  P  (2.1a) 

m2x2  -i  cx2  -t  k2x2  -  exj  -  k2xj  =  0  (2.1b) 


I j,  I  =  (kj-w  m2)  /{[(k2-w  mi)(k2-w  m2)-w  m2k2]  + 

w2c2[ki-to2(ni+m2 )  ]2}  (2.2) 

The  transfer  function  can  be  nondimensionalized  as 

IV.  1 9.  .  99  9  9  9.  ..  9  9  9  .  9  JZ.,2 

-p—  -  {HU  7  rp  -o  )  J  / 1  1/30  7  “17  -ini'  “©  )J  + 

I  ST  | 

4//V[72(l+/5)-l]2}  (2.3) 

where  the  notation  of  Reference  12  has  been  used: 

dST  =  P/kj  STRUCTURE’S  STATIC  DEFLECTION  (2.4a) 

pQ  =  (ki/jni  ) 1/2  DECOUPLED  STRUCTURAL  FREQUENCY  (2.4b) 

P  =•  (k2/m2)1/2  DECOUPLED  ABSORBER  FREQUENCY  (2.4c) 

6  -  pa/pQ  ABSORBER/STRUCTURE  FREQUENCY 

RATIO  (2.4d) 

7  =  u  /p  NONDIMENS IONAL  FOURIER 

FREQUENCY  RATIO  (2.4e) 

P  =  m2/mi  ABSORBER/STRUCTURE  MASS  RATIO  (2.4f) 

U  -  c/(2m2p0)  NONDIMENS IONAL  DAMPING  (2.4g) 

S  =  s/p0  NONDIME NS IONAL  LAPLACE  OPERATOR  (2.4h) 

The  characteristic  equation  for  the  poles  of  Eq.  2.3  is  given  by 

s4  +  2p(l+0)s3  +  (1+62(1+U))s2  4  2ps  4  62  .  (2.5) 


where  the  two  natural  frequencies  of  the  structure/absorber  system 


in  the  absence  of  a  damper  are 


i \,i\  =  j(62(1+/5)+1)  +  j[(62(H*)  +  1)2-462)1/2  (2.6)  | 

I 

i 

i 

Figure  2.3  shows  the  magnitude  of  the  response  Xj  to  the 
disturbance  P  versus  the  disturbance  frequency  for  the  arbitrary 
case  where  the  mass  ratio  (0)  is  0.02  (i.e.  the  actuator  is  2%  of 
the  structure’s  modal  mass),  the  frequency  ratio  (6)  is  1.0  (i.e. 
the  actuator  in  isolation  has  the  same  frequency  as  the  structure 
in  isolation)  and  the  damper  strength  is  as  indicated. 

Tritre  are  three  nondimeriEioral  absorber  parar^te’-s  in  the 
design  problem:  0  the  mass  ratio,  6  the  frequency  ratio  and  u  the 
damper  value.  For  now,  it  will  be  assumed  that  0  is  fixed  and  the 
other  two  parameters  can  be  varied  to  achieve  the  best 
steady-state  performance.  Figure  2.3  includes  several  important 
features  of  a  passive  absorber  system  with  a  fixed  mass  ratio  0. 

First,  for  any  value  of  the  frequency  ratio  6  there  exist  two 
frequencies,  indicated  by  points  S  and  T,  where  the  response 
magnitude  is  independent  of  the  passive  damping  in  the  system. 

Second,  for  this  value  of  6,  there  are  two  limiting  cases  of  the 
transfer  function:  no  damper  (u  =  0)  and  large  damper  (ji  -»  •).  If 
the  damper  is  removed  (u  =  0),  there  are  two  unbounded  resonant 
peaks  in  the  transfer  function  at  the  frequencies  given  by  the 
nondimens ional  frequency  ratios  7j  and  yj  (Eq.  2.6).  The 
corresponding  transfer  function  for  this  undamped  case  is 


(2.7) 


In  the  other  limiting  case,  an  infinitely  strong  damper  (//-»•), 
there  is  a  transfer  function  with  one  unbounded  resonant  peak 


I5M  * 

- V-l 

(2.8) 
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Since  all  choices  of  damper  ((/)  between  the  limiting  cases  yield  a 
system  with  bounded  response  (Fig.  2.3),  it  is  clear  that  the 
absorber  components  can  be  selected  to  yield  a  system  with  a 
minimum,  maximum  system  response. 

The  criterion  for  optimal  steady-state  response  of  a  system 
with  a  passive  vibration  absorber  leads  to  a  minimax  solution 
(i.e.  minimizing  the  maximum  system  response).  The  magnitudes  of 
the  transfer  function,  at  the  two  resonances  (Eq.  2.3),  are 
equated  and  minimized.  As  an  example,  the  curve  labled  n  -  0.1  in 
Fig.  2.3  illustrates  the  situation  where,  for  a  given  a  mass  ratio 
(fl)  of  0.02  and  a  frequency  ratio  6  of  1.0,  system  damping  is 
maximized  because  resonant  magnitude  is  minimized.  However,  if 
there  is  a  free  choice  of  the  frequency  ratio  6,  further 
reductions  in  the  peek  resonant  amplitude  might  be  obtained.  The 
procedure  outlined  below  derives  this  minimax  response  based  on 
free  choice  of  y  and  6  for  a  fixed  fi. 

The  first  step  of  the  procedure  for  maximizing  the  system 
damping  is  to  adjust  the  spring  of  the  passive  damper  (adjusting 
the  frequency  ratio  6)  to  equate  the  magnitudes  of  the  transfer 
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function  at  points  S  and  T  (i.e.  at  the  points  common  to  all 

transfer  functions).  Once  equated,  damping  is  added  through  the 
adjustable  viscous  damper  (p)  until  the  peak  resonant  response 
corresponds  to  the  magnitude  of  the  transfer  function  indicated  by 
points  S  and  T.  At  this  magnitude,  the  maximum  response  is  at  a 
minimum  and  the  minimax  criterion  is  satisfied.  The  frequency 
ratio  is  adjusted  first  because  the  tuned  damper  value  p^^ 

depends  on  the  frequency  ratio.  The  relation  which  determines 

frequencies  which  are  common  to  the  two  transfer  functions 
(denoted  by  i  and  7  )  can  be  found  by  equating  the  transfer 

O  I 

functions  for  the  limiting  cases  of  p  =  0  and  p  -♦  00  (Eqs.  2.7  and 

2.8), 


/  -  2 


'1  ♦  fi2(  1+^)  j  2  4  _26^ 

2  ■*  a  2  + 


(2.9) 


This  relationship  can  be  solved  for  7,,  and  7^.  Then,  the 
objective  is  to  find  the  frequency  ratio  (6)  which  makes  the 
magnitudes  of  the  transfer  function  at  the  frequencies  7^  and  7^, 
equal.  By  evaluating  Eq.  2.0,  at  7^  and  7^,  and  equating,  we  get 


2  2 
7s  7t 


(2.10a) 


where  the  magnitudes  are  related  by 


(2.10b) 


Since  the  solution  of  Eq.  2.9  is  of  the  form 
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7  -  (7S  +  7t)7  +  7s7t  =  0 


(2.11) 


equating  the  i  coefficient  in  Kq.  2.9  to  the  right  hand  side  of 
Eq.  2. 10a  gives 


h  =  +  Aw) 


(2.12) 


which  can  be  reduced  to  yield  the  expression  for  the  optimal 
actuator-to-modal  frequency  ratio  (6^^), 


OPT  1  +  fi 


(2.13) 


The  frequencies  associated  with  points  S  and  T  are  found  by 
placing  Eq.  2.13  into  Eq.  2.9  to  get 


-M1 ;  ns  ] 


(2.14) 


The  magnitude  of  the  transfer  function  associated  with  points  S 


and  T  is  found  from  Eq.  2.8  as 


|X,  I  .  [[  2  +  fil 

'^st's,t  ' 


(2.15) 


This  gives  the  minimum,  maximum  magnitudes  of  the  transfer 


function  using  a  passively  tuned  vibration  absorber  for  a  fixed 


mass  ratio  0.  As  an  example,  for  a  system  with  a  mass  ratio  of  0 
-  0.02,  the  optimal  frequency  ratio  6,  from  Eq.  2.13,  is  0.98. 
The  plot  of  magnitude  of  response  versus  frequency  for  this  system 
is  shown  in  Fig.  2.4. 

Now  that  the  optimal  frequency  ratio  has  been  found,  the 

optimal  value  for  the  nondimens ionalized  damping  (n)  is  found  by 

2 

solving  Eq.  2.3  for  n  and  substituting  for  6  from  Eq.  2.13,  y 

O 

and  -»T  from  Eq.  2.14  and  using  X,/^ST  from  Eq.  2.15.  This  gives 
an  expres  in  the  general  form 


2 

"opt 


N  -  Q(  XJ/dST  )2 
P(  X,//1ST  )2  -  M 


(2.16a) 


where 
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(2.16b) 


(2.16c) 


(2. 16d) 
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(2. 16e) 


which  is  indeterminant  because  both  the  numerator  and  denominator 
vanish  at  both  7,,  and  7^.  Therefore,  the  optimal  nondimensional 
damper  value  is  found  by  taking  the  limit  as  the  frequency 
approaches  either  or  7^.  An  approximate  optimal  value  of  y  can 
actually  be  found  by  considering  a  value  of  7g  or  7^  slightly 

giver.  Lr.  2.,.  2.14.  tui  ihe  example  of  a 
mass  ratio  0  -  0.02,  the  optimal  frequency  ratio  is  found  from  Eq. 
2.13  as  6  =  0.98  and  the  optimal  nondimensional ized  damping  is 
found  from  Eq.  2.16a  as  y  -  0.086.  The  corresponding  response  of 
the  structural  mass  is  shown,  in  Fig.  2.4,  by  the  curve  labled 
with  this  optimal  value  of  y.  This  curve  represents  the  optimal 
steady-state  response,  based  on  the  criterion  of  minimizing  the 
maximum  value  of  the  transfer  function. 

Figure  2.5  and  Table  2.1  give  some  insight  into  this  tuning 
scheme.  Figure  2.5  is  a  plot  of  the  natural  modal  frequencies  of 
the  system  versus  the  frequency  ratio  (6)  for  a  mass  ratio  of  4  = 
0.02.  The  vertical  line  is  located  at  the  optimal  frequency  ratio 
(60pt),  as  calculated  from  Eq.  2.13.  As  seen  in  the  figure,  the 
optimal  solution  is  that  frequency  ratio  (6)  that  places  the 
natural  modal  frequencies  as  close  together  as  possible.  This 
allows  the  greatest  ber/st  interact  <  on  and  maximized 
motion  across  the  passive  damper. 

Table  2.1  lists  eigenvectors  which  show  the  motion  of  the  two 
masses  for  several  values  of  frequency  ratio  (6).  For  high  values 
of  6  (10.00),  relative  motion  (the  difference  between  the  motions 


of  the  masses)  is  almost  zero  in  the  low  mode.  For  low  values  of 


6  (0.10),  relative  motion  in  the  high  mode  is  reduced.  Note  that 
all  values  of  6  yield  some  relative  motion  and  corresponding 
damping  in  the  high  mode.  Since  overall  damper  effectiveness  is 
maximized  by  maximizing  the  relative  motion  between  the  two  masses 
on  average  in  all  modes,  the  best  choice  is  if  the  absorber 
frequency  is  made  close  to  but  lower  than  the  decoupled 
structure’s  modal  frequency.  In  this  way  the  relative  motion 
across  the  damper  in  the  two  modes  is  approximately  equal. 

Thus  far  an  algorithm  has  been  determined  for  finding  the 
optimal  frequency  ratio  (6)  and  optimal  nondimensional  damper  (p) 
given  a  fixed  mass  ratio.  No  criterion  has  been  set  for 
determining  the  best  mass  ratio  (/}).  As  seen  in  Figure  2.4,  the 
minimum  achievable  value  of  the  maximum  amplitude  is  set  by  the 
frequency  separation,  subject  to  the  constraint  that  the 
magnitudes  at  points  S  and  T  be  equal.  By  comparing  the  maximum 
response  for  larger  mass  ratios  (0  =  0.2  in  Fig.  2.6)  with  that 
for  the  smaller  mass  ratios  (0  =  0.02  in  Fig.  2.4)  it  is  clear 
that  a  larger  mass  ratio  is  advantageous.  Bq.  2.15  gives  the 
relationship  between  mass  ratio  and  minimum,  maximum  magnitude  of 
the  transfer  function.  While  steady-state  response  improves  with 
larger  mass  ratios,  overall  system  mass  also  increases.  If  an 


optimal  mass  ratio  (0)  is  desired,  the  choice  of  0  must  be 
determined  through  a  systems  analysis  involving  a  tradeoff  between 
total  system  mass  and  damper  performance. 


2.1.2  POLE  PLACEMENT  DESIGN  CRITERION 


A  different  insight  into  the  effectiveness  of  a  passive 
vibration  absorber  can  be  gained  by  examining  the  location  of  the 
poles  in  the  S-plane  of  the  one  DOF  structure  plus  absorber 
system.  This  S-plane  interpretation  gives  a  better  understanding 
oi  tne  transient  characteristics.  as  previously  stated,  the 
optimal  pole  placement  design  usually  calls  for  placing  the  poles 
as  far  to  the  left  as  possible  in  the  S-plane,  thereby  decreasing 
the  system  response  time. 

Figure  2.7  shows  the  upper  half  of  the  S-plane  representation 
of  the  poles  of  this  one  DOF  structure  plus  absorber  system, 
plotted  for  the  case  of  optimal  frequency  ratio  (^Qp^,)  (derived 
fro*"  Fo  2.1^  »c  the  ’’elue  of  ooodimensioral  damping  n  is 
increased.  As  the  damper  strength  increases  the  system  poles  move 
to  the  left,  pass  through  the  poles  corresponding  to  i<QpT» 
approach,  meet  (equal  damping  ratios)  and  split.  Point  A 
corresponds  to  the  value  of  u  which  yields  the  case  where  the 
poles  are  the  farthest  left.  This  value  of  n  is  designated 
^OPT-PP  ^or  °Pti,na^  pole  placement  damping. 

The  optimal  damping  value  (^Qp^-pp)  for  the  pole  placement 
deeign  can  he  deriv-v*  fly  -oeni  piloting  the  characteristic 
equation  (Eq.  2.5),  the  S-plane  coordinates  and  nondimens ional 
damping  of  the  optimal  pole  placement  solution  can  be  determined. 
At  point  A  of  Fig.  2.7,  the  roots  of  Eq.  2.5  must  be  double  pairs 
of  complex  conjugate  poles.  Substituting  the  optimal  frequency 
ratio  (fiQpT)  from  Eq.  2.13  into  Eq.  2.5,  and  equating  the  pairs  of 
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roots,  four  equations  arise  from  equating  coefficients 


2p(l  +  0)  =  -4a 


2  +  0  „  2  2 
- - £  =  6a  4  2b 

1  +  4 


2#/  =  -4a(a2  +  b2) 


1  =  fl2 
i  +  fi 


(2.17a) 


(2.17b) 


(2.17c) 


(2. 17d) 


where  the  S-plane  pole  locations  are  S  =  a  +  bi.  Since  2.17c  and 
2.17d  can  be  combined  to  form  2.17a,  there  are  three  independent 
equations  with  which  to  solve  for  the  real  part  (a),  the  imaginary 
part  (b)  of  the  pole  location  and  the  nondimensional  damping  (p) 
in  terms  of  the  mass  ratio  (0) .  Solving  Eqs.  2.17  yields 


2  0 

"°pi-pp  -  77777s 


(2.18a) 


a  I  =  -0.5  — ~  =  Re(S) 

'A  J  1  +  0 


(2.18b) 
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4  -  0 
1  +  0 


lm(S) 


(2.18c) 


For  the  same  example  computed  in  the  last  section,  that  of  a  mass 
ratio  of  0.02,  the  optimal  value  of  the  nondimensional  damper  is 


=  0.137. 


^OPT-PI’ 

In  this  derivation  of  ^0pT_pp>  we  have  assumed  that  Eq.  2.13 
provides  the  optimal  value  of  6  even  though  it  was  defined  based 
on  the  steady-state  criterion.  The  rationale  for  this  choice  of  6 
=  6qPT  is  that  this  is  the  exact  value  of  the  frequency  ratio 
which  causes  the  two  branches  of  the  root  locus  (Fig.  2.7)  to 
exactly  converge  at  point  A.  Therefore,  the  damping  of  each  of 
the  two  poles  will  be  identical,  and  maximized. 

As  in  the  steady-state  design  procedure,  guidance  is 
available  in  choosing  6  and  p,  but  not  the  mass  ratio  p.  Again  it 
is  clear  that  a  larger  mass  ratio  gives  a  more  damped  system 
(contrast  Eqs.  2.18b  and  2.15),  and  so  again,  a  systems  analysis 
must  be  performed  to  determine  the  optimal  mass  ratio. 

T+  ""y  in<5+. -”«-+”■'>  steady— state  responses 

and  transient  responses  o?  the  two  optima]  damper  cases.  Figure 
2.8a  compares  the  steady-state  responses  of  the  optimal 
steady-state  (A<Qpj)  and  pole  placement  (^QpT_pp)  designs  for  a 
mass  ratio  of  p  =  0.02.  Note  that,  as  would  be  expected,  the 
optimal  steady-state  design  exhibits  better  steady-state  behavior. 
This  is  because  the  peak  responses  of  the  two  modes  in  the  optimal 
pole  placement  design  overlay.  Figure  2.8b  shows  the  transient 
r  tv.  i nf- i T*-,n  solution  also 

exhibits  slightly  better  transient  behavior  due  to  the  destructive 
interference  of  the  two  modes,  which  appears  as  a  beating.  Though 
the  two  design  procedures  result  in  similar  vibration  absorber 
parameters,  and  the  two  solutions  exhibit  very  similar  behavior, 
the  performance  of  the  steady-state  design  is  slightly  preferable. 
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Thus  far,  tuning  rules  have  been  developed  for  the  first  two 
cases  of  vibration  absorber  optimization.  For  the  steady-state 
case,  the  following  steps  are  followed: 

1.  Choose  the  mass  ratio  (fi)  through  a  tradeoff  between  mass 

and  performance,  as  dictated  by  Eq.  2.15. 

2.  Determine  the  frequency  ratio  ( ^OPT >  us*n^  2.13. 

3.  Determine  nondimens ional  damping  (**0pj)  through  Eq.  2.16. 
For  the  pole  placement  case,  the  following  steps  are  followed: 

1.  The  value  of  the  mass  ratio  (fi)  is  determined  through  a 

tradeoff  between  mass  and  performance,  as  dictated  by 
Eq.  2.18b. 

2.  The  frequency  ratio  (6^^)  is  determined  by  Eq.  2.13. 

3.  The  nondimens ional  damping  (^0pj_pp)  *s  fQund  through  Eq. 


2.1.3  WEIGHTED  QUADRATIC  COST  DESIGN  CRITERION 


A  third  possible  design  criterion  involves  selecting  absorber 
parameters  which  cause  the  system  response  to  minimize  a  weighted 
quadratic  cost  function  of  the  states.  The  first  section  of 
Appendix  A  details  this  analysis,  which  yields  the  set  of  optimal 


redetermined 


absorber-to-structure  mass  ratio.  The  resulting  set  of  optimal 
absorber  parameters  closely  matches  those  derived  by  the  previous 
two  criteria.  Thus  all  three  cases,  though  approaching  the 
absorber  tuning  problem  from  three  distinct  definitions  of  optimal 
performance,  result  in  similar  absorber  designs.  In  fact, 
considering  the  accuracy  to  which  the  absorber  parameters  can  be 
physically  adjusted,  these  three  solutions  might,  be  considered 
identical . 

Up  to  this  point,  the  interaction  of  the  passive  vibration 
absorber  has  been  studied  with  a  plant  with  only  a  single  degree 
of  freedom.  Yet  realistic  structures  have  higher  order  dynamics. 
Therefore,  the  next  step  is  to  derive  tuning  rules  for  selecting 
absorber  parameters  which  best  introduce  damping  into  higher  order 
structural  models.  Due  to  the  similarity  between  the  three  design 
cases  presented  in  damping  the  single  degree  of  freedom  structure, 
only  the  steady-state  design  solution  will  be  extended  in  Section 
2.2  to  the  damping  of  multiple  DOF  structures.  By  introducing  and 
tuning  a  vibration  absorber  to  a  structural  model  containing  two 
modes,  rules  can  be  derived  which  allow  the  extrapolation  of  the 
solution  derived  in  this  section  to  the  passive  damping  of  higher 
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2.2  TUNING  A  SINGLE  DOF  PASSIVE  VIBRATION  ABSORBER  TO  A  TWO  DOF 

STRUCTURAL  MODEL 

The  two  degree  of  freedom  (DOF)  system,  analyzed  in  the 
previous  section,  represents  the  simplest  model  describing  the 
•internet  ion  hetwppn  a  vihr«tinn  absorber  and  a  structure.  A  more 
accurate  structural  model  would  include  many  modes,  but  this  model 
could  not  be  dealt  with  as  completely  as  the  low  order  model.  The 
simplest  model  which  includes  an  absorber  and  more  than  one 
structural  mode  is  a  three  DOF  model  in  which  the  structure  is 
represented  by  two  modes,  or  masses,  and  the  absorber  by  one. 
This  section  studies  tuning  of  a  vibration  absorber  in  a  three  DOF 
model.  With  this  information,  the  tuning  rules  for  achieving 
maximum  system  damping  can  be  extrapolated  to  higher  order,  more 
realistic  structural  models. 

As  shown  in  Figure  2.9,  the  two  mode  structure  has  been 
modeled  using  two  lumped  masses  (M0)  and  two  interconnecting 
springs  (K0).  In  this  section,  two  approaches  will  be  taken  in 
devising  tuning  rules  for  this  system.  In  the  first  approach,  the 
selection  of  the  absorber  spring  (Kj)  will  be  set  based  upon  the 
structure’s  lower  modal  mass  in  a  fashion  similar  to  that  which 
resulted  in  Eq.  2.13.  Then,  the  damper  (C)  will  be  adjusted  to 
determine  its  effect  on  both  the  lower  and  the  higher  structural 
modes.  Secondly,  the  absorber  will  be  tuned  to  the  higher  mode 
and  the  same  process  repeated  to  determine  the  effect  of  the 
damper  on  both  modes.  These  two  analyses  will  be  compared  and  the 


best  approach  determined.  Again,  the  control  force  ( F_, )  is 
ignored  and  the  analysis  is  performed  assuming  a  fixed 
absorber/structure  mass  ratio. 

The  procedure  presented  here  is  similar  to  that  followed  in 
section  2.1.1  for  steady-state  optimization.  Again,  sinusoidal 
disturbances  and  motion  are  assumed.  Nondimensionalizing  the 
equations  of  motion 

m0xj  -*  2k0Xj  -  k0x2  =  0  (2.19a) 

m0x2  +  cx3  +  (k0  +  kj)x2-  k0Xi  -  cx3  -  k3x3  =  P  (2.19h) 

m3x3  +  cx3  -  cx2  -*•  k3x3  -  k3x2  =  0  (2.19c) 


results  in  the  transfer  function  relating  the  disturbance  on  the 
second  structural  mass  (P)  to  the  motion  of  that  mass  (Fig.  2.9) 


2222222  2 
=  0. 25(7  -2;  [  /  {[6^- 

-i2{  l+62(3+2/8)  )+74(3+62(  1+/J)  )-76]2+ 
47V(l-72(3+24)+74(l+/S)]2} 


(2.20) 


where 


dST  =  2P/ke  STRUCTURE’S  STATIC  DEFLECTION  (2.21a) 

p0  =  (k0/m0)i/2  CHARACTERISTIC  STRUCTURAL 


FREQUENCY 


(2.21b) 


p3  -  (k,/Bj)  DECOUPLED  ABSORBER  FREQUENCY  (2.21c) 

6  =  Pj/Pc  ABSORBER/STRUCTURE  FREQUENCY 

RATIO  (2.21d) 

•t  =  w  /p0  NONDIMENS IONAL  FOURIER 

FREQUENCY  RATIO  (2.21e) 

£  =  Bj / Dq  ABSORBER/STRUCTURE  MASS  RATIO  (2.21f) 

v  -  c/(la3p0)  NOmJIMENSlONAL  DAMPING  (2.21g) 

z  =  x  Mst  NONDIMENS  IONAL  DISPLACEMENT’  (2.21h) 


The  choice  of  having  the  force  (P)  act  on  the  second  mass  is 
arbitrary  since  its  location  only  affects  the  zeros  of  the  system. 
The  system  poles  which  most  directly  affect  resonant  amplitudes 
are  solely  a  function  of  the  homogeneous  dynamics. 

Thp  "fry  |  1  r\L^  i  r»  rf  or*  c,  1  yij  4  c  the*  th-»  ’absorber  mode’ 

(i.e.  tha*  mode  which  exhibits  most  of  the  absorber  motion)  is 
placed  near,  in  frequency,  to  a  structural  mode,  the  interaction 
between  these  two  modes  in  the  vicinity  of  their  resonances  can  be 
treated  similarly  to  the  two  DOF  analysis  previously  presented. 
If  this  assumption  is  correct,  the  relations  previously  presented 
can  be  readily  applied  to  this  three  DOF  system.  Therefore,  the 
tuning  rules  developed  for  a  single  structural  mass  will  be 
•*r,p !  i  t  '  r hp  ?  ■  ? -if  ? h  > wo  of  freedom  and 

the  resulting  response  of  the  full  three  degree  of  freedom  (DOF) 
system  will  be  evaluated.  If  these  rules  provide  an  optimal 
steady-state  solution  for  the  three  DOF  model,  using  the  minimax 
criterion,  then  it  will  be  assumed  that  they  can  be  used  to  design 
passive  absorbers  for  higher  order  structural  systems. 
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In  the  two  DOF  problem,  the  passive  absorber  components  were 
tuned  to  the  inner  mass  (Mt  in  Fig.  2.2)  using  Eqs.  2.13  and  2.16. 
For  the  three  DOF  case,  the  absorber  will  be  tuned  to  the  modal 
mass,  first  of  the  lower  and  then  of  the  higher  structural  mode. 

The  first  step  is  to  determine  the  modal  characteristics  of 
the  structural  system  when  the  absorber  is  removed  or  absent. 
Eliminating  the  absorber  parameters  from  Eqs.  2.19b  and  2.19c 
gives  the  nondimens ional  equations  of  motion  for  the  two  DOF 
structural  model  as 


(2.22) 


The  corresponding  eigenvectors  and  frequencies  are 


f 


T 

i 


1 


[1  1.618]  *r  i  =  0.618 


(2.23a) 


T 


T 

2 


[1  -0.618]  i7  -  1.618 


(2.23b) 


The  modal  masses,  as  perceived  from  the  actuator  attachment  point, 
are  found  by  normalizing  Eqs.  2.23a  and  2.23b  to  achieve  a 
displacement  of  unity  for  the  outer  mass.  By  assembling  the 
eigenvector  matrix  and  performing  the  following  operation,  the 
modal  masses  are 


O' 

m°  5k 

.0  m0. 
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(2.24) 


fi. 3o2  i) 
l  0  3.618 

where  fi  is  the  ratio  of  nodal  mass  to  the  reference  Bass  (M. ). 

Now  that  the  nodal  frequencies  and  Basses  have  been 
deternined,  the  tuning  rules  developed  in  Section  2.1.1  can  be 
applied.  At  this  point,  it  is  convenient  to  define  two  new 
variables:  tile  ratio  oi  accuacor-co-noaai  frequency  (6,^)  and  the 

ratio  of  actuator-to-modal  nass  ( fi .  );  as 

♦n 

6  -  p  /i  (2.25a) 

♦n  an 

fi%n  n.3/Mn  -  fi/fin  (2.25b) 

Note  that  three  njacc  ratios  have  been  defined.  As  in  Section  2.1, 

0  is  defined  as  the  ratio  of  absorber  mass  to  the  nass  of  one 

structural  element  shown  in  Figure  2.9  (i.e.  fi  =  Mj/M. ).  fi  , 

°  n 

defined  in  Equation  2.24,  represents  the  ratio  between  the  modal 
masses  of  the  two  degree  of  freedom  structure  and  the  nass  of  a 
structural  element  (M0).  The  third  mass  ratio  fi.  is  the  indexed 
ratio  of  absorber  mass  to  the  nodal  mass  of  the  structure  (Eq. 
2.25b).  The  same  variable  suffixes  have  been  used  for  the 
frequency  ra  t  i.cs  , 

In  the  selection  of  the  optimal  absorber  parameters  these 
ratios,  with  respect  to  the  nodal  frequency  and  mass  (6#n  and 
/J^),  will  be  used  in  lieu  of  the  ratios  to  the  physical  frequency 
and  mass  used  in  Section  2.1.  As  in  the  two  DOF  analysis,  let  the 
optimal  6^  be  defined  with  respect  to  the  modal  mass  ratio 
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*n 


OPT  1  +  0 


*n 


(2.26) 


Then,  the  optimal  value  of  the  frequency  ratio  (6),  defined  in  Bq. 

2.21d,  is  related  to  6.  by 

*nOPT 


60PT 


^a^n  „ 

pn  p‘  “  *nOPT,n 


(2.27) 


As  a  numerical  example,  assume  that  the  mass  ratio  0  =  0.02 
for  both  the  low  and  high  modes.  This  can  not  be  the  case  if  the 
same  absorber  mass  is  used,  but  this  assumption  makes  this  problem 
similar  to  the  numerical  cases  presented  in  Section  2.1.  Table 
2.2  lists  the  values  for  the  various  nondimens ional  parameters 
used  in  this  example.  The  first  column  corresponds  to  the  tuning 
of  the  actuator  to  the  low  frequency  structural  mode,  and  the 
second  column  corresponds  to  the  higher  mode.  The  next  step  is  to 
verify  that  the  absorber  parameters  derived  by  application  of 
these  rules  achieve  a  nearly  optimal  steady-state  response  for  the 
full  three  DOF  system. 

Recall  the  absorber  frequency  can  be  tuned  to  either  the 
frequency  of  the  lower  or  the  higher  structural  mode.  Up  to  this 
point  the  rule  for  deriving  the  optimal  damper  value  has  not  been 
extended  to  the  case  of  the  multiple  DOF  structural  model  (Figure 
2.9).  Fig.  2.10  illustrates  the  frequency  transfer  function,  or 
the  response  to  broad  band  disturbance,  of  the  three  DOF  system 
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when  the  absorber  frequency  is  tuned  to  the  frequency  of  the  lower 
structural  mode  according  to  Equation  2.26.  Two  curves  are  shown 
in  Figure  2.10.  For  Curve  1,  the  damper  has  been  set  to  optimally 
damp  the  lower  mode,  using  the  rule  for  optimal  damper  selection 
given  by  Eq.  2.16  in  Section  2.1.1.  This  damper  value  equates  the 
amplitude  of  the  actuator  and  structural  response  of  the  lower 
Th i c  ’.'"lid" ter  the  — «■  ■; thx.t  th~  tunieg  lew  for  a  one 
DOF  structure  governing  absorber  frequency  and  damping  can  be 
applied  to  this  two  DOF  structure,  provided  reasonable  separation 
in  structural  frequencies  exists.  However,  further  examination  of 
Curve  1  of  Fig.  2.10  reveals  that  while  the  two  resonant  peaks  in 
the  vicinity  of  the  lower  mode  have  been  equalized,  the  peak  of 
the  higher  frequency  structural  mode  is  four  times  greater  in 

magnitude.  Since  the  addition  of  more  damping  will  decrease  the 

response  of  this  higher  mode,  the  generalized  minimax  criterion 

has  not  yet  been  satisfied. 

Curve  2  in  Fig.  2.10  shows  the  transfer  function  for  the 

system  with  approximately  three  times  more  damper  strength  than 
that  given  in  Curve  1.  Notice  that  the  maximum  response  of  the 
upper  mode  has  decreased  at  the  cost  of  increasing  the  response  of 
the  lower  mode.  However,  the  minimax  criterion  has  been  satisfied 
in  that  any  change  in  '(a mo®1'  i'rora  this  value  will  increase 

the  maximum  response.  Thus,  in  order  to  achieve  minimum  system 
response,  the  damper  value  must  be  increased  from  the  optimal  two 
DOF  value  (Curve  1)  to  the  value  that  equates  all  of  the  maximum 
response  magnitudes  (Curve  2).  This  passive  absorber 

configuration  gives  the  optimal  steady-state  solution  for  the 
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three  DOF  system,  in  the  case  in  which  the  absorber  has  been  tuned 
to  the  lower  mode. 

A  similar  analysis  can  be  performed  in  which  the  objective  is 
to  tune  the  absorber  to  the  higher  structural  mode  using  the 
frequency  ratio  derived  from  the  simple  two  DOF  approach.  Curve  1 
of  Figure  2.11a,  b  and  c  shows  the  frequency  response  for  the  case 
in  which  the  absorber  damper  has  been  optimally  adjusted  for  the 
higher  structural  mode  only.  Note  that  while  the  high  mode  has  a 
peak  magnitude  on  the  order  of  unity,  the  low  mode  peak  magnitude 
exceeds  1000.  Curve  2  shows  the  minimum  low  mode  magnitude 
achievable  through  damper  adjustment.  The  maximum  amplitude  of 
the  low  mode  has  been  reduced  to  about  500,  at  the  cost  of 
increasing  the  r.uximum  response  of  the  high  mode  to  15.  The 
significant  result  of  this  case  is  that  no  damper  adjustment  can 
be  made  which  equates  all  of  the  peak  magnitudes.  Therefore, 
structural  modes  lower,  in  frequency,  than  the  mode  to  which  the 
absorber  is  tuned  exhibit  relatively  small  levels  of  damping. 
This  result  can  be  explained  using  the  modal  characteristics  of 
the  system. 

An  eigenvector  argument,  which  is  valid  for  small  damping, 
can  be  used  to  illustrate  the  lack  of  damper  effectiveness  in 
modes  lower  in  frequency  than  the  mode  to  which  the  absorber  is 
tuned.  Figure  2.12  shows  the  directions  of  motion  for  the 
possible  modes  of  the  system.  Table  2.3  lists,  eigenvectors  for 
various  frequency  ratios  (6)  with  the  motion  of  the  absorber  mass 
normalized  to  unity.  Together,  this  figure  and  table  illustrate 


how  the  absorber  and  structure  interact. 


Several  important  features  can  be  seen  in  Fig.  2.12.  The 
first  column,  containing  models  of  two  DOF,  lumped  mass  systems, 
illustrates  the  relative  directions  of  motion  between  the  two 
masses  for  the  two  structural  modes  when  the  absorber  has  been 
removed.  Columns  2  and  3  represent  the  three  modes  which  exist 
with  the  absorber  in  place.  The  absorber  in  Column  2  is  modeled 
ui  fK  a  frofinon^w  j'is t  be! o'.;  ihct  of  the  1  1  structural  mode  (a 
more  flexible  absorber;  whereas  column  three  shows  the 
eigenvectors  for  the  case  where  the  absorber  is  tuned  to  the 
higher  structural  mode  (a  stiff  absorber).  Since  the  modal 
damping  is  dependent  on  the  relative  motion  between  the  outer  two 
masses,  a  flexible  absorber  guarentees  significant  damping  in 
system  modes  2  and  3  due  to  opposing  motion  between  the  outer  two 
masses  in  these  modes.  The  system  response  corresponding  to  a 
stiff  absorber  only  exhibits  opposing  motion  between  the  outer  two 
masses  in  the  highest  mode.  In  fact,  the  stiffer  the  absorber 
(larger  6),  the  less  the  relative  motion  and  damping  in  the  lower 
two  modes. 

Table  2.3  verifies  this  argument.  This  table  lists  the  modal 
frequencies  and  eigenvectors  for  various  values  of  the  frequency 
ratio  6.  The  eigenvector  of  Mass  3  corresponds  to  the  absorber’s 
motion  The  lowest  ■node  e*.b.< hits  its  inrg*?'.  relative  motion  for 
small  values  of  frequency  ratio  (6).  As  the  frequency  ratio  is 
increased,  relative  motion  between  the  outer  two  masses  in  the 
lowest  mode  decreases  and  damper  effectiveness  is  lost.  Damper 
effectiveness  is  also  significantly  reduced  in  the  second  mode  at 
high  frequency  ratios.  The  third  mode  always  exhibits  opposing 


motion  between  the  outer  two  masses  and  as  a  result,  damped 
behavior.  Therefore,  to  most  effectively  introduce  damping  into  a 
structural  system  using  a  passive  vibration  absorber,  the  absorber 
should  be  tuned  to  the  lowest  mode  of  interest. 

When  designing  an  absorber  to  damp  several  modes  of  a 
multiple  DOF  system,  tune  the  absorber  frequency,  using  Eqs.  2.26 
and  2.27,  to  the  lowest  critical  mode.  This  might  be  the  lowest 
system  mode,  the  lowest  uncontrolled  mode,  or  the  lowest  frequency 
mode  which  can  be  most  affected  at  a  particular  location.  When 
adjusting  the  absorber  damper,  only  the  tuned  mode,  and  those  of 
higher  frequency,  can  be  significantly  damped.  In  other  words, 
when  attempting  to  passively  damp  several  modes,  tune  to  the 
lowest  and  adjust  the  damper  for  the  desired  response  of  it  and 
the  higher  modes. 

The  results  of  these  first  two  sections  can  be  summarized  as 
follows.  Thus  far,  three  passive  tuning  rules  have  been  derived 
for  the  two  DOF  model.  First,  Eqs.  2.13  and  2.16  size  the 
absorber  spring  and  damper  for  optimal  steady-state  response  based 
on  a  minimax  criterion.  For  the  optimal  pole  placement  response, 
Eqs.  2.13,  2.18a,  and  2.18b  are  used.  Thirdly,  a  weighted 
quadratic  cost  of  state  is  minimized  by  the  procedure  discussed  in 
Appendix  A.  The  three  DOF  problem  allows  the  extrapolation  of  the 
steady-state  tuning  rules  to  higher  order  models.  The  next  phase 
of  investigation  involves  the  introduction  of  a  control  force  into 
the  absorber,  transforming  it  from  a  passive  absorber  into  a 


control  actuator  with  passive  dynamics. 


2.3  SINGLE  STRUCTURAL  DOF  WITH  A  SINGLE  DOF  ACTIVE  CONTROL 

ACTUATOR 

Now  that  the  rules  for  tuning  a  passive  absorber  to  one  and 
two  DOF  structural  models  have  been  derived,  rules  for  tuning 
actuators,  whose  oassive  comonnen+s  reserohl**  wihmt.ion  absorbers, 
wiii  be  derived.  The  purpose  of  this  actuator  is  to  provide 
efficient  active  control  of  the  original  one  DOF  structure  (Fig. 
2.2)  using  state-dependent  feedback  (i.e.  F  included).  In  this 
section  the  passive  actuator  design  best  suited  for  an  active 
feedback  controller  is  determined.  It  is  assumed  that  the  control 
system  designer  has  the  freedom  to  select  the  passive  actuator 
characteristics,  and  therefore  the  actuator  components  can  be 
selected  in  such  a  way  as  to  best  make  use  of  the  control  effort 
available  to  achieve  a  certain  system  performance.  In  designing 
the  passive  aspects  of  an  actuator,  two  approaches  are  possible, 
the  passive  components  can  be  designed  first,  then  the  active 
gains  determined  (sequential  design);  or  the  passive  components 
and  active  gains  can  be  determined  simultaneously.  The  matrix 
presented  in  Table  2.4  summarizes  the  design  cases  analyzed  in 
this  chapter.  The  first  row  of  cases,  consisting  of  the  design  of 
passive  vibration  absorbers  using  the  steady-state,  pole  placement 
and  minimum  quadratic  cost  formulations ,  was  presented  in  Section 
2.1.  The  middle  row  includes  cases  of  sequential  designs  in  which 
the  passive  actuator  characteristics  have  been  designed  using  one 
of  the  techniques  in  Section  2.1  and  optimal  active  feedback  gains 
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were  derived  using  the  Linear  Quadratic  Regulator  (LQR) 
formulation.  The  third  row  involves  the  simultaneous  optimization 
of  the  passive  actuator  components  and  active  feedback  gains.  The 
objectives  of  this  section  are  threefold:  to  determine  the  passive 
components  of  the  control  actuator  that  maximizes  the  benefits  of 
the  available  control  effort;  to  determine  whether  the  passive  and 
active  portions  can  be  derived  sequentially  rather  than 
simultaneously;  and  determine  if  the  passive  components  of  the 
optimal  actuator  resemble  the  components  of  the  optimal  passive 
vibration  absorber. 

It  should  first  be  determined  whether  or  not,  for  this  simple 
structure,  active  control  provides  any  performance  improvements 
not  achievable  through  passive  means.  A  spring/damper  is 
effective  because  it  exerts  equal  and  opposite  forces  on  the  two 
masses,  between  which  it  is  attached,  proportional  to  their 
relative  motion.  For  this  reason,  the  stiffness  and  damping 
matrices  are  symmetric  in  the  passive  case.  This  corresponds  to 
self-adjoint  system  matrices.  Self-adjointness  does  not 
necessarily  exist  in  the  active  case.  Because  of  this,  active 
deunping  provides  more  versatility  in  dynamic  design  in  that 
nonsymmetric  state  feedback  terms  can  be  added  to  the  system 
matrices.  From  another  viewpoint,  it  can  be  shown  that  the  linear 
time-invariant  system,  shown  in  Fig.  2.2,  is  fully  controllable. 
This  means  that  the  system  poles  can  be  arbitrarily  placed  through 
full  state  feedback  gain  selection.  Arbitrary  pole  placement  is 
not  achievable  using  only  passive  components  capable  of 
functioning  in  the  space  environment.  The  procedure  followed  in 
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this  section  determines  the  effects  and  benefits  of  these  feedback 
terms. 

Feedback  is  introduced  through  the  control  force  labeled  F 

is 

in  Figure  2.2.  In  this  analysis,  the  control  force  F^  is 
determined  through  full  state  feedback 


1*1 


2*2 


x  > 


(2.28) 


Since  this  force  is  motion-dependent,  the  action  of  the  control 
force  produces  homogeneous  terms  which  can  be  moved  to  the  left 
hand  side  of  the  governing  equations  of  motion 
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+ 
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-k2+Fx  k2+F2 

(2.29) 


Notice  the  nonsymmetric  contributions  of  the  feedback.  The  new 
nondimensional  transfer  function  is  given  by 
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(2.30) 


where,  in  addition  to  the  terms  defined  in  Eqs.  2.4a  through  2.4g, 
there  are 


52 


STRUCTURE’S  RATE  FEEDBACK  GAIN 


(2.31a) 

(2.31b) 


n  =  Fg/(2m2p0 ) 

p  =  F4/(2m2p0)  ABSORBER’S  RATE  FEEDBACK  GAIN 
a1 2-  F1/(m2p^)  STRUCTURE’S  POSITION  FEEDBACK 

GAIN  (2.31c) 

C2=  F2/(m2pf)  ABSORBER’S  POSITION  FEEDBACK 

GAIN  (2. 31d) 

These  relations  are  used  in  the  subsequent  control  analysis. 

The  procedure  followed  in  Appendix  A  uses  LOR  design 
techniques  first  to  optimize  the  passive  absorber  design  (as 
already  discussed  in  Section  2.1.3)  and  secondly  to  optimize  the 
active  feedback  gains  and  passive  actuator  components  both 
sequentially  and  simultaneously.  These  passive  actuator 

components  are  then  compared  to  the  passive  components  derived  for 
the  vibration  absorber  in  Sections  2.1.1,  2.1.2  and  2.1.3. 

Appendix  A  and  Reference  10  give  the  theory  and  analysis 
involved  in  the  LQR  technique.  In  the  simultaneous  formulation, 
the  passive  actuator  parameters  and  active  feedback  gains  are 
simultaneously  optimized  based  on  the  minimization  of  a  consistent 
quadratic  state  plus  control  effort  cost  function.  The  following 
is  a  list  of  the  main  conclusions  drawn  in  Appendix  A. 

1.  This  formulation  yields  an  antisymmetric  active 
control  gain  matrix.  There  is  no  symmetric  part 
in  the  gain  matrix  since  all  symmetric  terms  are 
realized  through  the  passive  parameters  whose 
existence  is  not  penalized  in  the  cost  function. 

2.  The  passive  components  in  the  simultaneous 
design  are  similar  to  the  components  of  the 
passive  vibration  absorber  designs.  The 
sequential  and  simultaneous  design  of  the 
passive  components  and  active  feedback  gains 


give  tne  same  design. 

3.  The  optimal  active  feedback  gains  may  be 

determined  assuming  that  the  passive  components 
of  the  actuator  are  the  same  as  the  optimal 
passive  absorber  and  these  gains  may  be 
superimposed  upon  the  passive  vibration  absorber 
design  to  create  the  optimal  active  control 
actuator.  In  other  words,  the  simultaneous 
optimization  of  passive  components  and  active 
gains  is  not  necessary. 

Figure  2.13  compares  the  responses  of  the  system  in  Fig.  2.2 
subjected  breed  rend  disturbance  fer  c  variety  of  design 

cases.  A  mass  ratio  of  0  -  0.1  was  used  for  comparison  with  the 
results  of  Appendix  A  and  Reference  14.  The  design  cases 
corresponding  to  the  response  curves  in  this  figure  are  listed 
below: 


A  -  Optimal  passive  steady-state  response  from 

Section  2.1.1.  (Row  1,  column  1  in  Table  2.4.) 

2  OpwiL'al  p  1  aCcuicii l  response  from  Section 

2.1.2.  (Row  1,  column  2  in  Table  2,4.) 

C  -  Response  o*  passive  system  optimized  using  a 
weighted  quadratic  cost  function  of  the  states 
(Section  2.1.3).  (Row  1,  column  3  in  Table 
2.4.) 

Dl-  Response  of  passive  design  from  A  with  the 

superposition  of  optimal  active  control  based 
on  a  minimum  quadratic  cost  of  displacement, 
velocity  and  active  control  effort.  (Row  2, 
column  1  in  Table  2.4.) 

D2-  Response  of  passive  design  from  A  with  the 
superposition  of  optimal  active  control  based 
on  a  minimum  quadratic  cost  of  displacement  and 
active  control  effort  (i.e.  same  system  as  in 
case  Dl  but  with  the  velocity  penalty  removed. 
(Row  2,  column  l  in  Table  2,4.) 

El-  Response  of  the  passive  component  of  the 

simultaneous  optimization  detailed  in  Appendix 
A.  The  active  feedback  terms  have  been 
excluded.  (Row  3,  column  3  in  Table  2.4.) 

E2-  Response  of  system  that  has  undergone 

simultaneous  optimization  using  displacement, 
velocity  and  control  effort  penalties  in  the 
quadratic  cost  function.  Same  as  curve  El  but 
with  active  feedback  included.  (Row  3,  column 
3  in  Table  2.4. ) 


As  seen  in  curves  A,  B,  C  and  El  of  Fig.  2.13,  the  passive 
absorber  and  passive  component  of  the  active  actuator  are  roughly 
equivalent.  In  addition,  whether  the  passive  actuator  parameters 
and  active  feedback  gains  are  optimized  simultaneously  (Curve  E2) 
or  optimal  feedback  gains  are  superimposed  on  the  passive 
steady-state  design  (Curves  D1  and  D2),  the  closed-loop  responses 
are  roughly  equivalent.  Therefore,  the  optimal  active  feedback 
gains  and  optimal  passive  actuator  components  can  be  superimposed 
to  yield  the  optimal  control  actuator.  This  allows  the  passive 
actuator  components  to  be  selected  using  the  rules  derived  in  the 
vibration  absorber  analysis  (Sections  2.1  and  2.2).  It  should  be 
noted  that  these  conclusions  are  only  rigorously  valid  for 
collocated  feedback  control  using  the  actuator  states  and  the 
collocated  (actuator  mounting  point)  beam  states. 

It  must  be  emphasized  that  the  vibration  absorber  and 
simultaneous  LQR  techniques  differ  in  their  definitions  of  optimal 
and  in  their  design  techniques.  Yet,  the  behavior  of  the 
resulting  systems  are  roughly  the  same  for  a  mass  ratio  of  0.1, 
which  represents  a  practical  case  in  which  the  mass  of  the 
actuator  is  only  a  small  fraction  (10\)  of  the  total  system  mass. 

Additional  insight  can  be  gained  into  the  behavior  of  the  LOR 
design  formulation  by  analyzing  the  root  locus  of  the  optimized 
system  pole  locations  as  the  mass  ratio  (fi)  is  varied.  Figure 
2.14  shows  the  root  locus  for  the  design  cases  listed.  The 
optimal  steady-state  vibration  absorber  poles  (Section  2.1.1,  Case 
A)  are  in  the  proximity  of  the  LQR  passive  design  poles  (Section 


2.1.3,  Case  C)  for  low  mass  ratios  (0  =  0.01  to  0  =  0.1).  But, 
for  high  mass  ratios,  the  LQR  passive  poles  approach  the  poles  of 
the  system  whose  active  and  passive  portions  were  designed 
simultaneously  (Case  E  for  0  >  0.75).  This  implies  that  for  a 
given  amount  of  control  effort  available  the  active  aspects  of  the 
actuator  become  less  important  at  high  mass  ratios  where  the 
damping  achievable  through  passive  tuning  is  significant. 

2.4  CONCLUSIONS 


The  original  goals  of  this  investigation  have  been  achieved 
and  they  are  summarized  here  along  with  the  tuning  rules.  Three 
definitions  of  optimal  passive  design  were  derived  and  their 
resulting  solutions,  compared.  The  steady-state  solution 
equalizes  and  minimizes  peak  magnitudes  of  the  modal  response. 
The  pole  placement  solution  places  the  system  poles  as  far  left  in 
the  S-plane  as  possible.  The  weighted  quadratic  cost  solution 
minimizes  the  cost  of  the  system  subjected  to  a  set  of  quadratic 
state  penalties.  All  three  techniques  resulted  in  almost 
identical  vibration  absorber  configurations. 

The  optimal  passive,  steady-state  solution  was  derived  for 
two  and  three  DOF  systems.  The  solution  methods  were  found  to  be 
similar  for  each  case.  The  method  calls  for  passively  tuning  to 
the  lowest  of  a  set  of  'critical'  modes  and  adjusting  the  damper 
setting  for  the  desired  response.  'Critical'  modes  are  defined  as 
the  lowest  modes,  those  outside  the  bandwidth  of  the  active 


controller  or  those  most  effectively  damped  at  that  absorber’s 
location. 

The  steady-state  passive  design  and  the  passive  component  of 
the  LQR  design  were  compared  and  were  in  rough  agreement  for  small 
values  of  mass  ratio  (/5) .  Therefore,  the  LQR  active  feedback 
gains  are  superimposable  with  the  vibration  absorber  design 
allowing  the  passive  and  active  designs  to  be  performed 
sequentially:  passive  followed  by  active  design.  A  practical 
advantage  is  realized  from  the  superposition  of  the  optimal  active 
and  passive  designs.  In  the  event  of  a  failure  of  an  active 
component,  the  structure  still  exhibits  maximum  passively  damped 
behavior. 

The  typical  design  steps  for  an  inertial  control  actuator  are 
listed  below.  An  actuator  model,  as  seen  in  Fig.  2.1,  is  assumed. 


1.  Select  the  actuator-to-modal  mass  ratio  using  a 
tradeoff  between  mass  and  performance. 

A.  For  steady-state  response,  Bq.  2.15 
determines  performance. 

B.  For  pole  placement,  Eq.  2.18b  determines 
performance. 

2.  Determine  actuator  stiffness  using  Eq.  2.26. 

3.  Determine  actuator  damping  using  either  Eq.  2.16 
or  Eq.  2.18a.  The  optimal  damper  value  will 
exceed  these  values  in  the  multimode  case. 

4.  Determine  a  suitable  quadratic  cost  index  for 
the  required  system  performance  and  calculate 
the  optimal  LQR  feedback  gains  (Appendix  A). 


3.0  CONTROL  PROCESSOR  ARCHITECTURE 


The  processor  is  that  control  system  component  that  is  used 
to  implement  the  closed-loop  compensation.  These  compensation 
alogrithms  can  be  realized  through  analog  circuitry  or  digital 
hardware  and  software,  or  a  hybrid  of  the  two.  Though  analog 
computers  are  able  to  track  and  process  higher  frequency  signals, 
without  discretization  effects  such  as  aliasing,  analog 
compensation  is  not  as  versatile  as  digital  compensation. 
Algorithms  can  be  more  easily  altered  if  implemented  in  software, 
and  in  addition,  digital  compensation  reduces  the  problems  of 
component  drift  and  noise  which  are  common  among  analog 
components.  digital  realization  also  enables  compensation  of  low 
frequency  signals  for  which  analog  circuits  often  contain 
unreliable,  high  drift  elements.  For  these  reasons,  a  digital 
control  implementation  was  chosen. 

Controller  bandwidth  is  the  limitation  encountered  when  using 
a  digital  computer  to  implement  control  algorithms.  Controller 
bandwidth  is  directly  related  to  controller  sample  rate,  which  is 
defined  as  the  number  of  times  the  control  algorithm  is  repeated 
per  second.  The  higher  the  control  rate,  the  larger  the 
achievable  bandwidth  since  the  Nyquist  frequency,  one  half  the 
sample  frequency,  effectively  marks  the  high  frequency  limit  of 
the  controller. 

Structures  are  dynamic  systems  of  almost  infinite  order  and 


the  objective  of  this  research  is  to  devise  high  authority,  high 


bandwidth  controllers  for  these  structures.  As  the  bandwidth  of  a 


structural  compensator  is  increased,  not  only  must  the  cycle  time 
of  the  fixed  number  of  calculations  in  the  digital  processor  be 
reduced,  but  the  order  of  the  estimator  model  (and  therefore  the 
number  of  calculations  per  cycle)  must  be  increased  to  reflect  the 
increase  in  bandwidth.  Due  to  the  limited  speed  of  a  given 
processor  and  the  virtually  unlimited  number  of  structural  modes 
some  method  is  needed  to  increase  the  effective  control  rate. 

Several  options  are  viable  for  increasing  the  controller 
bandwidth.  First,  as  advances  are  made  in  digital  technology, 
faster  processors  can  be  acquired.  Secondly,  parallel  processing 
can  be  used  to  perform  simultaneous  operations  in  parallel 
hardware.  Finally,  distribution  of  the  processing  task  among 
several  processors  enables  faster  processing  of  system  states 
which  increases  the  bandwidth  of  the  controller.  If  this 
distribution  of  tasks  is  done  in  a  hierarchic  sense  to  reflect 
structural  behavior,  further  benefits  may  occur. 

This  chapter  discusses  the  advantages  of  applying  distributed 
mul t iprocessor  control  to  modify  elastic  structural  behavior.  As 
advances  are  made  in  the  first  two  technological  areas,  faster  and 
more  parallel  processors,  they  can  be  readily  incorporated  into 
this  distributed  hardware  format  which  represents  a  long  term 
processing  solution.  Section  3.1  presents  a  proposed 
microprocessor  architecture  based  upon  distributed  processing  and 
control  hierarchy.  Section  3.2  details  two  digital  hardware 
configurations  which  support  this  architectural  format.  The 
present,  first  generation  processor  system  and  planned  expansion 


to  the  second  generation  system  are  detailed  in  Sections  3.3  and 
3.4,  respectively. 

3.1  ARCHITECTURE  FOR  A  DISTRIBUTED  HIERARCHY  OF  CONTROL 

The  processor  architecture  discussed  in  this  section  centers 
around  two  important  principles  which  potentially  enhance  the 
performance  of  high  speed  digital  control.  These  principles  are 
distributed  processing  and  hierarchy  of  control.  Distributed 
processing  allows  simple  multiplication  of  data  throughput 
capacity.  Since  the  processing  at  all  localities  is  performed 
simultaneously,  the  achievable  control  rates  are  significantly 
faster  man  tnac  aemevable  with  control  using  a  single, 
equivalent  processor.  This  increases  both  the  bandwidth  of  the 
controller  and  the  number  of  structural  modes  that  can  be 
controlled.  But,  there  are  many  options  for  distribution  of 
processing.  Parallel  processing  is,  in  a  sense,  a  scheme  of 
distributed  processing,  as  are  dynamic  task  allocation  schemes. 
One  potentially  useful  distribution  scheme  is  coupled  to  the 
concept  of  hierarchic  control. 

Control  hierarchy  can  be  used  to  reduce  the  order  of  the 
structural  model  required  in  each  processor.  In  such  a  scheme, 
each  processor  is  responsible  for  controlling  a  particular  set  of 
modes  or  range  of  wavelengths.  The  command  signal  sent  to  any 
given  actuator  is  the  sum  of  command  signals  from  the  set  of 
processors  whose  combined  model  fully  represents  the  controlled 
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modes  of  the  structure.  Figure  3.1  illustrates  the  processor 


network  for  distributed,  hierarchic  control.  Distributed 
processing  increases  control  speed  by  reducing  the  number  of 
sensors  and  actuators  for  which  each  of  the  lowest  level 
processors  is  responsible.  Hierarchic  processing  reduces  the 
number  of  modes  or  range  of  wavelengths  for  which  each  processor 
is  responsible.  Both  principles  tend  to  distribute  the  processing 
task,  increase  the  control  rate,  and  expand  the  controller’s 
bandwidth. 

In  the  following  discussion,  several  terms  are  used  which 
require  definition.  The  terms  ’master’  and  ’slave’  refer  to  any 
pair  of  processors  which  interact  directly  with  each  other.  The 
master  processor  maintains  some  level  of  authority  or  control  over 
its  respective  slave(s)  in  the  areas  of  control  synchronization 
and  data  transfer.  The  terms  ’local’,  ’ component ’  and  ’global’ 
refer  to  the  levels  of  hierarchic  processing  illustrated  in  Fig. 
3.1.  The  local  processors  interact  most  closely  with  the 
structure  through  sensor  data  retrieval  and  processing  and  the 
issuance  of  actuator  commands.  A  component  level  processor 
interacts  with  several  local  level  processors  (its  slaves).  A 
global  processor  is  in  the  highest  level  and  interacts  with 
several  component  level  processors  in  a  master/slave  arrangement. 
Interaction  between  a  component  and  local  level  processor  is  a 
second  example  of  master/slave  interaction. 


With  each  processor  controlling  a  certain  set  of  modes  at  a 
particular  location,  the  idea  of  hierarchic  control  naturally 
follows.  Assign  to  the  zero  or  local  level  processor  (Fig.  3.1) 


the  highest  frequency  modes  or  waves  which  lie  within  its 
bandwidth.  This  processing  level  has  the  highest  bandwidth 
because  most  of  the  processing  load  is  distributed  throughout  the 
higher  levels.  Space  these  processors,  and  their  corresponding 
actuators,  along  the  length  of  the  structure  as  is  appropriate. 
This  processing  level  would  typically  perform  local,  collocated 
state  feedback  using  high  speed  digital  hardware.  In  this 
situation,  each  processor  obtains  information  about  a  small 
portion  of  the  structure  and  controls  these  high  frequency,  short 
wavelength  modes  over  a  limited  section  of  the  structural  domain. 

At  this  point  a  one  or  component  level  processor  is  assigned 
to  preside  over  several  adjacent  zero  level  units  which  are  its 
slave  processors.  By  obtaining  all  of  the  sensor  data  from  its 
distributed  <m. o  level  units,  the  one  level  processor  can  observe: 
longer  wavelength  vibrations  better  than  its  associated  zero  level 
units  because*  it  possesses  sensor  information  over  a  large:r 
portion  of  the  structure.  This  enables  the  component  level 
processor  to  perform  feedback  between  multiple  sensors  and 
actuators  at  different  locations  (noncol located  feedback).  This 
processor  task  versus  level  is  illustrated  in  Fig.  3.2. 

If  the  lower  frequency,  longer  wavelength  modes  are 
considered  the  most  critical  modes  to  control  (i.e.  most  sensitive 
to  a  performance  index  and  most  easy  to  excite),  the  advantage  of 
this  distributed  arrangement  is  obvious.  The  component  level  can 
process  more  information  on  the  lower  frequency,  more  critical 
modes  at  a  slower  speed  since  these  lower  modes  lie  within  the 
bandwidth  of  processors  with  slower  control  rates.  If  need  be, 


more  component  levels  may  be  added  in  a  pyramid  fashion.  The 
highest,  global  level  of  processing  would  perform  the  highest 
level  of  control  on  the  lowest  structural  modes  and  have 
coprocessors  to  perform  such  tasks  as  component  failure  detection, 
control  reconfiguration  and  system  identification.  Figures  3.3a 
and  3.3b  show  two  possible  control  authority  arrangements  in  the 
frequency  domain  for  these  three  levels  of  processing. 

This  scheme  provides  several  advantages  over  single  processor 
control.  First,  portions  of  the  processing  task  can  be  performed 
simultaneously  allowing  the  control  of  higner  freqency  modes  due 
to  the  increase  in  controller  bandwidth.  Secondly,  control 
activation  can  parallel  the  construction  phase  in  a  building  block 
fashion.  As  the  structure  becomes  larger,  additional  processors 
ai  levels  may  be  added.  Thirdly,  local  level  processors 
manipulate  a  small  amount  of  information  on  a  small  portion  of  the 
structure  enabling  them  to  work  at  high  speeds.  Higher  level 
units  process  more  information  on  the  more  critical  modes  without 
reducing  the  bandwidth  of  the  entire  control  system. 

There  remain  a  number  of  important  issues:  The  applicability; 
appropriateness  or  effectiveness  of  such  a  concept;  the  design  and 
analysis  scheme  for  such  a  distributed  hierarchy;  and  its 
implementation  in  hardware.  The  former  two  issues  are  the  subject 
of  Reference  14.  The  following  sections  discuss  how  this 
architecture  can  be  implemented  in  hardware  and  Appendix  R 
discusses  how  structural  models  influence  the  distribution  of 


control  tasks  in  a  hierarchic  controller. 


3.2  COMPATIBLE  HARDWARE  DESIGN  OPTIONS 


Several  hardware  configurations  can  be  used  to  interface 
processors  in  a  distributed  and  hierarchic  architecture.  Two 
configurations  will  be  discussed  here:  a  single  bus  and  a  multibus 
system.  A  single  bus  system  could  support  all  of  the  control 
processor"  or  the  same  bus.  A  ’housekeeping*  processor  would  be 
used  to  synchronize  the  data  transfer  between  control  processors 
on  adjacent  levels.  In  Figure  3.4,  data  transfer  between 
processors  G  and  Cl,  or  between  C2  and  L2,l,  constitutes 
communication  between  master  and  slave.  A  second  approach  would 
allow  each  level  of  processing  to  communicate  directly  with 
adjacent  levels  in  a  multibus  system  (Fig.  3.5).  This 
configuration  will  be  refered  to  as  multibus  and  resembles  a 
pyramidal  communications  network.  The  merits  of  each  arrangement 
will  tie  discussed,  followed  by  a  discussion  of  the  designs  chosen 
for  the  first  and  second  generation  control  experiment. 

In  the  single  bus  system,  interprocessor  communication  takes 
place  through  shared  memory  space.  In  other  words,  there  exists  a 
section  of  memory  which  can  be  accessed  by  two  different  CPU’s. 
All  memory  in  the  single  bus  system  is  addressable  by  the 
housekeeping  processor  while  different  portions  of  that  same 
memory  are  accessible  by  each  of  the  other  processors.  The 
control  processors  cannot  address  each  other’s  memory  space 
directly.  Thus,  the  housekeeper  is  required  to  handle 
synchronization,  avoid  bus  contention  and  perform  data 
transmission. 
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A  typical  interrupt  routine  for  this  single  bus  system  is 
shown  in  Fig.  3.6.  The  local  level  and  housekeeping  processors 
cycle  through  their  respective  interrupt  handlers  faster  than  the 
higher  levels.  The  local  level  processors  retrieve  sensor 
information  via  analog-to-digital  converters  and  place  this 
information  in  memory.  Every  n  local  interrupts  the  housekeeper 
passes  this  information  to  the  appropriate  master  processors  on 
the  component  level.  The  same  is  done  for  the  global  level  every 
nxm  local  interrupts  In  each  case,  n  and  m  are  positive  integers 
whose  values  depend  oi;  the  relative  control  rates  of  the  component 
aid  global  levels.  Control  commands  are  processed  by  till  levels 
and  placed  in  the  common  memory.  After  the  appropriate  number  of 
local  interrupts,  the  housekeeper  sums  and  passes  the  commands 
from  the  higher  levels  to  the  appropriate  local  level  processors 
which  create  the  actuator  commands  through  digital-to-analog 
converters.  The  control  commands  issued  by  the  local  levels  occur 
more  often  than  those  issued  by  the  higher  levels.  This  causes 
the  lower  levels  to  have  higher  control  rates  and  corresponding 
bandwidth.  The  programmable  clocks  keep  all  of  the  processors 
synchronized. 

This  scheme  has  two  important  attributes.  First,  the 
housekeeper  has  the  ability  to  monitor  the  performance  of  all  of 
the  processors.  Secondly,  the  local  level  processors  issue 
control  commands  at  much  higher  rates  than  the  higher  levels. 

This  configuration  has  two  major  drawbacks.  First,  the 
number  of  control  processors  is  limited  by  the  amount  of  memory 
addressable  by  the  housekeeper.  Secondly,  the  rate  of  data 
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transmission  is  reduced  because  all  communication  is  performed 
over  a  single  bus. 

The  multibus  approach  allows  data  transmission  to  be 
performed  directly  by  the  processors  involved.  This  eliminates 
the  extra  step  of  transfer ing  data  to  the  housekeeper.  In  this 
design  the  lower  levels  are  interfaced  to  the  higher  levels 

through  asultiple  porta  and  utilize  shared  memory  for  communication 
(Fig.  3.5). 

Sensor  information  is  passed  to  higher  levels  in  a  fashion 
resembling  a  communication  pyramid  while  control  commands  are 
passed  down  to  the  local  level.  The  typical  interrupt  sequence  is 
illustrated  in  Fig.  3.7.  Data  acquisition  is  performed  on  the 
local  level  as  is  command  output.  Sensor  information  is  pushed  up 
the  pyramid  and  made  available  to  the  higher  levels.  The  higher 
processors  issue  commands  based  on  the  available  information  at 
rates  inversely  proportional  to  the  respective  level.  In  other 
words,  the  higher  the  level,  the  more  processed  information  and 
the  longer  the  processing  delay  resulting  in  slower  control  rates. 
Commands  from  each  level  are  then  passed  back  down  through  the 
pyramid,  summed  with  lower  level  commands  and  issued  to  the 
actuators  from  the  zero,  or  local  level. 

The  multibus  approach  eliminates  some  of  the  limitations 
encountered  with  the  single  bus  configuration.  This  scheme  is  not 
limited  by  the  amount  of  memory  addressable  by  the  global 
processor.  More  information  can  be  passed  at  a  time  because 
several  busses  are  available.  Because  of  more  efficient 
communication  and  limitless  expandability,  the  multibus 


66 


configuration  is  used  in  the  experimental  setup. 


3.3  PRESENT  PROCESSING  CAPABILITY 

The  initial  objective  of  the  microprocessor  development 
program  was  to  implement  the  processors  for  one  element,  i.e.  for 
one  component  and  several  local  processors  (Fig.  3.2).  The 
present  digital  hardware  therefore  consists  of  three  processors 
and  support  peripherals.  A  component  and  local  level  are  achieved 
through  one  master  and  two  slave  processors,  respectively,  which 
communicate  through  shared  memory  space  on  the  slave  boards  (Fig. 
3.8).  The  master  contains  an  INTEL  8088  CPU  which  performs 
sixteen  bit  operations,  synchronized  by  a  five  megahertz  system 
clock,  and  an  8087  NDP,  numeric  data  processor.  The  CPU  board  is 
a  model  8759  built  by  JF  Microsystems.  Each  Ziatech  ZT8830  slave 
board  contains  an  INTEL  8088  CPU.  Both  slave  and  master  memory  is 
addressable  by  the  master  across  an  eight  bit  STD  bus.  But,  each 
slave  can  only  address  its  respective  memory.  In  the  shared 
memory  space,  memory  access  is  performed  through  bus  arbitration 
to  avoid  bus  contention.  Support  peripherals  include 

analog-to-digital  and  digital-to-analog  conversion  boards  for  the 
master  (Analog  Devices:  RTI-1260  and  RTI-1262,  respectively)  and 
an  Intel  iSBX-311  Multimodule  analog  input  board  for  each  slave. 
This  system  is  housed  in  a  Labtech  70  built  by  Laboratory 
Technologies,  Inc.  and  runs  Fortran  86  and  Assembler  86. 

As  will  be  discussed  in  the  next  chapter,  four  actuators  were 


built  and  each  outfitted  with  three  collocated  sensors.  This 
arrangement  calls  for  processing  twelve  measurements  and  producing 
four  command  signals.  To  achieve  this  capability,  each  slave  was 
equiped  with  eight  A/D  channels  and  the  master  was  equiped  with 
four  A/D  and  four  D/A  channels.  This  configuration  deviates 
slightly  from  that  detailed  in  the  previous  section  in  that  the 
all  D/A  functions  take  place  on  the  component  level  through  the 
master  processor.  This  temporary  deviation  from  the  planned 
implementation  was  made  due  to  the  initial  difficulty  of  providing 
for  D/A  capability  at  the  local  slave  processor  level. 

The  interrupt  sequence  used  in  the  control  experiments  is 
displayed  in  Fig.  3.9.  The  interrupt  cycle  starts  with  a  clock 
interrupt  on  the  master  board  which  causes  the  instruction  pointer 
co  exit  Keyboard  communications  and  execute  the  master’s  interrupt 
handler.  The  first  operation  of  the  master  interrupt  handler  is 
to  interrupt  both  slave  boards  by  setting  a  flag  in  the  maskable 
interrupt  control  address  on  each  slave  board.  This  address  is 
located  in  the  shared  memory  space.  The  slaves  digitize  the 
sensor  signals,  pass  some  of  the  digitized  information  to  the 
master,  process  the  remaining  information  and  then  pass  the 
resulting  command  data  to  the  master.  Then  the  slaves  are  idle 
until  the  next  interrupt,  from  the  master.  After  initiating  the 
slaves,  the  master  processes  the  first  set  of  incoming  data,  reads 
the  second  set  of  slave  data,  D/A’s  the  actuator  command  signals 
and  returns  to  its  interrupted  routine  which  performs  keyboard 
communications  and  monitors  the  interrupt  clock.  This  cycle 
repeats  every  five  milliseconds  resulting  in  a  two  hundred  hertz 
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control  rate.  The  second  generation  system  is  presently  under 


development  and  will  provide  greater  processing  power  and  better 
resemble  the  multibus  system  of  Fig.  3.5. 


3.4  SECOND  GENERATION  PROCESSING  NETWORK 

In  the  second  generation  arrangement,  several  new  features 
will  be  added  *  the  original  experimental  hardware.  A  new 
processing  level  is  being  introduced  in  addition  to  a  second  level 
of  communication  links,  or  busses.  Each  current  slave  processor 
will  also  gain  a  high  speed  floating  point  processor  (MATH)  and 
additional  D/A  capability.  In  Figure  3.10,  the  master  now 
represents  a  global  level  processor  with  the  capability  of 
communication  with  the  keyboard  and  the  three  component  level 

processors  (Slaves  C,  D  and  E).  The  three  component  level 

processors  now  contain  A/D,  D/A  and  floating  point  math 

capabilities.  These  processors  can  communicate  with  the  global 
processor  and  the  local  level  processors;  denoted  by  SPH  in  Figure 
3.10.  The  local  level  processors  (SPH)  correspond  to  signal 

processing  hardware  which  perform  real-time,  high  speed  signal 
processing  on  various  sensor  signals. 

This  new  system  will  contain  one  global,  three  component  and 
three  local  level  processors  resulting  in  seven  processors 
operating  simultaneously.  A  total  of  fourty  channels  of  A/D, 
twenty  lght  channels  of  D/A  and  four  independent  communication 
busses  will  be  available  for  control  purposes.  This  hardware 


implementation  more  closely  resembles  the  distributed,  hierarchic 
processing  of  Fig.  3.5  in  that:  the  local  level  processors  operate 
at  rates  much  higher  than  the  component  or  global  levels;  each 
component  level  processor  can  monitor  up  to  four  actuators  within 
its  domain;  and  digital  communication  flows  from  level  to  level  in 


the  processing  pyramid. 


4.0  EXPERIMENTAL  HARDWARE  AND  SOFTWARE  DESIGN 


This  chapter  discusses  the  design  considerations  involved  in 
the  construction  of  both  the  experimental  test  bed  and  the  various 
active  control  components.  A  schematic  of  the  experimental 
hardware  is  shown  in  Fig.  4.1.  First,  the  selection  of  the  test 
article,  or  model  space  structure,  will  be  discussed.  This  will 
be  followed  by  a  discussion  of  the  design  issues  dealt  with  during 
actuator  construction  and  sensor  selection.  Finally,  the 
development  of  digital  and  analog  signal  processing  is  presented. 
In  each  case,  component  development  is  guided  by  the  desire  to 
develop  a  robust  controller  using  components  that  are  capable  of 
functioning  in  the  space  environment  (space-realizable).  In 
addition,  the  experimental  test  bed  is  designed  with  enough 
versatility  to  support  numerous  types  of  control  experiments. 


4. 1  TEST  ARTICLE 

A  free  free  beam  was  selected  as  the  first  test  article. 
These  boundary  conditions  provide  the  most  realistic  model  of  a 
free  floating  space  structure.  A  uniform  beam  was  chosen  as  a 
test  article  because  it  has  a  linear  response  and  can  be  easily 
modeled.  A  twenty  four  foot  brass  beam  was  suspended  from  cables 
arranged  in  a  swinging  parallelogram  fashion  (Fig.  4.2).  This 
parallelogram  suspension  helps  to  suppress  torsional  response  in 
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the  beam  which  may  arise  from  gravity  moment  couples  created  by 
the  motion  of  the  beam  and  actuators.  Except  for  the  suspension 
cables  and  power  cords,  this  test  article  is  free  of  the 
laboratory  frame  and  is  considered  a  realistic  representation  of  a 
large  flexible  space  structure. 

The  desire  was  to  design  a  test  article  whose  standing  and 
propagating  wave  oenavior  were  distinguishable  to  the  control 
processor.  Table  4.1  lists  the  first  sixteen  experimentally 
determined  modal  frequencies  of  this  beam  along  with  its  geometric 
and  material  properties.  The  parameters  of  this  beam  were  chosen 
for  several  reasons.  The  beam  was  made  very  long  and  flexible  to 
test  the  frequency  range  of  the  actuators  and  controller  and  to 
provide  a  test  article  representative  of  a  space  structure.  For  a 
standing  wave  controller  the  beam  dimensions  were  chosen  to  result 
in  many  closely  spaced  modes  at  low  frequencies. 

For  a  propagating  wave  controller,  this  beam  provides 
seemingly  finite  wave  propagation  speeds,  phase  velocities  and 
energy  transport  speeds,  or  group  velocities.  With  a  processor 
control  rate  of  two  hundred  hertz,  waves  associated  with  the  first 
ten  modes  can  be  observed  over  one  hundred  times  during  a  single 
traversal  of  the  beam.  Depending  on  the  wave  speeds  and  the 
closed  loop  damping  in  the  system,  vibrational  motion  may  be 
better  modelled  as  propagating  wave  behavior.  The  theory  behind 
the  wave  mechanics  based  controller  will  be  discussed  briefly  in 
Appendix  B  and  may  be  implemented  in  later  tests.  The  next  step 
is  to  outfit  this  structure  with  the  controller  components:  its 
actuators:  sensors;  and  signal  processing  equipment. 
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4 . 2  ACTUATOR  IMPLEMENTATION 


From  amongst  the  four  types  of  space-realizable  actuators: 
linear  momentum  exchange  devices;  angular  momentum  storage  devices 
(CMG’s);  relative  motion  devices  (hydraulic  actuators);  and 
expelled  mass  momentum  devices  (thrusters);  the  linear  momentum 
devices  were  chosen  for  development.  These  devices  were  modeled 
and  optimized  for  this  application  in  Chapter  2,  but  there  remains 
the  issue  of  hardware  implementation.  The  objective  of  the 
implementation  is  to  create  a  proof-mass  actuator  capable  of 
producing  smooth  and  authoritative  command  forces  or  moments  on 
the  structure  while  also  incorporating  adjustable  passive 
properties. 

Several  types  of  inertial  actuators  have  been  developed  at 
all  having  DC  motors  as  their  electromechanical  driving 
elements.  Translating  proof-mass  actuators,  resembling  audio 
speakers  or  solenoids,  have  been  built.  Here  the  motor  is  a 
translating  DC  motor  with  field  coils  or  permanent  magnets.  These 
devices  proved  to  be  conceptually  simple  but  their  translating 
suspensions  introduced  considerable  friction  and  stiction 
problems.  Rack  and  pinion  devices  incorporate  a  rotational  DC 
motor  and  a  linearly  translating  mass,  but  the  rack  and  pinion 
introduce  backlash  problems.  The  lowest  friction  and  sturdiest 
suspension  proved  to  be  a  pivoting,  proof-mass  arrangment.  In 
this  design,  a  direct  current  motor  rotates  a  reaction  arm.  This 
device  produces  a  torque,  if  the  center  of  mass  (CM)  of  the  arm  is 
placed  at  the  pivot.  If  the  CM  is  offset  from  the  pivot,  a  force 


proportional  to  commanded  motor  torque  is  produced  for  small  arm 


angle  motion.  This  pivoting  configuration  was  chosen  for 

t 

experimental  development.  •  I 

1 

Four  pivoting,  proof-mass  actuators  were  built  each  capable  j 

l 

of  generating  control  forces  on  the  beam  of  up  to  3  Nts.  Figure  j 

4.3  shows  a  typical  actuator,  shown  in  both  the  force  and  torque  1 

cont igurat ion .  1'nis  actuator  design  resembles  a  device  developed 
at  Lockheed  Missiles  and  Space  Company  (Reference  1).  Several 
improvements  have  been  made  over  the  original  Lockheed  design. 

M.I.T.’s  actuators  are  built  from  stock  components  including 
samerium  cobalt,  permanent  magnet  motors.  The  same  actuator  can 
be  adjusted  to  produce  either  torque  or  force,  based  on  the 
arrangement  of  the  masses  on  the  reaction  arm.  Adjustable  passive 
components  were  added  to  the  devices  to  alter  the  passive  dynamics 
of  the  actuator  as  discussed  in  Chapter  2  The  appropriate 
sensors  were  added  to  obtain  full  measurements  of  the  actuator 
states,  in  order  to  provide  full  actuator  observability  and 
compensation  for  the  active  actuator  dynamics. 

The  conceptual  actuator  models  developed  in  Chapter  2  can  be 
directly  applied  to  this  pivoting,  proof-mass  actuator  when  it  is 
used  either  in  the  torque  or  force  configurations.  Fig.  4.4a 
shows  the  model  of  the  torque  configured  actuator  (i.e.  the  CM  of 
the  arm  is  located  at  the  pivot).  In  this  case  a  pure  torque  is 
applied  to  the  beam  and  the  model  of  Chapter  2  can  be  applied  with 
just  a  change  in  notation.  The  following  are  the  governing 
equations  of  motion  for  the  model  shown  in  Fig.  4.4a,  based  upon 
Eq.  2.1. 
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I,«;  4  C80  +  Ko®o  -  C8|  -  Koe,  -  -r  (4.1a) 

IiSi  +  CB,  +  (K,+K0)8,  -  C80  -K08e  -  rp  +  r  (4.1b) 

where 

I0  Rotary  inertia  about  the  center  of  mass  of  the  pivoting 
arm. 

I,  Rotary  inertia  of  the  beam  mode. 

C  Rotary  damping  of  actuator. 

K0  Rotary  stiffness  of  actuator. 

K,  Rotary  stiffness  of  beam  mode, 

r  Control  torque. 

Tp  Disturbance  torque. 


The  optimal  frequency  ratio,  from  the  analysis  in  Section  2.1.1 
using  the  two  step  minimax  criterion,  is 


OPT  "  1  +  P 


'K0/V 

1/2 

K,/I,. 

(4.2a) 


OPT 


P  I0/Ii 


(4.2b) 


Figure  4.4b  shows  the  model  of  the  actuator  in  force 
configuration  (i.e.  the  CM  of  the  arm  is  offset  from  the  pivot  by 
the  length  b>.  The  following  are  the  relations  corresponding  to 
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the  force  configuration,  based  again  on  Eqs.  2.1  and  2.13. 


(lo+mb2)0-  4  C®  -t  K®  +  mbxi  -  r 


(4.3a) 


(Mi4-m)Xi  4  KiXi  4  mb®  =  P 


(4.3b) 


_  f  1  4  (i-p;* 

^APT  i 

~*  *  4  1+0 


(4.4a) 
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2 

mb'" 

I 0  4  mb2 


(4.4b) 


The  new  parameters  for  this  model  are  listed  below. 


I0  Rotary  inertia  about  the  center  of  mass  of  the 
pivoting  arm. 

b  Location  of  the  center  of  mass  with  respect  to  the 
pivot . 

m  Mass  ol  pivoting  arm.  m  -  b,  +  mi 
C  Actuator  damper. 

K  Actuator  stiffness. 

r  Actuator’s  motor  torque.  Control  torque. 

M!  Modal  mass  of  structure. 

In  the  force  configuration,  the  variable  p  accounts  for  the  rati 
of  offset  inertia  to  total  inertia.  Notice  the  similarity  to  Eq 
2.13  for  values  of  p  approaching  unity. 


Table  4.2  lists  the  actuator  specifications  and  components  as 
built.  The  only  actuator  component  that  is  not  space-realizable 
is  the  adjustable  air  damper.  This  is  presently  being  replaced  by 
an  adjustable  electromagnetic  device. 

The  permanent  field  magnet,  DC  motors  of  these  pivoting, 
proof-mass  actuators  have  their  armature  windings  driven  by 
current  amplifiers.  Current  amplifiers  output  a  current  that  is 
proportional  to,  and  in  phase  with,  the  input  voltage  signal 
supplied  by  the  control  processors.  Since  motor  torque  is 
proportional  to,  and  in  phase  with,  armature  current,  the  current 
amplifiers  allow  the  control  designer  to  ignore  the 
electrodynamics  of  the  motor  coils.  Thus,  the  electromagnetic, 
proof-mass  actuator  coupled  to  a  current  source  provides  a  simple, 
easily  modeled,  space-realizable  actuator. 


4.3  SENSOR  IMPLEMENTATION 


The  actuators  and  beam  are  equiped  with  the  necessary  sensors 
for  making  state  measurements  of  the  system.  Note  that  the  states 
include  the  states  of  the  beam  and  the  relative  motion  between  the 
actuators  and  beam.  This  state  information  can  be  obtained  in  two 
ways.  First,  a  state  estimator  can  be  used  to  transform  a 
required  minimum  set  of  measurements  into  estimates  of  all  of  the 
states  of  the  system.  This  is  done  by  using  a  plant  model  in  the 
control  processor  to  relate  sensor  measurements  to  state 
variables.  Secondly,  the  test  article  can  be  outfitted  with  a 
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sufficient  number  of  sensors  to  sense  all  of  its  states  directly 
and  these  state  measurements  can  be  used  to  perform  state  feedback 
control.  In  the  interest  of  implementing  a  simple  robust 
controller  in  the  first  set  of  control  experiments,  the  latter 
approach  was  chosen  and  simple  collocated  state  feedback  was  used. 

The  state  measurements  necessary  to  provide  full  state, 
output  i eeuoocti  t'oi  inis  beam/ accuai-ur  system  are  the  velocity  and 
position  of  the  actuator  relative  to  the  beam  and  the  absolute 
velocity  and  position  of  the  beam  at  each  actuator  location.  To 
maintain  a  space-realizable  test  bed,  only  relative  and  inertial 
motion  sensors  were  considered.  The  measurement  of  the  relative 
motion  between  the  actuator  and  the  beam  was  straight  forward. 
Tachometers  and  potentiometers  are  mounted  on  the  actuator  to 
provide  measurements  of  the  relative  rotational  velocity  and 
position  with  respect  to  the  beam.  Analog  sensors  were  chosen  to 
enhance  low  level  measurements  of  actuator  motion  associated  with 
low  level  vibrations.  This  avoids  the  loss  in  resolution 
encountered  with  discrete  measurement  devices.  The  major 
difficulty  in  state  measurement  stems  from  obtaining  the  beam 
states  in  absolute  coordinates  in  a  way  other  than  measuring 
against  the  reference  frame  of  the  laboratory.  Inertial 
measurements  would  preserve  the  space-realizable  nature  of  the 
experiment . 

Inertial  piezoresistive  accelerometers,  whose  frequency 
response  is  constant  from  zero  frequency  to  2  KHz,  were  used  to 
sense  the  inertial  acceleration  of  the  beam.  The  control 
processors  were  used  to  obtain  velocity  and  position  estimates 


through  a  stabilized  double  integration  algorithm.  This  is  done 
without  the  use  of  a  system  model  and  will  be  discussed  in  Section 
4.4.  Table  4.3  lists  the  specifications  for  each  of  these 
sensors . 

The  final  sensor  configuration  called  for  three  sensors 
(tachometer,  potentiometer  and  accelerometer)  at  each  actuator 
location.  This  provides  four  collocated  state  measurements,  after 
double  integrating  the  accelerometer  signal,  for  use  in  generating 
a  full  collocated  state  feedback  control  command  for  each 
actuatoi . 


4.4  SIGNAL  PROCESSING  IMPLEMENTATION 

The  digital  control  processors  generate  control  commands 
based  upon  the  measured  sensor  data.  In  the  initial  set  of 
experiments,  the  processor  performed  two  functions:  the  stable 
double  integration  of  the  accelerometer  to  yield  the  beam  states; 
and  the  computation  of  force  feedback  commands  from  the  measured 
actuator  states  and  integrated  beam  states.  This  section  will 
describe  the  processing  algorithms  implemented  in  the  processors 
in  the  initial  experiments. 

To  obtain  beam  state  information  (i.e.  velocity  and  position) 
from  a  measured  acceleration  signal,  two  temperal  integrations  are 
required.  The  beam  response,  as  measured  by  the  accelerometers, 
contains  the  true  acceleration  signal  as  well  as  two  noise 


components . 


The  first  noise  component  can  be  modeled  as  a  zero 
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mean,  broad  band  measurement  noise.  Such  noise  is  limited  to  band 
limited  white  noise  after  passing  through  analog,  low  pass, 
antialiasing  filters  (Fig.  4.1)  and  diminishes  in  importance  after 
integration.  The  second  noise  component  consists  of  nonzero  mean 
DC  offsets  in  the  sensors  and  slow  signal  drift  due  to  sensor 
heating  or  deterioration.  Since  the  first  integral  of  nonzero 

vyffseLc  i3  Ixiicai  in  Lime,  anu  toe  second  integral  is  second 
order  in  time,  their  presence  wiii  destabilize  an  exact 
integration  of  the  beam  acceleration.  A  minor  loop  stabilization 
technique  is  required  to  keep  the  displacement  and  velocity 
integrals  within  the  computer’s  digital-to-analog  converter 
windows  (-5  to  +5  volts)  and  avoid  actuator  saturation. 

This  stabilization  of  the  integration  can  be  achieved  by 
output,  (i.e.  displacement  and  velocity)  feedback  and  can  be 
implemented  in  either  analog  hardware  or  digital  hardware  or 
software.  Digital  feedback  stabilization  is  useful  when  working 
at  these  low  beam  frequencies  where  analog  components  become 
unreliable.  Reference  10  discusses  the  derivation  of  appropriate 
state  feedback  algorithms  for  discrete  control. 

Fig.  4.5a  shows  the  open  loop  system  representing  the  simple 
double  integration  which  must  be  stabilized.  The  state  space 
representation  of  the  governing  equations  and  output  of  a  double 
integrator  are 
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where  u  is  the  system  input,  the  accelerometer  signal,  and  y  is 
the  output,  the  position  and  velocity  signals.  Using  any  of  the 
techniques  for  transforming  a  continuous  representation  into  a 
discrete  time  representation  (e.g.  matrix  exponential:  Reference 
10)  gives  the  discrete  time  equations 
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The  sample  period  is  T,  the  processing  delay  is  b  and  state 
feedback  is  implemented  through 
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This  arrangement  is  illustrated  in  Fig.  4.5b. 

The  closed  loop  characteristic  equation  of  the  double 
integrator  is 


z2  i  ( Ti  0 . 5T1 1  ■+  F2)  2 )z  +  T(0.5TF,  -  F2 )  ■+  1 


(4  .b) 


where  '  is  the  discrete  Laplace  operator.  This  second  order 
characteristic  equation  can  be  altered  to  meet  the  objective  of 
providing  stable  integrations  of  the  frequencies  of  interest  while 
rejecting  the  nonzero  mean  noise.  The  solution  is  to  perform  high 
pass  filtering  at  low  frequencies  and  integrations  at  the  higher 
frequencies  of  interest,  in  a  closed-loop  fashion.  The  resulting 
fic4uci^>  transfer  functions  are  shown  in  Fig.  4.6. 

For  a  second  order  system  the  corner  frequency  and  damping 
ratio,  characterizing  the  performance  of  the  filter- to-integrator 
transition  zone,  can  be  defined.  These  characteristics  determine 
the  s -plane  poles  which,  after  being  converted  to  z-plane  poles 
through  Eq.  4.9,  determine  the  desired  roots  of  Eq.  4.8. 

z  ^  es  (4.9) 

Assuming  a  sample  period  of  T  =  0.005  sec,  a  negligible  processing 
delay,  a  corner  frequency  of  1.0  rad/sec  and  a  damping  ratio  of 
0.7071,  the  algorithm  for  the  stabilized  double  integrator  is 
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These  parameters  were  chosen  because  this  is  the  fastest  sample 
period  presently  achievable,  the  stabilization  corner  frequency  is 
half  the  frequency  of  the  lowest  mode  of  interest  (allowing  proper 
integration  of  this  lowest  mode)  and  this  damping  ratio  provides 
minimum  overshoot.  For  a  negligible  processing  delay,  the  output 
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equation  is  identical  to  the  state  equation.  Therefore,  only  the 
state  equation  is  computed  by  the  processor  and  the  velocity  and 
position  values  are  output  from  this  integration  computation 
during  each  interval. 

Fig.  4.6  shows  the  frequency  response  of  this  double 
integration  algorithm  (Eq.  4.10).  Low  frequency  components  of  the 
accelerometer  signal  are  filtered  out  of  the  velocity  output  term 
and  appear  in  the  position  output  term,  multiplied  by  a  constant. 
The  presence  of  this  low  frequency  component  in  the  position 
output  term  is  undesirable.  Therefore  a  single  pole,  high  pass 
digital  filter,  given  in  Equations  4.11a  and  4.11b,  is  used  to 
filter  the  position  term  after  the  integration  routine.  This 
removes  the  nonzero  mean,  low  frequency  error  in  the  position 
estimate. 

xi  -  *i  4  xi  (4.11a) 

z.  .  -  0.99843  z  -  0.00157  x.  (4.11b) 

l  +  l  j  i 

This  algorithm,  Eqs.  4.10  and  4.11,  yields  a  low  frequency 
range  in  which  filtering  of  the  various  input  and  output  terms 
takes  place,  a  middle  frequency  range  in  which  the  integrations  of 
the  accelerometer  signals  take  place  and  a  high  frequency  range  in 
which  the  Nyquist  Frequency  correspondes  to  a  system  zero  due  to 
the  aliasing  caused  by  digitization  of  the  accelerometer  signal. 
The  total  configuration  of  the  integration  process,  consisting  of 
the  analog  antialiasing  filter,  A/I)  and  digital  integrator,  is 
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shown  in  Fig.  4.7.  The  accelerometer  signal  is  low  pass  filtered, 
integrated  twice  and  the  estimated  position  signal  is  high  pass 
filtered. 

Figure  4.8  gives  the  frequency  response  in  magnitude  and 
phase  for  the  estimated  velocity  and  position  signals.  At  low 
frequencies,  the  signals  are  filtered  to  eliminate  nonzero  mean 
noise.  The  transition  region  exhibit4,  the  response  near  the 
selectable  poles  of  the  system.  The  integrating  regions  contain 
tiie  frequencies  corresponding  to  the  beam  modes  of  interest.  The 
highest  corner  frequency  corresponds  to  the  filtering  of  the  zero 
mean  noise  bv  the  antialiasing  filter  and  aliasing  effects 
introduced  by  the  discretization  of  the  accelerometer  signal.  In 
the  current  implementation,  the  poles  are  located  at  1.0  Hz,  for 
Lilt:  integrator  and  high  pass  filter  routines,  and  at  700  Hz,  for 
the  low  pass,  antialiasing  filters.  They  are,  of  course, 
adjustable  i oi  easy  manipulation  of  cornel'  frequencies,  putting 
the  frequencies  of  the  state  signals  of  interest  in  the 
integrating  regions. 

As  a  final  task  of  the  processor,  these  estimated  beam  states 
are  combined  with  the  measured  actuator  states  to  generate  the 
full  collocated  state  feedback  command  to  each  actuator.  The 
computational  task  is  evenly  distributed  among  the  three 
processors.  This  arrangraent  deviates  slightly  from  the  original 
plan  for  hierarchic  control,  as  discussed  in  Section  3.3.  The 
master,  or  component  processor  digitizes  and  double  integrates  two 
accelerometer  signals.  Each  of  the  two  slave,  or  local  processors 
digitizes  five  signals  (including  one  accel erometer  signal), 
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perforins  one  double  integration,  multiplies  states  by  gains  and 
passes  the  information  to  the  master.  Then  the  master  converts 
the  four  command  signals  to  analog  signals  via  its  D/A.  This 
routine  is  repeated  every  5  msec,  yielding  a  200  Hz  sample  and 
control  rate  for  a  twelve  sensor,  four  actuator  system  being 
processed  on  three  commercially  available  microprocessors. 

The  final  hardware  components  in  the  overall  control  system 
are  the  analog  signal  processing  equipment  shown  in  Fig.  4.1  as 
signal  amplifiers  and  filters.  The  preprocessor  signal  processing 
equipment  consists  of  signal  amplifiers,  for  boosting  sensor 
signal  levels  followed  by  low  pass,  antialiasing  filters;  just 
discussed.  Table  4.4  lists  the  specifications  of  the  amplifiers 
and  filters.  Antialiasing  filters  prevent  the  occurence  of  low 
frequencies  which  result  from  digital  misrepresentation  of  high 
frequency  signals.  These  filters  are  second  order  Sallen  and  Key 
circuits.  After  the  command  signals  are  digitally  processed,  they 
pass  from  the  D/A  through  low  pass,  smoothing  filters.  This 
smooths  the  two  hundred  hertz  step  generated  by  the  D/A.  If  not 
filtered,  this  200  Hz  step  component  causes  the  actuators,  and 
high  frequency  beam  modes,  to  buzz.  These  analog  command  signals 
then  pass  to  the  current  amplifiers  to  drive  the  actuators. 

These  hardware  and  software  components  comprise  the 
closed  loop  controller.  Every  feature  is  space-realizable  with 
only  two  exceptions.  First,  the  air  dampers  are  not  functional  in 
space  and  are  being  replaced  by  passive  electrodynamic:  dampers. 
Secondly,  current  amplifier  power  is  not  generated  or  dissipated 
locally.  This,  of  course:,  could  be  done  with  appropriate  space 


based  power  generating  equipment. 

With  this  controller  hardware,  several  tests  were  performed 
to  demonstrate  the  use  of  some  collocated  state  feedback 
techniques  and  verify  that  this  equipment  functions  properly. 
These  experiments  will  be  described  in  the  next  chapter. 
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5.0  CONTROL  EXPERIMENTS 


A  preliminary  set  of  simple  collocated  state  feedback 
experiments  was  deemed  essential  to  verify  that  the  closed-loop 
controller  functions  properly.  The  overall  goals  of  the 
experiments  detailed  in  this  chapter  were  twofold:  to  verify  the 
controllability  of  structural  modes  using  an  inertially  actuated, 
space-realizable  control  i  'stem;  and  to  establish  the 
fund ioruil i t y  of  the  versatile  testbed  upon  which  future  testing 
of  the  full  hierarchic  controller  can  be  performed. 

These  tests  had  several  specific  objectives.  First,  they 
were  used  to  identify  any  necessary  hardware  changes.  This 
includes  evaluation  of  the  digital  processing  performance,  analog 
signal  processing  equipment  and  sensor  performance.  Secondly,  the 
experiments  established  the  degree  to  which  the  actuators  can  be 
tuned,  according  to  the  rules  derived  in  Chapter  2.  Third,  the 
test  would  serve  as  a  verification  of  the  controllability  of 
structural  modes  using  inertial  actuation  and  collocated  state 
feedback.  Lastly,  these  experiments  help  identify  some  of  the 
experimentally  encountered  limitations  involved  in  state  feedback, 
structural  control. 

Four  series  of  tests  were  planned  to  achieve  these  specific 
objectives.  The  tests  were  chosen  in  a  sequence  that  would  first 
demonstrate  the  effectiveness  of  strictly  active  and  then  strictly 
passive  damping  and,  then  demonstrate  the  effectiveness  of  a 


combination  of  both  damping  techniques. 


In  the  conduct  of  each 


test,  any  limitations  on  the  system’s  performance  were  noted  and 
appropriate  solutions  or  explanations  are  presented. 

The  first  experiment  tests  a  single  input,  single  output 
(SISO)  active  controller.  In  this  test,  one  actuator  was  used  to 
control  structural  modes  using  only  collocated  structural  rate  in 


the  feedback  loop.  The  second  test  improves  on  this  technique  by 
adding  co u oca tea  actuator  rate  to  the  feedback  loop.  In  the 
second  test,  two  actuators  were  used  to  control  several  modes, 
each  using  both  collocated  structural  and  actuator  rate 
measurements.  Thirdly,  the  passive  components  of  the  actuators 
were  tuned  to  different  modes  following  the  tuned  vibration 
absorber  methods  detailed  in  Chapter  2,  and  the  performance  of  the 
passive  system  was  evaluated.  Lastly,  both  active  and  passive 
techniques  were  employed  to  damp  the  same  structural  modes. 

In  the  following  section,  the  preliminary  design  of  the 
experiment,  is  detailed  (Section  5.1).  Sections  5.2  through  5.5 
present  the  experimental  results.  In  each  of  these  sections,  the 
expected  results,  actual  results  and  degree  of  achievement,  or  the 
limitations  encountered,  are  discussed. 


5.1  DESIGN  CONSIDERATIONS  FOR  EXPERIMENTAL  PHASE 

There  are  several  questions  bearing  on  the  experimental 
design  that  must  be  addressed  by  analysis  prior  to  testing.  These 
questions  are  concerned  with  how  the  controller  hardware  should  be 
arranged  with  regard  to  the  test  structure  and  what  stability 
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The  four  such 
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limitations  can  be  expected  during  testing, 
questions  discussed  in  this  section  are: 

-  Under  what  hardware  and  configuration  conditions  are  the 
beam  modes  controllable  (i.e.  arbitrary  pole  placement 
through  feedback  selection  given  a  linear,  time-invariant 
system)? 

-  What  are  the  primary  energy  dissipation  mechanisms  in  this 
system? 

-  Where,  on  the  beam,  should  the  four  actuators  be  placed  in 
order  that  maximum  control  performance  is  achieved? 

-  What  stability  limitations  exist  when  using  only  beam  state 
feedback  (i.e.  not  compensating  for  actuator  dynamics)? 

An  uridei standing  of  the  answers  to  these  questions  helps  determine 
an  appropriate  sequence  of  experiments. 

The  question  of  prime  importance  is  that  of  controllability. 
The  dynamics  of  the  beam  will  be  modeled  as  a  single  structural 
mode.  Can  a  S1S0  controller,  using  an  inertial  actuator  and  a 
collocated  beam  rate  measurement,  exert  control  to  the  extent  that 
there  is  an  increase  in  that  mode’s  structural  damping? 
Furthermore,  does  full  collocated  state  feedback  (including 
actuator  states)  provide  the  system  with  the  capability  to  perform 
arbitrary  pole  placement? 

The  question  of  controllability  is  answered  by  examining  the 
'Controllability  Matrix’  (Reference  10).  If  the  matrix  has  full 
rank  and  the  system  is  linear  and  time-invariant,  the  system  is 
fully  controllable  and  the  closed-loop  poles  may  be  arbitrarily 
placed  through  the  feedback  of  all  of  the  system  states.  If  the 
determinant  of  this  matrix  is  nonzero,  for  a  single  input  system, 
it  is  of  full  rank. 


In  order  to  examine  the  controllability,  the  system  matrices 


for  a  pivoting,  proof-mass  actuator,  coupled  to  a  structural  mode, 
must  be  derived  in  appropriate  form.  The  reaction  equations  of 
motion  for  a  pivoting  proof-mass  actuator,  used  as  a  force 
actuator  coupled  to  a  single  structural  mode,  are  repeated  from 
Chapter  4  as 

tl0+mb')8  +  C0  +  K8  +  mbXi  =  r  (5.1a) 
(Mj+m)Xj  +  CjXi  +  Ki Xj  +  mb®'  =  P  (5.1b) 


Equivalent  viscous  structural  damping  (Ct)  is  introduced  to 
determine  its  effect  on  controllability.  Since  pole  placement 
represents  only  homogeneous  dynamics,  P  can  be  set  to  zero.  The 
control  torque  (r)  is  a  function  of  state  variables  and  is 
therefore  homogeneous.  Transforming  these  equations  into  state 
variable  torn:, 

z  -  Az  -*  Bu  (5.2a) 
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yields  the  Fundamental  Matrix  (A)  and  the  Control  Effectiveness 
Matrix  (B)  where 


90 


a  -  m  -  KMj+m)  (5.2c) 

I  =  I0  +  mb2  (5. 2d) 


The  controllability  matrix  is  defined  as 
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A  linear,  time-invariant  system  (LTI),  as  defined  in  Eq.  5.2a,  is 
controllable  if  and  only  if  the  columns  of  P  span  the  dimension  of 
the  state  space  representation  (i.e.  the  determinant  of  P  is 
nonzero  for  this  single  input  system).  The  determinant  of  P  is 


det[Pj 


(mbT / a 
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(5.4) 


Thus,  the  following  conclusions  can  be  made  for  a  single  DOF 
damped  structural  mode  acted  on  by  a  pivoted,  proof-mass  actuator. 
To  be  controllable,  the  structural  mode  must  be  elastic  (T  >  0),  a 
react  on  mass  must  exist  (m  >  0)  and  that  mass  must  be  offset  from 
the  pivot  to  generate  force  (b  #  0).  The  denominator  (a)  requires 
that  I  ,  M,  m  and  b  be  finite.  The  latter  is  a  trivial  result  in 
that  the  beam’s  modal  mass  must  not  be  immovable  and  the  actuator 
components  must  allow  it  to  pivot.  No  restrictions  are  placed  on 


actuator  stiffness  and  damping  or  on  passive  modal  damping  in  the 


Since  this  single  structural  mode  is  controllable,  even  in 
the  absence  of  passive  energy  dissipation  mechanisms,  a  mechanism 
associated  with  the  active  actuator  must  exist  which  enables 
dissipation  to  be  introduced. 

To  identify  the  source  of  dissipation,  examine  the  work 
balance  at  the  actuator.  The  work  performed  on  a  system  equals 
that  system's  change  in  energy.  The  work  on  the  structure’s  modal 
mass  *M’  is 


Work  -  |  Fdx  -  }  mb^8  d8  (5.5) 

The  work  is  done  as  a  result  of  relative  motion  between  the 
pivoting  part  of  the  actuator  and  structure.  This  relative  motion 
causes  the  motor  coils  to  travel  through  the  magnetic  field  of  the 
permanent  magnet  thur  generating  a  back  EMF  proportional  to  the 
relative  velocity  (8). 

In  steady-state  oscillatory  motion,  the  product  of  this 
voltage  and  the  motor's  command  current  gives  the  power 
transmitted  to  the  system  as 

P  =  \-l-cosr  (5.6) 

o 

If  the  phase  (f)  between  voltage  and  current  is  other  than  90  or 
© 

270  ,  some  energy  is  supplied  to,  or  drained  from,  the  structural 
system.  Thus,  the  presence  of  the  actuator  effects  a  change  in 
system  damping  by  this  energy  transfer  process.  When  increasing 


the  damping  in  the  system,  the  controller  acts  as  a  generator 


converting  mechanical  energy  into  electricity  which  is  eventually 

c  o 

dissipated  through  heat.  If  the  phase  is  exactly  90  or  270  ,  no 
energy  transfer  occurs  and  the  system  experiences  only  a  change  in 
stiffness. 

Since  the  work  on  the  structure  is  equal  to  the  integral  of 
force  times  the  structural  displacement,  the  structure  must 
exhibit  motion  at  the  actuator  location  in  order  to  effect  energy 
transfer.  The  structural  modal  mass  must  therefore  be  finite;  the 
same  requirement  derived  to  insure  controllability.  Since 
apparent  modal  mass  is  determined  through  Eq.  5.7,  where  the  mode 
shape  O')  is  normalized  to  the  position  where  the  modal  mass  is 
desired  (actuator  location),  if  the  actuator  is  placed  at  a  node, 
■f  cannot  be  normalized  to  that  location  and  ’M  ’  is  infinite. 


M  -  J  m  ( x )  ( x  )  dx 

1  0  1 


(5.7; 


This  line  of  reasoning  raises  the  question  of  allowable  and 
optimal  actuator  placement.  Clearly,  Eq.  5.7  suggests  the 
actuator  should  not  be  placed  at  a  node.  To  maximize  the  work 
done  on  the  structure,  the  force  from  the  actuator  should  be 
exerted  through  the  greatest  displacement.  Since  the  mode  shape  t 
represents  the  open-loop,  spacial  variation  of  displacement  of  the 
mode,  this  function  might  be  used  to  determine  the  most  effective 
actuation  locations.  More  rigorous  solutions  to  optimal  actuator 
placement  are  discussed  in  Reference  11  and  alternative  placement 


criteria  are  presented  in  Reference  13.  For  a  single  mode,  the 
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actuator  would  clearly  be  placed  at  the  location  of  maximum  modal 
displacement .  lor  multiple  mode  control,  both  the  sum  and  product 
of  the  squares  of  the  modes  might  be  considered.  Figure  5.1a 
shows  the  sum  of  the  mode  shapes  squared  for  a  uniform  free-free 
beam,  defined  by 


(11  -  1  i‘‘ 
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and  Fig.  5.1b  shows  the  product  of  the  mode  shapes  squared, 
defined  bv 


(2)  17.  ,  ? 

l-l  l 


(5.9) 


Figure  5.Ja  indicates  positions  on  the  beam  where 
displacement  is  at  a  minimum  and  maximum  for  the  combined  motion 
of  the  modes  listed.  The  zero  levels  in  Figure  5.1b  indicate 
positions  on  the  beam  where  at  least  one  mode  is  uncontrollable. 
Note  the  relative  difficulty  in  locating  actuators  for  control  of 
the  higher  modes.  For  an  unrestrained  free-free  beam,  both 
criteria  indicate  that  for  ail  modes  the  tip  is  the  best  place  to 
locate  an  actuator,  so  two  of  the  four  actuators  will  be  placed  at 
the  tips.  Other  relative  maxima  in  the  actuator  effectiveness 
curves  for  the  higher  modes  occur  c  approximately  the  0.08  and 
0.25  length  fractions  from  the  tip.  For  these'  experiments,  the 
remaining  two  force  actuators  a*~e  placed  at  a  quarter  of  t  he 
length  i n  from  each  tip.  The  actuator  locations  are  indicated  by 
asterisks  at  the  bottom  of  Figures  5.1a  and  5.1b. 
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Once  a  controllable  system  is  built  and  actuators  are  placed, 
an  investigation  can  be  made  into  stability  limits  for  various 
state  feedback  methods.  Several  assumptions  are  made  for  this 
analysis.  First,  these  stability  limitations  will  only  be  derived 
for  a  single  actuator  controlling  a  single  mode  using  collocated 
beam  and  actuator  states.  Secondly,  state  estimator  dynamics  are 
not  considered  since  states  are  assumed  to  be  measured  directly. 
Finally,  the  actuator  is  assumed  to  have  a  frequency  that  is  much 
lower  than  the  frequency  of  the  mode  of  interest.  In  other  words, 
it  is  a  ’weak’  actuator  not  tuned  to  the  frequency  of  the  mode  of 
interest.  The  model  shown  in  Fig.  2.2  is  used.  The  feedback  laws 
that  will  be  investigated  are  of  the  form 

t  -  -FjXj  -  F2S  -  F3X,  -  F4»  (5.10) 

where  Xi  and  \t  are  the  beam  states  and  ®  and  8  are  the  actuator 
states.  The  discussion  will  be  limited  to  the  rate  feedback 
terms,  F^  and  F^  (i.e.  F^  =  Fg  =  0). 

The  most  obvious  approach  to  increasing  the  beam  damping  is 
to  feed  the  actuator  the  collocated  velocity  of  the  beam.  As  seen 
in  the  root  locus  of  Figure  5.2a  for  increasing  beam  rate 
feedback,  this  introduces  damping  into  the  high,  predominantly 
beam  mode  while  destabilizing  the  lower,  actuator  mode.  The 
opposite  is  true  when  only  absolute  actuator  rate  is  fed  back 
(i.e.  the  upper  beam  mode  is  destabilized  and  the  actuator  mode 
becomes  damped).  If  relative  actuator  rate  is  fed  back,  the 
system  responds  as  if  a  passive  damper  was  placed  between  the  beam 


and  actuator.  Thus  in  the  absence  of  passive  damping,  feedback  of 
either  the  absolute  beam  or  absolute  actuator  rate  alone  is 


destabilizing. 

Two  options  are  therefore  available  to  stabilize  the 
actuator/beam  system.  First,  if  passive  damping  is  present  in  the 
actuator,  both  the  lower  and  upper  modes  are  initially  damped  as 
shown  oy  tne  locations  ot  the  open-loop  poles  in  Figure  5.2b. 
Upon  feeding  back  small  amounts  of  beam  rate,  the  low  (actuator) 
mode  will  be  destabilized,  but  still  have  net  positive  stability 
due  to  the  passive  damping.  With  further  increases  in  beam  rate 
feedback,  the  system  will  become  unstable  (point  B).  An 
experiment  based  on  this  single  rate  feedback  will  be  discussed  in 
Section  5.2.  As  a  second  approach,  both  beam  and  actuator  rate 
can  be  fed  back,  in  effect  compensating  for  actuator  dynamics. 
Actuator  rate  feedback  stabilizes  the  dynamics  of  the  actuator 
which  were  destabilized  by  the  use  of  beam  rate  feedback,  as  shown 
in  Fig.  5.2c.  Therefore,  a  combination  of  both  beam  and  actuator 
velocity  feedback  is  needed  to  maintain  system  stability  while 
generating  significant  increases  in  beam  damping. 

With  the  four  actuators  placed  at  the  tips  and  the  quarter 
length  locations,  a  test  was  performed  using  one  of  these  four 
distinct  S1S0  controllers  to  demonstrate  the  stability  limitation 
inherent  in  single  collocated  beam  rate  measurement  feedback.  The 
actuator  dynamics  are  compensated  through  active  means  in  the  MISO 
tests  '.iscussed  in  Section  5.3. 


5.2  SINGLE  INPUT,  SINGLE  OUTPUT  (SISO)  COLLOCATED  CONTROLLER 


The  first  test  involved  controlling  individual  bean  modes  by 
feeding  back,  to  one  tip  actuator,  the  corresponding  collocated 
beam  velocity.  This  test  was  performed  to  investigate  the  onset 
of  the  instability  discussed  in  the  previous  section.  The  tuned 
frequency  of  the  force  actuator  was  set  well  below  the  frequency 
of  the  seventh  mode.  This  situation  in  which  the  actuator 
frequency  is  much  less  than  the  beam  frequency  will  hereafter  be 
refered  to  as  a  ’weak*  actuator  case.  The  closed-loop  dynamics  of 
the  seventh  mode  were  chosen  for  investigation  so  as  not  to 
require  the  weak  actuator  frequency  to  be  set  so  low  that  friction 
effects  in  the  mechanism  dominate  its  behavior. 

As  discussed  above  in  Section  5.1,  increasing  beam  rate 
feedback  should  monotonically  increase  the  damping  in  the  higher, 
beam  mode  and  destabilize  the  lower,  actuator  mode.  At 
sufficiently  high  gains,  the  feedback  will  overcome  the  passive 
damping  in  the  low  mode  and  an  instability  should  appear.  This 
test  was  run  at  various  negative  beam  velocity  feedback  gains. 
Figure  5.3  shows  the  time  response  of  the  beam  acceleration  as  a 
result  of  an  initial  displacement.  The  gain  is  set  to  that  value 
which  first  exhibited  the  low  mode  instability.  In  this  plot  of 
the  accelerometer  response,  the  lightly  damped  5.2  Hz  beam 
response  initially  decays.  The  very  slightly  unstable  lower  mode, 
which  primarily  contains  actuator  motion,  eventually  drives  the 
beam  unstable.  This  verifies  the  trend  of  the  root  locus  analysis 
performed  in  Section  5.1  for  feedback  of  beam  rate  only. 
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The  factor  that  limited  the  achievable  damping  in  this  single 
rate  feedback  case  was  the  occurance  of  this  instability.  This 
can  be  corrected  in  two  ways:  a  passive  damper  can  be  mounted  on 
the  actuator  to  stabilize  the  low,  ’actuator*  mode,  or  actuator 
rate  feedback  can  be  introduced  through  feedback  of  the  actuator 
tachometer  signal.  The  latter  method  was  chosen  and  a  second  set 
of  tests  was  pert'ormeo. 

5.3  MULTIPLE  INPUT,  SINGLE  OUTPUT  (MISO)  CONTROLLER 

The  next  set  of  tests  used  collocated  beam  and  actuator  rate 
feedback  to  stabilize  the  low  mode  instability  and  thus  enable 
higher  levels  of  active  beam  mode  deaping.  In  these  tests,  two 
force  actuators  were  used;  one  at  each  end  of  the  beam.  Note 
that  although  a  total  of  eight  state  variables  were  available, 
beam  and  actuator  rate  and  position  at  each  tip,  each  actuator  was 
only  fed  back  its  two  collocated  state  measurements.  Thus,  two 
independent  MISO  processors  were  used.  Again,  ’weak’  actuators 
were  used  to  control  individual  beam  modes. 

Results  of  tests  run  on  modes  seven  through  fourteen  are 
summarized  in  Table  5.1.  These  modes  were  chosen  for 
investigation  since  their  frequencies  were  high  enough  that  the 
actuator  friction  was  again  not  dominant,  and  due  to  a 
signal-to-noise  difficulty  in  the  accelerometer  which  will  be 
discussed  below.  Figures  5.4a  through  5.11b  show  the  open-loop, 
free  decays  and  closed- loop,  controlled  decays  of  these  modes. 
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Table  5.1  lists  the  open  and  closed-loop  frequencies  and  damping 
ratios,  measured  using  logarithmic  decrement  methods,  for  each 
mode.  Figures  5.12a  and  5.12b  show  the  results  of  a  frequency 
sweep  of  the  uncontrolled  and  controlled  beam,  respectively.  The 
sweeps  were  performed  at  constant  feedback  gains.  These  results 
show,  on  the  average,  a  factor  of  4  to  6  decrease  in  the  forced 
vibration  response  of  the  beam.  This  is  in  substantial  agreement 
with  the  data  taken  from  the  transient  responses  which  show  an 
average  improvement  in  damping  ratio  (?)  of  7.25. 

Two  major  difficulties  arose  in  these  tests.  First,  at  high 
feedback  gains,  a  system  mode  at  42  Hz.  was  driven  unstable.  With 
the  computers  sampling  two  hundred  times  a  second,  42  Hz.  lies 
near  the  frequency  where  the  processors  must  be  considered  as 
discrete  time  controllers.  Furthermore,  phase  shifts  due  to 
filter  rolloffs  also  occur  in  this  frequency  range  causing  a 
lurther  reduction  in  phase  margin  vsee  Table  4.4). 

The  second  difficulty  prevented  the  control  of  the  first  six 
beam  modes.  At  these  modal  frequencies,  the  accelerometer  signal 
contained  too  low  of  a  signal-to-noise  ratio.  In  order  to 
maintain  a  constant  level  of  accelerometer  signal  in  low  modes, 
the  displacement  must  increase  by  the  square  of  the  frequency 
ratio  between  the  higher  and  lower  modes.  Shaker  power  and 
suspension  limitations  prevented  such  large  amplitude  motions  of 
these  lower  modes.  In  addition,  the  actuator  .frequency  of  the 
’weak’  actuator  was  in  close  proximity  to  the  modal  frequencies  of 
these  lower  modes.  This  would  require  a  more  complicated  state 
feedback  algorithm  employing  additional  displacement  feedback. 
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5.4  TESTS  USING  PASSIVE  VIBRATION  ABSORBERS 


As  described  in  Chapter  2,  passive  damping  of  the  beam  can  be 
augmented  by  adding  spring-mass-damper  passive  vibration 
absorbers.  The  proof-mass  actuators  were  designed  with  variable 
springs  and  dampers  so  that  they  could  also  fulfill  this  function. 
Each  actuator  was  equiped  with  a  variable  spring  and  adjustable 
damper  to  allow  passive  tuning.  Tuning  refers  to  the  adjustment 
of  passive  absorber  parameters  to  achieve  some  desired  system 
performance.  In  these  tests,  absorbers  were  tuned  to  various 
individual  beam  modes.  The  performance  of  the  absorbers  was 
assessed  in  free  decay  tests  and  in  steady-state  sine  sweep  tests. 
This  section  discusses  the  expected  and  actual  test  results  along 
with  the  probable  causes  of  disagreement  and  the  limitations 
encount  ered. 

The  seventh,  eighth  and  nineth  modes  of  the  beam  were  damped 
using  passive  vibration  absorbers  (Table  4.1).  These  modes  were 
selected  to  enable  a  comparison  with  results  from  the  active 
control  tests  presented  in  Section  5.3.  The  first  step  involves 
determining  the  achievable  amount  of  damping  using  a  single 
absorber  tuned  to  a  single  mode.  From  the  analysis  of  Chapter  2, 
predictions  can  be  made  as  to  the  best  theoretically  achievable 
performance.  The  dynamic  magnification  factor  for  the  magnitude 
of  the  beam  response  at  resonance,  as  a  function  of  modal  damping 
ratio  and  static  response,  is 
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Equation  2.15,  repeated  below,  presents  this  sane  magnification 
factor  as  a  function  of  absorbei — to-modal  mass  ratio  (0)  for  the 
optimally  tuned  absorber. 
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(2.15) 


lor  these*  passive  tuning  techniques,  these  two  equations  can  be 
combined  to  predict  the  highest  achievable  modal  damping  ratio  for 
an  otherwise  undamped  beam 
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The  model  of  the  pivoted,  proof-mass  actuator  must  be  put  in 
terms  of  the  vibration  absorber  of  Chapter  2.  For  this  purpose, 
the  absorber  arms  are  assumed  massless  (1=0)  and  the  fraction  of 
the  absorber  mass  associated  with  the  nonpivoting  base  is  added  to 
the  modal  mass  of  the  beam.  For  the  modes  tested,  Table  5.2  lists 
the  absorber  position  on  the  beam,  effective  modal  mass,  mass 
ratio  (0),  optimal  frequency  ratio  (6QPT)  and  predicted  maximum 
modal  damping  ratio.  In  each  case,  the  proof  mass  is  0.165  kg. 


An  attempt  was  made  to  experimentally  tune  the  absorbers  to 
the  values  shown  in  Table  5.2,  but  this  attempt  encountered 
certain  limitations.  In  tuning  each  individual  absorber  to  a 


particular  mode,  a  two  step  iteration  was  repeated  until  the  best 
results  were  obtained.  The  tuning  process  started  with  a  'weak* 
absorber  (low  absorber/beam  frequency  ratio  6).  For  this  weak 
absorber  spring,  the  damper  was  adjusted  to  yield  the  highest 
achievable  modal  damping  ratio.  Figure  5.13a  shows  the  achieved 
damping  in  the  eighth  mode  for  such  a  weak  actuator  (i.e.  for  a 
frequency  ratio  of  5  -  2.3  Hz/5.711  Hz  =  0.387).  In  this  case,  a 
tip  aosorber  was  tuned  to  the  eighth  mode.  The  next  step  was  to 
increase  the  frequency  ratio  and  readjust  the  damper.  This  two 
step  iteration  was  repeated  until  maximum  damping  was  obtained. 
Table  5.3  lists  the  intermediate  results  of  this  tuning  process 
(i.e.  the  largest  damping  ratios  achieved  at  each  frequency  ratio 
for  tuning  one  tip  absorber  to  the  eighth  mode).  The  free 
vibration  decays  are  shown  in  the  corresponding  figures. 

The  two  other  absorbers,  located  at  the  one  quarter  length 
and  at  the  other  tip,  were  individually  tuned  to  modes  seven  and 
nine,  respectively.  The  achieved  modal  damping  values,  for  each 
individual  absorber  acting  alone,  are  listed  in  Table  5.4  and 
Figures  5.14  and  5.15  show  the  modal  free  decays.  Figure  5.16 
shows  a  frequency  sweep  through  modes  seven,  eight  and  nine  with 
one  absorber  tuned  to  each  mode.  This  can  be  compared  to  Fig. 
5.12a  for  the  baseline,  open-loop  beam  response  and  to  Fig.  5.12b 
for  the  closed- loop  actively  controlled  beam  response.  By 
comparing  the  analytic  predictions  of  Table  5.2  and  the  measured 
results  of  Table  5.4,  it  is  clear  that  overall,  the  tuned 
absorbers  were  able  to  achieve  about  40*  of  the  analytically 


predicted  level.  Mode  eight  was  as  high  as  84*  while  mode  nine 


was  the  lowest  at  34%. 


Unmodeled  friction  present  in  the  absorbers,  the  proximity  of 
the  absorber  to  beam  nodes  and  the  difficulty  of  tuning  the 
absorbers  close  to  resonance  presented  the  major  problems  in 
performing  this  portion  of  the  experiment.  The  effects  of 
friction  were  twofold.  As  the  frequency  ratio  was  increased 
towards  the  optimal,  defined  in  Eq.  2.13,  the  damper  setting 
needed  to  create  the  optimal  results  decreased  (Fig.  5.17a). 
Above  a  frequency  ratio  of  6  =  0.9,  the  damping  associated  with 
the  friction  in  the  absorber  mechanism  was  greater  than  that 
needed  to  yield  the  maximum  amount  of  beam  damping.  The 
horizontal  lines  represent  the  minimum  achievable  damper  levels. 
These  levels  are  set  by  the  amount  of  friction  in  the  absorbers. 
The  equivalent  nondimensional  damping  due  to  the  friction  was 
determined  by  inferring  the  necessary  nondimensional  damping 
needed  to  produce  the  highest  beam  damping  ratios  achieved  in  the 
experiments.  The  horizontal  lines  correspond  to  the  absorber 
tuned  to  the  eighth  and  nineth  modes.  Figure  5.17b  shows  that  up 
to  this  friction  limit,  the  experimentally  achieved  damping  is  in 
fairly  close  agreement  with  the  analytic  prediction. 
Discrepencies  are  due  to  measurement  errors  and  the  fact  that  the 
beam  possesses  some  level  of  material  and  aerodynamic  damping.  To 
compensate  for  the  latter,  the  open-loop  damping  ratio  from  Table 
5.1  was  added  to  the  analytic  prediction  curve.  Figure  5.18  shows 
the  value  of  the  nondimensional  damper  (u)  and  equivalent  absorber 
damper  strength  (nP)  versus  mass  ratio.  These  parameters  are 
defined  in  Eqs.  2.4a  through  2.4h.  Equivalent  damper  strength 


increases  with  mass  ratio.  This  indicates  that  the  absorber 
dampers  are  more  easily  tuned  at  high  mass  ratios.  For  the  low 
mass  ratios  used  in  the  experiment,  as  6  approached  the 
damping  required  fell  below  the  minimum  set  by  friction  and  the 
experimental  results  and  analytic  predictions  departed. 

Fig.  5.15  shows  the  free  vibration  decay  for  the  case  where 
the  damping  in  the  vibration  absorber  is  in  fact  dominated  by  the 
friction.  The  linear  decay  of  the  response  shown  is 
characteristic  of  friction  damping.  A  second  noticable  feature  of 
the  frictional  damping,  as  opposed  to  linear  viscous  damping,  is 
that  at  low  vibration  levels,  the  friction  would  seize  the  motion 
of  the  absorber.  This  eliminates  motion  across  the  damper  and 
therefore  stops  any  further  damping  of  the  beam  motion  (Fig. 
5.13J1.  These  difficulties  can  be  reduced  through  the  use  of  a 
lower  friction  suspension.  However,  since  all  rotational 
mechanisms  contain  some  friction,  this  represents  a  lower  limit  of 
the  motion  amplitude;  for  which  such  absorber  designs  are  useful. 

The  proximity  of  beam  nodes  to  the  absorber  locations 
strongly  affect  the  amount  of  achievable  damping  (Reference  11). 
When  absorbers  are  located  other  than  at  the  tip,  the  apparant 
modal  mass  of  the  beam  is  greater  causing  smaller  absorber- to-beam 
modal  mass  ratios  for  a  fixed  absorber  mass.  The  smaller  the  mass 
ratio,  the  lower  the  modal  damping  achievable  (Eq.  5.12).  A 
second  effect  occurs  due  to  the  change  in  the  open-loop  beam 
dynamics  caused  by  the  addition  of  the  absorber  mass.  When  the 
concentrated  mass  of  the  absorber  is  added  to  a  beam,  nodes  in 
proximity  to  that  mass  are  drawn  towards  it  causing  a  reduction  in 
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modal  amplitude  at  the  location  of  the  absorber,  and  a  subsequent 
further  increase  in  the  apparant  modal  mass.  For  these  two 
reasons,  actuators  and  absorbers,  which  were  placed  in  the 
interior  of  the  beam,  were  less  effective  than  those  that  were 
placed  at  the  tips.  The  former  is  due  to  the  fundamental  physics 
of  the  mode  shape  while  the  latter  can  be  diminished  by  reducing 
the  nonpivoting  mass  of  the  absorbers,  or  by  using  rotational 
torque  absorbers.  These  arguements  as  to  why  absorbers  have 
diminished  effectiveness  are  equally  applicable  to  the  selection 
of  locations  lor  active  control  actuators,  which  also  rely  on 
local  beam  motion  to  be  effective  (Reference  11). 

The  third  difficulty  encountered  in  the  experimental 
imp] ementat 1  on  of  the  absorbers  stei^>  from  the  fact  that  tuned 
absorbers  are  very  close  in  frequency  to  the  beam  mode.  This, 
along  with  the  fact  that  the  optimal  damper  setting  required  (as 
the  frequency  ratio  6  approaches  optimal1!  decreases,  causes 
absorber  motion  to  become  extremely  large.  In  many  cases,  this 
response  violates  small  angle  approximations  in  the  analysis. 
Furthermore,  the  physical  constraints  of  small  angle  absorber 
motion  are  encountered  and  absorber  saturation  occurs.  Increasing 
the  offset  inertia  (mb)  of  the  absorber  will  lower  this  response 
amplitude  to  some  extent. 
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5.5  DAMPING  USING  ACTIVE  ACTUATORS  AND  VIBRATION  ABSORBERS 


A  final  test  was  performed  using  a  combination  of  both  active 
feedback  and  passively  tuned  absorbers.  In  this  experiment,  the 
actuators  located  one  quarter  of  the  length  in  from  the  ends  of 
the  beam  were  tuned  to  the  ninth  mode  as  passive  vibration 
absorbers.  At  the  same  time,  ’weak’  actuators,  located  at  the 
tips  of  the  beam,  were  used  to  actively  control  the  same  mode. 
Figure  5. 19d  shows  the  time  response  results  to  initial 
displacement  of  this  configuration.  Figures  5.19a,  5.19b  and 
5.19c  are  riineth  mode  free  decays  for  the  open-loop,  actively 
controlled  and  passively  damped  systems  repeated  from  the  previous 
figures  for  comparison. 

As  these  two  forms  of  damping  were  used  simultaneously,  the 
problems  associated  with  each  were  compounded.  Limited  response 
amplitude  due  to  insufficient  excitation  and  actuator  stiction, 
prevented  accurate  analysis  of  the  overall  effectiveness  of  this 
combined  damping  scheme.  The  initial  response  in  Fig.  5.19a  is 
well  over  an  order'  of  magnitude  larger  than  that  in  Fig.  5.  19d. 
In  actuality,  the  beam  has  become  too  heavily  damped  for  the 
exciting  shaker-  to  create  a  significantly  large  response 
displacement.  In  other  words,  the  beam  has  become  so  damped  that 
its  response  is  within  the  measurement  noise  and  the  deadband  of 


the  absorbers. 


6.0  CONCLUSIONS 


The  work  presented  represents  the  present  state  of 
development  of  the  finite  element  controller.  The  following  is  a 
summary  of  the  goals  achieved  in  this  phase  of  research  and  the 
experiments  planned  for  the  immediate  future. 

Throughout  this  phase  of  research,  all  components  of  the 
active  control  system  were  chosen  based  upon  the  criterion  that 
they  be  capable*,  in  principle,  of  functioning  in  the  space 
environment.  Several  types  of  electromagnetic  actuators  were 
designed  and  built.  Each  is  capable  of  inducing  force  or  torque 
on  a  structure  through  acceleration  of  captive  masses.  By 
operating  on  inertial  principles,  these  devices  function 
independent  of  the  laboratory  (earth)  frame.  Inertial  and 
relative  sensors  were  used  to  eliminate  the  need  for  state 
measurements  relative  to  the  laboratory  reference.  These 
space-realizable  components  were  used  to  control  the  elastic 
behavior  of  a  quasi  free-free  beam. 

Active  control  actuators  were  built  with  variable  passive 
dampers  and  springs  to  incorporate  into  these  devices  the  passive 
damping  characteristics  typical  of  passive  vibration  absorbers. 
Optimal  tuning  laws  were  developed  to  provide  maximum  passive 
damping  in  a  single  beam  mode  and  several  beam  modes.  In  the 
latter  case,  to  introduce  maximum  damping  into  several  modes,  the 
stiffness  of  the  actuator  should  be  tuned  to  the  lowest  frequency 
mode  and  the  actuator  damper  adjusted  for  desired  response. 


Two  optimizations  were  used  to  design  the  actuators.  First, 
a  sequential  optimization  was  performed  which  used  the  vibration 
absorber  tuning  laws  to  maximize  passive  damping  followed  by  a 
Linear  Quadratic  Regulator  optimization  to  select  the  active 
feedback  gains.  Secondly,  a  simultaneous  optimization  was 
performed  which  used  the  LQR  formulation  to  derive  gains  for 
various  passive  actuator  configurations  and  an  iterative  search 
was  used  to  determine  which  combined  set  of  passive  and  active 
actuator  characteristics  minimized  the  quadratic  state  and  control 
effort  cost  of  the  system.  Both  the  sequential  and  simultaneous 
optimizations  resulted  in  similar  actuator  designs  implying  that 
the  vibration  absorber  tuning  laws  can  be  used  to  select  the 
passive  characteristics  of  the  optimal  active  control  actuator. 
Finally,  optimal  actuation  locations  were  determined  and  simple 
collocated  feedback,  distributed  control  was  performed. 

An  experimental  test-bed  was  designed,  for  this  portion  of 
the  research,  which  exhibits  many  of  the  detrimental 
characteristics  of  large,  flexible  space  structures.  The  highly 
flexible  free- free  beam  possesses  a  low  wave  propagation  to 
computation  speed  ratio  in  addition  to  low  frequency,  large 
amplitude  vibrations.  The  geometric  dimensions  resulted  in  low 
frequency,  closely  spaced  modes.  Yet,  the  beam  has  linear 
properties  which  were  beneficial  in  developing  this  first  phase  of 
tests. 

The  combined  use  and  interaction  of  two  zero  level  and  a  one 
level  processor  was  demonstrated  in  the  framework  of  a  finite 


control  element.  A  robust  collocated  controller  was  demonstrated 


which  performed  independently  of  a  system  model.  Theoretically, 
vibration  absorber  damping  should  achieve  5 %  critical  damping 
through  the  addition  of  0.5%  beam  mass.  Experimentally,  4.2% 
critical  damping  was  achieved  through  the  addition  of  2.3%  of  the 
beam  mass.  The  active  control  actuators,  without  vibration 
absorber  characteristics,  achieved  2%  to  3%  critical  damping. 
These  actuators,  along  with  the  two  zero  level  processors  and  one 
component  level  processor,  comprize  the  present  state  of  the 
finite  control  element. 

There  are  three  sets  of  experiments  planned  for  the  near 
future.  The  first  set  involves  improving  the  performance  of  the 
tests  detailed  in  Chapter  5.  This  includes  improving  on  some  of 
the  techniques  used  in  signal  processing,  measurement  feedback  and 
beam  excitation  and  solving  some  of  the  control  and  friction 
limitations  encountered. 

The  second  set  of  tests  will  use  the  feedback  of  more 
collocated,  state  information  to  more  actuators.  In  these  tests, 
the  actuators  will  also  be  tuned  as  passive  vibration  absorbers. 
The  microprocessor  hardware  will  be  expanded  to  include  seven 
processors  distributed  on  three  levels  of  control.  Tests  will  be 
performed  which  include  free  decays,  frequency  sweeps,  impact 
testing  and  multimode  control. 

The  third  set  of  tests  will  attempt  to  exploit  the  wave 
propagation  behavior  of  the  beam.  Controllers  will  be  devised  to 
simulate  elastic  and  viscous  foundations  to  generate  propagation 
cutoff  frequencies  and  damped  behavior.  The  transmissivity  and 
reflectivity  of  the  beam  will  be  passively  and  actively  altered. 
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Finally,  disturbance  rejection  schemes  may  be  employed  to  suppress 
propagating  behavior. 

Thus  far,  this  research  has  demonstrated  the  feasibility  of 
incorporating  machine  intelligence  into  the  design  of  space 
structures.  Currently  available  technology  has  been  used 

effectively  in  a  space-realizable  experiment.  Though  passive 
techniques  currently  exhibit  the  same  performance  as  active 
feedback,  the  problems  experienced  in  the  active  tests  are 
surmountable.  Active  control  still  provides  a  versatility  in 
design  not  available  through  passive  methods.  Inevitably,  the 
optimal  form  of  structural  static  and  dynamic  control  will  be 
achieved  through  a  combination  of  various  techniques  of  which 
active  control  will  be  the  most  innovative. 
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A. 0  APPENDIX  A:  OPTIMIZING  THE  PASSIVE  PARAMETERS  OF  PASSIVE 

VIBRATION  ABSORBERS  AND  ACTIVE  CONTROL  ACTUATORS 

This  appendix  addresses  selection  of  optimal  spring  and 
damper  values  for  inertial-reaction  devices.  Two  device  types  are 
considered:  a  strictly  passive  inertial-reaction  vibration 
absorber,  and  an  active  control  actuator,  which  contains  passive 
elements  implementing  feedback  control.  The  optimal  values  for 
the  passive  components  are  determined  for  both  devices  when 
mounted  to  one  and  two  degree  of  freedom  structures. 

Three  optimizations  can  be  performed.  First  the  passive 
components  of  the  vibration  absorber  can  be  determined  based  upon 
the  minimization  of  a  weighted  quadratic  cost.  This  cost  is  the 
temporal  integral  of  state  penalties  times  the  response  of  the 
system  (Case  C  of  Table  2.4).  The  second  optimization  involves 
the  simultaneous  design  of  the  passive  actuator  components  and 
active  feedback  gains  using  Linear  Quadratic  Regulator  Theory  to 
determine  the  gains  for  each  set  passive  actuator  components  and 
performing  a  search  through  component  sets  to  determine  which 
control  actuator  components  and  gains  minimize  a  given  cost  (Case 
E  of  Table  2.4).  The  third  case  involves  using  LQR  theory  to 
sequentially  derive  a  set  of  feedback  gains  for  the  optimal 
passive  design  developed  in  Section  2.1.1  (Case  D  of  Table  2.4). 
The  passive  and  active  portions  of  each  design  can  then  be 
compered.  Section  A.l  presents  the  system  models  used  in  this 
analysis.  Section  A. 2  presents  parameter  optimization  using  a 


quadratic  cost,  and  Section  A.  3  details  the  results  and 
conclusions . 


A. 1  SYSTEM  MODELS 

Linear  models  were  used  for  both  the  control  actuator  and  the 
passive  vibration  absorber  paralleling  those  used  in  Chapter  2. 
The  passive  absorber  model  is  shown  in  Figure  A. la  and  the  control 
artuj'.or  model  is  shown  in  Figure  A.  lb.  The  actuator  or  passive 
absorber  was  mounted  to  two  separate  spring  and  mass  structures. 
The  first  structure  was  a  spring-mass  oscillator  that  had  a  single 
mode  of  vibration.  The  second  structure  was  a 
spring-mass-spring-mass  oscillator  with  two  vibration  modes.  The 
frequency  ratio  between  the  low  and  high  vibration  modes  of  the 
second  structure  was  2.618.  This  ratio  was  thought  to  be  small 
enough  to  represent  two  adjacent  modes.  The  mounting  location  on 
the  single  mode  structure  for  the  actuator  or  passive  absorber  was 
at  the  mass  (Figure  A. 2a).  The  actuator  or  passive  absorber  was 
mounted  to  the  two  DOF  structure  at  both  the  outer  mass  (Figure 
A. 2b)  and  at  the  inner  mass  (Figure  A. 2c).  Both  mounting 
locations  were  considered  because  of  different  modal  influence 
levels  at  the  two  masses.  The  influence  is  greater  on  the  low 
mode  at  the  outer  mass  and  the  influence  is  greater  on  the  high 
mode  at  the  inner  mass.  For  ease  of  reference,  the  systems  in 
Figures  A. 2a,  A. 2b,  and  A. 2c  will  be  referred  to  as  System  1, 


System  20  (outside  mount),  and  System  21  (inside  mount), 


respectively.  System  numbers  may  have  a  postscript  indicating  the 


type  of  vibration  damper  being  used.  A  postscript  of  "a" 
specifies  an  active  control  actuator  and  a  postscript  ”p" 
specifies  a  passive  absorber. 

The  state  differential  equation  describing  the  motion  of 
System  1  is 


x  -  Ax  +  Bu 
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and 


u  =  0  {for  passive  absorber} 


u  =  -Fx  {for  active  actuator}. 
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where  F  are  the  optimal  regulator  gains.  The  following 


substitutions  were  made  to  Equations  A.l  through  A. 6 


P0  -  K/M  DECOUPLED  STRUCTURAL  FREQUENCY 

PQ  -  k/m  DECOUPLED  ABSORBER  FREQUENCY 

a 

6  =  Pa/P0  ABSORBER/STRUCTURE  FREQUENCY  RATIO 

V  -  c/(2mp0)  NONDIMENS IONAL  DAMPING 
P  =  m/M  AflSORBER/STRUCTURE  MASS  RATIO 

d  =  1/K  STRUCTURE’S  STATIC  DEFLECTION 

z(n)  -  q(n)/d  NONDIMENS IONAL  POSITION 
i.iw  -  q(n)/'4p0  NONDIMENS  IONAL  VELOCITY 
z(n)  =  q(n)/dp,  NONDIMENS IONAL  ACCELERATION 
u’  =  u/ (Mdpp )  NONDIMENS IONAL  CONTROL  FORCE 

The  resulting  demens ionless  state  equations  are 
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where  F  are  the  nondimensiorial  optimal  regulator  gains. 

Now  the  fundamental  matrix  A’  (Equation  A. 11)  and  control 
effectiveness  matrix  B’  (Equation  A. 12)  are  functions  of 
dimensionless  quantities  6,  n,  and  /3,  which  represent  spring, 
damper,  and  mass  values  of  the  absorber  or  actuator.  For  this 
investigation  the  absorber/structural  mass  ratio  /}  was  set  to  0.05 
and  6  and  n  were  left  as  the  free  variables  to  be  used  in  the 
optimization.  At  this  value  the  absorber  mass  is  much  lighter 
than  the  structure  but  large  enough  to  assure  noticeable  damping 
in  the  systems.  Substituting  the  value  of  P  =  0.05  into  Equations 
A. 11  and  A. 12  yield 


-(H62/20)  P62 


(A. 14) 


which  is  the  final  nondimensional  form  of  the  equations  for  System 
1 . 

Nondimensionalizing  the  state  equation  for  System  20  and 
System  21  and  substituting  (i  -  0.05  yields  the  following  state  and 
control  effectiveness  matrices: 
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(A. 17) 


(A. 18) 


(A. 19) 


A. 2  QUADRATIC  COST  OPTIMIZATION  OF  PARAMETERS 


Values  for  the  frequency  ratio  and  nondimensional  damper  will 
be  derived  in  this  section  for  the  models  developed  in  the 
previous  section.  The  values  which  minimize  a  cost  based  on 
system  response  will  be  chosen  as  the  optimal  absorber  component 
settings.  In  order  to  evaluate  the  effectiveness  of  a  particular 
selection  of  6  and  p  (spring  and  damper  parameters)  for  a 
vibration  damper  design,  a  technique  is  needed  to  analyze  the 
structure  plus  damper  system  performance.  A  commonly  accepted 
method  of  measuring  performance  is  the  quadratic  cost  of  bringing 
the  system  to  rest  from  a  perturbed  state.  The  quadratic  cost  is 
defined  by 


J  =  /  (zTW  z  +  2zTW  u  +  uTW  u]  dt 

zz.  ru  uu 

1  o 

where 

W  -  state  penalty  matrix 
zz 

W  -  cross  cost  penalty  matrix 
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W  -  control  penalty  matrix 
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(A. 20) 


Optimal  full  state  feedback  control  (optimal  regulator  theory)  is 

/ 

based  on  minimization  of  the  cost  represented  by  Eq.  A.  20.  Thus  / 

by  using  the  cost  formulation  to  select  the  passive  device 


parameters,  the  active  feedback  gains  and  passive  parameters  could 


be  optimized  in  a  consistent  manner.  Now  this  cost  formulation 
will  be  used  to  simultaneously  optimize  both  the  passive 
parameters  and  active  gains  of  the  active  actuator,  then  secondly 
to  optimize  the  passive  parameters  of  the  vibration  absorber. 

Simultaneous  optimization  of  the  active  actuator  parameters 
for  Systems  la,  20a,  and  21a  involved  an  iterative  procedure.  The 
iteration  had  two  steps.  Lu  the  first  step  a  value  was  selected 
for  6  and  n,  and  the  gain  matrix  F*  for  optimal  state  feedback 
control  was  calculated  for  these  fixed  values  of  6  and  u.  The 
feedback  gain  matrix  F’  was  actually  calculated  for  a  discrete 
time  version  of  the  System  Equations  9  and  10b.  In  order  that  the 
discrete  time  feedback  gain  matrix  closely  match  the  continuous 
time  matrix  F ’  ,  a  period  of  O.Olw^  seconds  was  used  for  equation 
discretization.  The  second  step  was  to  calculate  the  quadratic 
cost  for  a  given  initial  condition  using 


}'  Hiccat  i  solution  matrix 
z0  -  initial  condition 

The  Hiccat i  solution  matrix  is  a  output  of  step  one,  the  feedback 
optimization  step,  while  the  initial  condition  was  a  predetermined 
vector  value.  This  quadratic  cost  calculation  was  tabulated. 
This  procedure  was  repeated  over  a  range  of  6  and  a  range  of  /i  and 
the  parameter  value  pair  that  resulted  in  the  lowest  quadratic 
cost  was  defined  to  be  the  optimal  passive  parameter  values. 
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Parameter  optimization  for  the  passive  absorber  parameters  in 
Systems  lp,  20p,  and  2Ip  followed  the  same  process  as  the  above 
actuator  parameter  optimization  described  above.  Unlike  the 
actuators,  passive  dampers  do  not  have  to  go  through  feedback 
optimization.  Yet  the  simplest  method  for  calculating  quadratic 
cost  was  from  the  Riccati  solution  matrix  (Equation  A. 21).  To 
take  advantage  of  this  fact  the  control  effectiveness  matrix  B’ 
(Equation  A.  17  or  A.  19)  was  set  to  zero.  With  the  system 
unaffected  by  control  action,  the  optimal  feedback  gain  matrix  F’ 
becomes  zero,  satisfying  Equation  A. 10a.  The  Riccati  solution 
matrix  will  still  be  generated  so  the  quadratic  cost  could  be 
calculated  as  before  (Eq.  A. 20).  Again,  the  quadratic  costs  were 
tabulated  for  a  range  ol  values  for  both  6  and  and  the  value 
pair  that  resulted  in  the  lowest  cost  was  defined  to  be  the 
optimal . 

Two  important  factors  in  the  quadratic  optimization  process 
need  further  discussion:  the  weighting  matrices  that  were  used  in 
the  quadratic:  cost  equation  (Equation  A. 20);  and  the  initial 
conditions  used  in  cost  calculations  (Equation  A. 21). 

The  criteria  used  to  determine  the  values  for  the  three 
weighting  matricies  were  as  follows.  First,  for  all  three  systems 
the  energy  storage  components  of  the  structure  were  penalized 
equally.  For  System  1  the  structural  energy  Es  is 

Eg  =  Mv?/2  +  Kx?/2  (A. 22) 
Nondimens ional izing  Equation  A. 22  by  Equations  A.7a-j  results  in 
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The  structural  energy  for  Systems  20  and  21  in  demensional  and 
nondimensional  forms  are 


Mvf/2  +  Mv2 /2  +  Kx?/2  +  K(x2  -  xx)2/2 


(A. 25) 
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The  coefficients  in  Equation  A. 23  (System  1)  or  in  Equation  A. 26 
(Systems  20  and  21)  were  incorporated  into  the  state  penalty 
matrix.  To  avoid  extra  restrictions  in  selection  of  optimal 
passive  parameters  or  feedback  gains,  the  energy  in  the  vibration 
damper  was  not  penalized.  This  is  reasonable  only  if  the 
resulting  damper  design  does  not  have  zero  spring  and  damper 
values  or  exhibit  extreme  resonant  peaks.  The  cross  cost 
state-control  penalty  matrix  was  set  to  zero.  Lastly, 
requirements  had  to  be  balanced  in  the  selection  of  the  control 
penalty  had  to  balance  two  conflicting  desires.  On  one  hand  the 
control  penalty  should  be  low  enough  to  allow  significant  control 
levels  that  will,  presumably,  affect  optimal  passive  parameter 
selections.  On  the  other  hand  the  control  penalty  should  be  high 
enough  to  keep  control  forces  small  compaired  to  structural  forces 


(a  practical  limit  for  space  applications). 


The  pentil ty  matrices  that  were  used  to  satisfy  the  above 
criteria  were 
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(A. 27) 

(A. 28) 

(A. 29) 


(A. 30) 

(A. 31  ) 

(A. 32) 


(In  addition  to  the  above  state  weighting  matrix  for  Systems 
20  and  21  a  state  weighting  matrix  with  1,  2,  1  in  the  first 
four  diagonal  elements  and  0  elsewhere  was  also  used.  The  results 
using  either  state  weighting  matrix  were  very  similar  so  the 


results  of  this  latter  case  will  not  be  presented.) 


The  second  factor  which  requires  consideration  is  the  choice 


of  initial  conditions.  Passive  parameters  were  optimized  for  the 
systems  under  several  sets  of  initial  conditions.  The  initial 
conditions  were  all  chosen  to  represent  an  impulsive  energy 
imparted  to  the  structure  when  it  was  at  rest.  The  initial  state 
vector  therefore  had  zero  positions  and  nonzero  velocities.  The 
impulse  was  assumed  to  have  excited  only  homogeneous  structural 
modes  and  not  the  vibration  absorber.  Table  A.l  lists  the  ten 
parameter  optimization  runs  that  were  done  and  their  corresponding 
initial  condition  sets.  For  optimizations  1  and  2  the  single 
structure  mass  had  an  initial  velocity  and  all  other  states  were 
zero.  For  optimizations  3  through  10  three  types  of  initial 
conditions  were  used.  The  initial  energy  for  optimizations  3  and 
4  was  in  the  low  frequency  mode  of  the  structure.  Similarly,  the 
initial  energy  was  in  the  high  frequency  mode  of  the  structure  for 
runs  5  and  6.  For  optimizations  7  through  10  the  initial 
condition  was  the  sum  of  the  previous  two  initial  conditions 
resulting  in  equal  energy  in  both  structural  modes. 


A. 3  RESULTS  AND  EVALUATION 

Following  the  procedure  outlined  in  the  previous  section, 
optimal  values  of  the  passive  parameters  were  determined  using  an 
iterative  search  procedure.  The  optimal  passive  parameter  values 
were  selected  after  generating  tables  of  quadratic  cost  versus 
values  of  6  (spring)  and  p  (damper).  Initially,  "coarse"  cost 


tables  were  formed  that  spanned  a  wide  range  of  6  and  p.  The  6 
and  p  value  pair(s)  that  produced  a  cost  minima  in  this  table  was 
noted.  "Fine"  cost  tables  were  formed  "around"  minima  over  a 
narrower  range  of  parameter  values  to  determine  more  accurately 
the  optimal  values  for  6  and  p.  Table  A. 2a  is  a  "coarse"  cost 
table  for  System  ip,  the  single  degree  of  freedom  passive 
vibration  absorber  only  case  (optimization  1  in  Table  A.l),  with 
one  minimum  at  a  6  and  p  of  0.968  and  0.106.  Table  A. 2b  is  a 
"fine"  cost  table  for  the  same  system  where  the  final  optimal 
parameter  selection  for  6  and  p  is  0.970  and  0.107  respectively. 
Similar  "coarse"  (Tables  A. 4  13)  and  "fine"  (Tables  A.  14  23)  cost 
tables  were  fornu  1  for  all  tne  optimizations  in  Table  A.l.  Only 
one  minimum  was  found  for  each  optimization.  The  results  of  the 
optimizations  are  summarized  in  Table  A. 3.  The  optimal  values  of 
6  and  the  damping  ratio  p/6  are  recorded  next  to  the  system  and 
initial  condition  they  were  optimized  for.  Table  A. 3  is  split  in 
two  so  the  results  for  passive  absorbers  and  active  actuators  can 
be  compared  side  by  side.  The  initial  condition  column  of  Table 
A. 3  has  three  symbols;  "1”,  for  initial  energy  in  the  low 
stuctural  mode  only,  "H" ,  for  initial  energy  in  the  high 
structural  mode  only,  and  for  initial  energy  equally 
distributed  between  both  modes. 

The  passive  component  selections  for  the  optimal  actuator  and 
passive  damper  designs  have  some  interesting  trends.  The 
resulting  natural  frequency,  pQ  -  6p^  was  always  near  a  mode 
frequency  of  the  structure.  The  natural  frequency  for  both  damper 
types  in  System  1  is  just  below  the  structural  frequency  w  . 


System  20  has  two  structural  nodes;  a  low  frequency  node  at 
0.618w^  and  a  high  frequency  mode  at  1.62w&.  When  initial  energy 
was  only  in  the  low  mode  the  natural  frequencies  for  both  actuator 
and  passive  damper  were  just  below  the  low  mode  frequency. 
Likewise,  the  natural  frequencys  were  just  below  the  high  mode 
frequency  for  initial  energy  in  the  high  node  only.  The 
similarity  between  this  low  or  high  mode  tuning  and  the  tuning  in 
System  1  was  probably  due  to  the  small  damper  to  structural  mass 
ratio  P>.  This  mass  ratio  keeped  the  coupling  between  the  two 
structural  vibration  modes  in  System  20  and  21  low.  The  trend  of 
the  damping  ratio  p/6  was  toward  light  damping  when  the  initial 
energy  was  in  one  structural  mode.  The  net  result  of  these 
component  selections  are  vibration  dampers  that  were  resonant  at  a 
frequency  near  the  structural  mode  with  contained  the  initial 
system  energy. 

Actuator  and  passive  damper  designs  that  resulted  from 
optimizations  with  equal  initial  energies  in  both  structural  modes 
had  somewhat  different  characteristics.  The  natural  frequencies 
of  both  devices  were  near'  the  low  frequency  structural  mode.  For 
both  systems  20p  and  21p  the  passive  absorber  Systems  20p  and  2lp 
the  natural  frequencies  were  close  to,  but  above,  the  low 
frequency  mode.  The  natural  frequencies  of  the  two  actuators  were 
below  the  low  modal  frequency  but  above  the  natural  frequency  of 
the  atuator  optimized  with  just  low  mode  initial  energy.  The  more 
noticable  effect  of  optimizing  with  both  modes  initially  excited 
is  the  substantia]  increase  in  the  required  damping  ratios.  Also 
note  the  difference  in  damping  ratios  between  System  20p  and 


System  21p.  This  difference  was  the  most  dramatic  effect  of 
different  damper  mounting  locations. 

After  the  ten  passive  parameter  optimizations  (see  Table  A.l) 
were  completed  the  pole  locations  for  the  resulting  closed-loop 
systems  were  calculated.  The  pole  locations  were  originally 
calculated  in  the  Z-plane  from  the  discrete  time  system  equations. 
The  Z-plane  poles  were  translated  into  S-plane  poles  by  the 
equat ion 


S  :  lnCzl/l  (A. 33) 

Because  of  the  short  sample  period  of  0.01wq  sec.,  these  pole 
locations  should  not  differ  appreciably  from  poles  calculated  from 
the  continuous  time  system. 

Five  additional  actuator  systems  were  designed  using  a 
sequential  design  procedure.  First  the  optimal  passive  parameters 
were  determined.  Then  with  the  spring  and  damper  parameters  6  and 
p  constrained  to  optimal  passive  absorber  values,  optimal 
regulator  feedback  gams  were  calculated.  This  will  be  refered  to 
as  the  sequential  or  suboptimal  active  actuator  design.  The 
closed-loop  poles  were  then  calculated  as  they  were  for  the 
optimal  active  actuator  and  optimal  passive  absorber  systems. 

The  closed-loop  poles  for  the  optimal  passive  absorber, 
optimal  active  actuator,  and  corresponding  suboptimal  actuator 
systems  of  each  row  in  Table  A.  3  are  plotted  in  Figures  A.  3  to 
A. 7.  For  instance,  Figure  A. 3  shows  the  pole  locations  for  System 
1  (first  row  of  Table  A. 3)  with  either  an  optimal  active  actuator, 
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an  optimal  passive  absorber,  or  a  suboptimal  actuator  with  6  and  n 


constrained 

to 

0.970 

and 

0. 107  (see 

Table 

A. 2b). 

It  is 

interesting 

t  o 

observe 

the 

closed- loop 

pole 

locations 

of  the 

optimal  and  suboptimal  actuator  systems.  When  the  initial  energy 
of  the  structure  is  in  one  mode,  the  pole  locations  for  both 
optimal  and  suboptimal  actuator  systems  are  nearly  the  same 
(Figures  A. 3,  A. 4  and  A.^i.  Similar  results  are  evident  for  two 
oi  me  poies  oi  systems  20a  and  Zia  when  optimizations  were  done 
with  both  strucural  modes  having  equal  initial  energies  (Figures 
A. 6  and  A.  7 , .  The  one  pole  that  was  significantly  different  in 
tin  two  actuator  systems  is  dominated  by  motion  of  the  actuators. 


A.  4  CONCLUSION’ 

Two  important  insight;-  can  be  drawn  from  these  result.-.. 
First,  actuator  design  might  be  done  sequentially  by  optimizing 
the  passive  components  and  then  optimizing  the  feedback  control 
law.  The  differences  in  system  performance  between  sequential  and 
parallel  optimization  of  passive  parameters  and  control  (as  judged 
by  pole  locations)  appear  to  be  small.  Second,  passive  parameters 
might,  be  selected  to  damp  the  lowest  mode  of  significant  influence 
at  the  damper  mounting  point.  Improving  damping  over  higher  modes 
is  achieved  by  larger  damping  ratios  with  the  same  natural 
frequency.  If  both  these  insights  prove  correct,  a  damper  design 
procedure  could  be  fairly  straightforward. 
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B.O  APPENDIX  B:  INFORMATION  TRANSPORT  IN  STRUCTURES 


Explicitly  or  implicitly,  all  active  control  schemes  require 
an  accurate  plant  model  to  assure  effective  high  authority 
control.  When  using  classical  or  output  feedback,  this  model  is 
used  implicitly,  in  that  it  was  necessary  for  the  design  of  the 
compensation.  In  estimator  based  control,  the  model  is  used  not 
only  in  design,  but  is  explicitly  used  in  the  real-time 
processing.  When  attempting  to  devise  algorithms  for  controlling 
flexible  structures,  several  difficulties  arise  in  the 
implementation  of  these  control  techniques.  These  difficulties 
stem  from  the  large  number  of  degrees  of  freedom  exhibited  by 
elastic  systems. 

For  practical  purposes,  an  elastic  structure  is  a  complex 
dynamic  system  of  almost  infinite  order.  The  demand  created  by 
traditional  approaches  to  high  authority  control  challenge  the 
modeling  of  this  complex,  infinite  order  system  in  two  ways. 
First,  the  process  requires  a  system  model  which  includes  many 
degrees  of  freedom.  Secondly,  the  accuracy  of  the  model  must  be 
very  good,  even  in  higher  order  modes.  For  these  reasons,  this 
chapter  investigates  alternative  methods  by  which  the  structure 
can  be  modeled  accurately  without  the  use  of  high  order  models. 

Specifically,  attention  will  be  focused  on'  the  transport  of 
information  in  structures.  The  two  questions  that  will  be 
addressed  are:  over  what  regions  of  a  structure  is  the  information 
that  a  force  has  been  applied  instantaneously  transmitted  (i.e. 


what  is  the  instantaneous  local  sphere  of  influence  of  a  force) 
and  secondly;  once  that  force  has  established  a  displacement 


disturbance;  over  what  range  is  that  disturbance  eventually 
transported  in  a  damped  structure. 

This  chapter  contains  two  discussions  of  structural  modeling 
to  address  these  two  questions  of  instantaneous  and  eventual 
information  transport.  The  first  is  a  discretized  elasticity 
model  which  would  lead  to  a  traditional  standing  mode  model.  The 
second  is  a  travelling  wave  or  wave  propagation  treatment. 

Several  options  are  available  to  discretize  the  elastic 
behavioi  of  a  structure.  If  the  structure  is  sufficiently  simple 
to  describe  exactly  with  a  partial  differential  equation  then 
exact  flexibility  and  stiffness  influence  coefficients  can  be 
derived.  Alternatively,  finite  difference  operators  can  be  used 
to  discretize  the  partial  differential  equation.  If  the  structure 
is  sufficiently  complex,  finite  elements,  which  exactly  or 
approximately  model  the  behavior  of  the  structure  in  a  piecewise 
manner,  can  be  employed  to  discretize  the  structure.  Using  the 
example  of  a  simple  uniform  beam,  these  three  approaches  will  be 
used  to  delineate  the  sphere  of  influence  of  static,  elastic 
structures . 

The  second  model  formulation  is  related  to  the  propagating 
transportation  of  information  in  a  structure.  Structural  response 
can  be  treated  in  a  standing  wave  fashion,  where  the  lag  between 
disturbance  and  response  is  infinitesimal,  or  treated  using  a  wave 
propagation  approach.  The  factor  that  determines  which  approach 
is  more  applicable  is  the  relative  propagation  to  decay  times  of 


a 


the  disturbance.  In  other  words,  are  the  group  velocities  slow 
enough,  or  the  structure  large  enough,  that  the  controller  can 
dissipate  the  propagating  waves  before  they  reflect  and  reinforce 
(or  constructively  interfere)  to  generate  a  standing  wave  motion? 

The  motivation  for  studying  the  information  transport 
arguments  is  to  understand  the  analytic  structural  mechanics  basis 
for  specially  discretizing  the  control  problem  into  actively 
controlled  elements  or  components,  much  like  the  structural 
analysis  problem  is  discretized  into  finite  elements.  If  the 
control  of  large,  flexible  structures  can  be  performed  using  an 
assemblage  of  controlled  elements,  several  advantageous  features 
would  result.  First,  the  controller  would  exhibit  robust  behavior 
because  the  complex  system  model  would  simply  be  a  combination  of 
many  simpler,  theoretically  more  accurate  models  which  reflect  the 
principles  involved  in  basic,  local  structural  behavior. 
Secondly,  implementation  of  the  structural  control  lei  could 
parallel  the  construction  of  the  space  structure  in  a  building 
block  fashion.  Thirdly,  the  real-time  computational  task  would  be 
significantly  reduced  since  the  control  response  can  be  computed 
on  either  a  local  or  component  level,  and  a  higher  or  global 
level,  and  the  matrices  involved  in  the  global  level  would  be 


expected  to  be  sparse. 
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B.l  DISCRETIZATION  OF  LOCAL  ELASTIC  BEHAVIOR 


In  considering  the  extent  of  instantaneous  transport  of 
information  in  a  discretized  structure,  one  could  consider 
coupling  in  the  three  structural  dynamic  matrices.  However,  it  is 
well  known  in  finite  element  modeling  that  simple  diagonal,  lumped 
mass  models  give  results  comparable  to  more  fully  populated, 
consistant  mass  matrices.  Furthermore,  from  first  principles,  it 
is  the  localized  inertia  which  effects  the  dynamics  of  a 
structure.  Thus,  if  the  generalized  coordinates  are  chosen  in  an 
appropriate  way,  the  coupling  in  a  structure  can  be  largely 
eliminated  from  the  mass  matrix.  Since  the  coupling  in  the 
structural  damping  matrix  is  a  reflection  of  that  in  the  stiffness 
matrix,  a  study  of  the  real  stiffness  matrix  will  be  sufficient  to 
identify  the  sphere  of  influence  and  coupling  patterns  present  in 
a  structure. 

In  studying  this  coupling,  finite  elements,  finite 
differences  and  flexibility  influence  coefficients  will  be 
examined.  See  Reference  3  for  basic-  theory  behind  finite  element 
analysis.  The  governing  equation  for  the  static  response  of  a 
structure,  using  any  of  these  discretization  methods,  is 

{q  f  =  fC]  (Q I-  (B.la) 

{01  -  j  K  !  { q  *  1 B . lb) 


In  these  matrix  equations,  C  is  the  flexibility  matrix,  K  is  the 


stiffness  matrix,  q  is  the  generalized  displacement  vector  and  Q 
is  the  generalized  force  vector.  When  using  a  discrete  method, 


the  structure  is  conceptually  divided  into  segments,  areas  or 
elements  whose  boundaries  occur  at  nodes.  The  generalized 
displacements  and  forces  are  defined  to  occur  at  these  nodes  (Fig. 
B.l).  The  fineness  of  the  discretization  and  the  number  of 
allowable  nodal  displacements  determine  the  number  of  degrees  of 
freedom  and  therefore  the  order  of  the  matrices  in  Eqs.  B.l. 

This  stiffness  matrix  can  be  constructed  using  several 
methods.  Construction  by  flexibility  influence  coefficients, 
finite  difference  operators  and  finite  elements  will  be  discussed. 
The  charat.  t.eristics  of  each  formulation  will  be  compared  and  the 
implication  on  information  transport,  analyzed. 

The  flexibility  influence  coefficient  matrix  is  formed  for  a 
constrained  structure  by  combining  vectors  which  contain  the 
displacements  at  all  of  the  i^  locations  due  to  a  unit  load  at 
the  location.  By  computing  the  vector  for  a  load  at  each 

location,  a  populated  flexibility  matrix  results.  The  flexibility 
influence  coefficient  matrix  will,  in  general,  be  fully  populated, 
because  all  i ^  locations  displace  when  a  force  is  placed  at  any 

location.  However,  the  flexibility  matrix  is  not  the 
appropriate  representation  of  the  structure’s  elasticity  to 
consider  when  addressing  the  assessment  of  spheres  of  influence 
for  the  control  problem.  The  control  problem  is  not  posed  as: 


given  a  force  at  A,  how  will  it  affect  the  structure.  Rather,  the 
control  problem  is  posed  as:  given  a  required  response  at  B,  what 
forces  must  be  applied  to  achieve  it.  This  more  closely  parallels 


the  function  of  the  stiffness  influence  coefficient  matrix. 


The  stiffness  matrix  is  the  inverse  of  the  flexibility 
matrix.  Eq.  B.2  gives  the  Bernoulli  Euler  differential  equation 
governing  static,  elastic  beam  bending  where  F(x)  is  the 
distributed  force,  w  is  the  transverse  displacement,  E  is  the 
modulus  of  elasticity  and  I  is  the  bending  moment  of  inertia. 


F(x)  =  — r 
dx" 


EI(x) 


d2w(x) 

dx2 


(B.2) 


Using  this  governing  differential  equation,  the  stiffness  matrix, 
derived  by  inversion  of  the  flexibility  matrix  for  a  uniform 
clamped-clamp^d  beam  with  nine  internal  translational  degrees  of 
freedom  (ten  "elements")  is 
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A  moderately  large  number  of  "elements"  was  chosen  to  reduce  the 
effect  of  boundary  conditions  at  the  Lot.ernai  nodes.  This 
influence  coefficient  method  yields  an  exact  discrete  solution, 
within  the  limitations  of  beam  theory,  but  is  restricted  to  simple 
structures . 

Finite  diffence  operators  provide  an  alternative  approach  to 
approximate  discretization  of  the  governing  differential  equation 


(Eq.  B.2). 


Assuming  constant  El,  finite  difference  operators  can 


be  used  to  approximate  the  fourth  derivative  spatial  dependence  of 
transverse  displacement.  The  lowest  order  difference  operator  for 
this  equation  is  second  order  accurate 

A  ,  , 

--4  =  q_2  -  4 . Oq_j  +  6.0qo  -  4-Oqj  +  q£  (B.4) 
dx  "  ' 

where  the  q’s  are  the  w(x)  displacements  at  the  nodes.  Table  B.l 
lists  appropriate  fourth  order  operators  of  increasing  truncation 
orders  and  illustrates  the  relative  magnitude  of  the  terms  in 
each.  These  FD  operators  are  approximate  discrete  representations 
of  the  fourth  derivative  using  displacements  at  internal 
locations.  The  effects  of  boundary  conditions  are  not  yet 
included.  A  portion  of  a  resulting  stiffness  matrix,  using  the 
second  order  accurate  FD  operator,  is 
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(B.5) 


Note  that  the  exact  stiffness  matrix  (Eq.  B.3),  although 
diagonally  dominant,  is  fully  populated.  The  approximate 
stiffness  matrix  is  banded.  Thus,  for  identical  discretization 


mesh  sizes,  stiffness  matrices  derived  using  FD  operators,  because 


of  their  sparseness,  require  less  computation  for  inversion  than 
those  derived  using  flexibility  influence  coefficient  methods. 
But,  the  FD  operators  do  not  provide  the  same  degree  of  accuracy. 

In  both  the  stiffness  influence  coefficients  and  finite 
difference  operator  formulations,  the  relative  magnitudes  of  the 
stiffness  matrix  entries  in  a  given  column  indicate  several 
important  trends.  Recall  that  a  column  of  the  stiffness  matrix  is 
the  set  of  forces  necessary  to  create  a  unit  displacement  at  the 


diagonal ,  and 

zeroes  elsewhere. 

Fig.  B.2a 

compares 

the 

distribution  of 

magnitude  of  the 

terms  in  one 

column  in 

the 

interior  of  the  stillness  matrix  derived  using  fourth  derivative 
spatial  finite  difference  operators  of  different  truncation 
errors.  The  entries  are  shown  with  reference  to  the  diagonal 
entry.  Fig.  B.2b  overlays  the  central  column  of  the  stiffness 
matrix,  derived  using  flexibility  influence  coefficients  (Eq. 
B.3i,  with  a  sixth  order  accurate,  fourth  derivative  finite 
difference  operator’. 

In  each  curve  of  both  figures  (Figs.  B.2a  and  B.2b),  the 
symbols  indicate  the  relative  elastic  influence  of  a  force  at  each 
nodal  location  on  the  location  represented  on  the  diagonal.  The 
influence  of  the  two  adjacent  forces  is  obviously  important.  At  a 
position  three  node  points  away,  the  influence  of  a  force  is  less 
than  ten  percent  of  the  influence  of  a  force  at  the  diagonal 
position.  At  a  position  five  node  points  away,  the  influence  of  a 
force  is  almost  insignificant.  This  rolloff  is  a  measure  of  the 
relative  effectiveness  of  noncollocated  forces  in  influencing 
static,  elastic  behavior. 
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Thus,  the  following  observations  can  be  made  on  the 
bandedness  of  the  stiffness  matrix  when  transverse  displacements 
are  the  only  deflections  represented: 


A.  When  the  elasticity  is  represented  exactly  (by  stiffness 
influence  coefficients)  the  stiffness  matrix  is  not 
banded,  but  is  diagonally  dominant.  At  least  two  off 
diagonal  terms  must  be  maintained,  but  certainly  not  more 
than  four  or  five,  in  any  truncation. 

B.  When  the  elasticity  is  represented  approximately  (by  FD 
approximations  to  the  governing  differential  equation)  the 
stiffness  matrix  is  banded.  The  width  of  the  band  depends 
on  the  order  of  the  FD  approximation.  At  a  minimum,  a 
second  order  accurate  approximation  is  necessary  (two  off 
diagonal  terms).  Rarely  in  finite  difference  calculation 
is  more  than  a  fourth  order  accurate  scheme  (three  off 
diagonal  terms)  used. 

A  third  approach  to  the  discretization  of  Eq.  B.2  uses  finite 
elements.  A  two  node  finite  element  is  show  in  Fig.  B.3.  For  a 
uniform  beam,  the:  static  equations  of  Bernoulli  Euler  beam  theory 
can  be  satisfied  on  the  interior  of  the  element  by  a  third  order 


polynomials  used  as  interpolation  functions  for  this  element. 
Thus,  this  formulation  is  exact  and  is  the  same  as  would  result  by 
using  flexibility  influence  coefficients  if  translations  and 
rotations  were  the  generalized  displacements.  The  corresponding 
element  stiffness  matrix  is 
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Assembling  these  elements,  as  shown  in  Fig.  B.l,  yields  the  highly 


banded  stiffness  matrix 
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0 

-12 
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0 
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-6L 

24 

The  reason  fhe  finite  element  stiffness  matrix,  which  is 
exact  in  this  case,  is  strictly  banded,  as  opposed  to  the  exact 
stiffness  influence  coefficient  matrix  which  is  only  diagonally 
dominant,  is  that  the  FE  formulation  includes  rotational  degrees 
of  freedom  at  the  nodes.  The  terms  in  a  column  of  the  stiffness 
influence  coefficient  matrix  are  the  forces  required,  at  each 
location,  to  create  a  unit  displacement  at  the  i*"*1  location  and 
zero  displacement  at  ail  other  nodal  locations.  Figure  B.4 
graphically  represents  the  distribution  of  nodal  forces.  At  all 
nodes  other  than  that  exhibiting  unit  displacement,  displacement 
is  zero  but  slope  is  nonzero.  Thus,  the  vector  is  fully 
populated. 

Figure  B.5a  and  B.5b  illustrate  the  same  situation  for  the 
finite  element  formulation  with  generalized  nodal  displacements  of 
displacement  and  rotation  and  generalized  nodal  forces  of  force 
and  torque.  By  the  definition  of  the  columns  of  a  stiffness 
matrix,  displacement  and  rotation  are  zero  at  locations  adjacent 
to  the  unit  displacement  position.  In  a  flexural  structure,  if 
the  displacement  and  slope  are  zero  at  the  adjacent  node,  the 
structure  is  stress  free  beyond  these  adjacent  locations  and  the 


corresponding  stiffness  terms  are  zero.  Therefore,  if  a  certain 
nodal  static  shape  is  desired  between  a  pair  of  nodes,  an  exact 
formulation,  which  includes  only  translational  displacements, 
requires  that  the  structure  be  influenced  at  all  nodal  locations. 
However,  with  a  formulation  that  includes  both  translational 
displacements  and  rotations,  the  availability  of  control  forces 
and  torques  reduces  the  required  actuation  effort  to  the  node  of 
interest  and  the  two  adjacent  nodes.  The  following  observation 
can  therefore  be  made  on  the  bandedness  of  the  stiffness  matrix 
when  transverse  displacements  and  rotations  are  the  deflections 
represented: 

C.  Even  when  the  elasticity  is  represented  exactly,  the 
stiffness  matrix  is  banded. 

There  are  certain  implications  of  the  effect  of  this 
bandedness  on  controller  formulations.  Collocated  control  uses 
sensors  to  create  control  commands  for  actuators  located  at  the 
same  position.  In  other  words,  local  information  is  used  locally. 
The  argument  above  regarding  the  bandedness  can  be  interpreted  in 
the  following  way:  the  further  away  an  actuator  is  from  a 
particular  location,  the  less  ability  it  has  in  influencing 
elastic  behavior  at  that  location.  Thus,  a  slight  generalization 
of  collocated  control  to  include  feedback  amongst  adjacent  nodes 
may  be  a  natural  way  of  dividing  the  control  problem. 
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B . 2  WAVE  PROPAGATION  MODEL  FORMULATION 


An  understanding  of  wave  propogation  models  can  provide 
further  insight  which  may  be  exploitable  in  the  formulation  of  a 
structural  controller.  A  wave  propogation  model  is  simply  an 
alternate  view  of  a  standing  wave  representation.  If  the 
transverse  motion  of  a  beam  is  characterized  by  an  isolated  wave 
packet  propagating  with  some  finite  phase  velocity,  sensors  on 
those  segments  of  the  beam  not  moving  need  not  be  monitored.  The 
existance  of  an  isolated  wave  packet  is  particularly  noticeable 
during  the  early  stages  of  the  generation  of  standing  wave  motion. 
A  discussion  of  wave  propagation  can  be  found  in  Reference  6  and 
wave  propagation  applied  to  the  modeling  and  control  of  large 
space  structures  can  be  found  in  Reference  5. 

Both  standing  and  propagating  wave  formulations  are  obtained 
from  the  same  mathematical  model  of  the  structure.  The  equation 
governing  the  transverse  motion  of  a  thin  Bernoulli  Euler  beam  is 


1 

2 

a 


— —  F(x, t) 
m(x) 


where 


(B.8a) 


2  EI(x) 

a  -  - 

m(x) 


(B.8b) 


For  the  homogeneous  case  an  allowable  solution,  corresponding  to 
rightward  and  leftward  travelling  waves,  is 
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wix,t)  -  A  (  sin(kx-ut)  +  sin(kx+wt)  ] 


(B.9a) 


Through  angle  summation  identities,  Eq.  B.9a  can  be  transformed 
into  Eq.  B.9b  which  gives  the  separation  of  variable  solution 
typical  of  standing  wave  motion. 

w(x,t)  =  2  A  sinkx  coswt  (B.9b) 


For  further  analysis,  the  more  generalized  Timoshenko  beam 
model  is  used.  tijs.  L.iGd  and  B.lOb  are  the  governing  equations 
where  G  and  k  are  the  shear  modulus  and  shape  factor, 
respectively. 
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The  section  coordinates  are  transverse  displacement  (w)  and  cross 
sectional  rotation  (lp).  Several  interesting  conclusions  can  be 
drawn  from  the  analysis  of  propagating  disturbances  in  a  damped 
Timoshenko  beam  presented  in  Reference  8. 

In  the  analysis  performed  in  Reference  8,  damping  is  modeled 
as  complex  extension  and  shear  modulus: 


E  --  EM  i  \y ) 


(B. 11a) 


G  -  G( 1  +  i? ) 


(B. 11b) 


The  damping  contribution  from  extensional  and  shear  effects  is 
assumed  identical  and  small.  With  the  introduction  of  damping,  a 
characteristic  length  for  dissipation  can  be  defined  which 
corresponds  to  the  length  over  which  a  wave  has  traveled  when  its 
amplitude  has  decayed  to  about  37%  (1/e)  of  the  original  amplitude 
(i.e.  the;  length  traveled  in  one  damping  time  constant).  This 
dissipation  characteristic  length,  as  a  function  of  wavelength  and 
beam  damping  at  resonance,  is 


L  - 


A 

4rr$ 


1  +  4? 


■*  1 


(B. 12) 


where 


/  is  wavelength 
$  is  critical  damping  ratio 

Consider  the  implication  of  this  relationship  on  the 
propagation  of  a  disturbance  in  a  damped  structure.  Immediately 
after  a  force  impulse  is  exerted  on  a  structure,  waves  travel 
along  the  structure,  reflect  and  interfere  with  waves  traveling  in 
the  opposite  direction  until  a  steady-state  motion  is  achieved. 
If  the  structure  is  damped,  these  waves  can  be  significantly 
dissipated  before  a  large  steady-state  motion  occurs.  Fig.  B.6  is 


a  plot  of  the  ratio  of  the  dissipation  characteristic  length  to 


wavelength  as  a  function  of  damping  ratio,  defined  in  Eq.  B.12. 


If  a  particular  motion  exhibits  a  damping  ratio  of  $  =  0.707 
(Log{$)  =  -0.15)  during  free  decay,  a  traveling  wave,  at  the  same 
frequency,  will  decay  to  within  21%  of  its  original  magnitude  in 
less  than  0.31  wavelengths  (At  this  value  of  damping  ratio,  Fig. 
B.6  gives  the  distance  traveled  in  one  time  constant  to  be  L/A  = 
0.3075  (Log(0.3075)  =  -0.5122)).  To  be  less  than  5 X  of  original 
amplitude,  the  wave  must  have  propagated  for  three  time  constants 
and  three  dissipation  characteristic  lengths,  or  0.93  wavelengths. 
Thus,  in  less  than  one  wavelength,  a  wave  propagating  in  a 
structure  with  $  -  0.707  will  effectively  disappear! 

The  wavelengths  of  interest  in  a  free- free-  beam  are*  those  of 
the  free-free  normal  modes,  since  their’s  is  the  motion  which  is 
most  easily  generated.  fig.  B.7  shows  the  fraction  of  beam  length 
travelled  for  a  modal  wave,  starting  at  one  end,  to  decay  to 
within  fit  of  its  original  magnitude  (i.e.  three  dissipation 
characteristic  lengths).  These  curves  are  functions  of  modal 
damping  ratio.  Modes  1,  2,  3,  10,  20  and  50  are  represented.  The 
portion  of  a  curve  below  the  horizontal  reference  line  (i.e.  below 
Log(3L/Beam  Length)  -  0.0)  represent  damping  ratios  for  which  that 
mode  fails  to  travel  one  length  of  the  beam  before  its  amplitude 
has  decreased  to  less  than  5 %  of  its  original.  One  length  of 
travel  occurs  at  (  -  0.955  for  the  first  mode,  $  =  0.477  for  mode 
two,  ?  -  0.318  for  mode  three,  0.095  for  mode  ten,  $  =  0.047  for 
mode  twenty  and  $  =  0.019  for  mode  fifty.  It  is  therefore 


reasonable  to  expect  that  high  frequency  disturbances  will  not 
propagate  the  length  of  the  structure  when  a  reasonable  amount  of 


Table  2.1;  Eigenvectors  for  Two  DOF  Model  for  a  Mass  Ratio  of  0 
0.02 


Frequency 

Mode 

Inner 

Outer 

Relative  Motion 

Ratio _ 

Number 

Mass 

Mass 

Across  Damper 

0.1000 

1 

1.0 

4950.00 

4949.00 

2 

1.0 

-  0.01 

1.01 

0.9804 

1 

1.0 

8.74 

7.74 

2 

1.0 

-  5.72 

6.72 

10.0000 

1 

1.0 

1.01 

0.01 

2 

1.0 

-49.51 

50.51 

Table  2.2;  Parameter 

Values  for  Tuning  to 

the  Low  and  High  Mode 

of  a  Two 

Mode  Structural  Model 

Mode  XN] _ 

1 

2 

0  =  M  /M 

n  no 

1 . 3820 

3.6180 

0  ■-  Ma/M^ 

0.0276 

0.0724 

*tr,  '  M>/K„ 

0.0200 

0.0200 

6*"oi'i 

0.9804 

0.9804 

Optimal  Frequency 
Ratio  based  on  £ 

7 

n 

0.6180 

1.6180 

Modal  Frequency 

6opt 

0.6060 

1.5863 

Optimal  Frequency 
Ratio  Based  on  0 

Table  2.3; 

Eigenvectors  and  Frequencies  for 

Different 

Absorber 

F  requency 

Rat ios . 

(Three  DOF  Model) 

Absorber 

Mode 

Modal 

Eigenvector  of 

Mass 

Frequency 

Frequency 

Rat  in 

1 

2 

3 

0.  1000 

1 

0.100 

0.0003 

0.0006 

1.0000 

2 

0.618 

-  23.0044 

-  37.2171 

1.0000 

3 

1.618 

421.9374 

-260.8053 

1.0000 

0.6060 

1 

0,570 

0 . 0567 

0 .  i  i  50 

1.0000 

o 

n  r,-*.  7 

(i  lino 

n 

1.0000 

3 

1.619 

9.8783 

-  6.1374 

1.0000 

1 . 1000 

1 

0.609 

0.4256 

0.6932 

1.0000 

J. 

1.112 

-  0.0289 

-  0.0221 

1.0000 

3 

1.623 

1.8538 

-  1.1780 

1.0000 

1.600" 

1 

0.611 

0.5257 

0.8552 

1.0000 

•') 

1.553 

-  0.1582 

0.0650 

1.0000 

3 

1 . 693 

0.1287 

-  0.1114 

1.0000 

2. 1000 

1 

0.611 

0.5628 

0.9152 

1.0000 

2 

1.604 

-  0.7262 

0.4164 

1.0000 

3 

2.  141 

0.0153 

-  0.0396 

1.0000 

20.0000 

1 

0.612 

0.6146 

0.9991 

1.0000 

2 

1.612 

-  1.6599 

0.9935 

1.0000 

3 

20.275 

0.0001 

-  0.0277 

1.0000 

Tab  I< 


!.4;  Actuator  Design  Cases  Studied 
Appendix  A. 


Steady  State 
Mj  nimax 

Pole 

Placement 

Passive 
Vibration 
Absorber  onlv 

Section 

2.1.1  (A) 

Section 
2.1.2  (B) 

Sequential 

Choice  of 

Passive 

Actuator 

Design  and 

Active 

Feedback 

Gains  * 

Section  2.3 

(D) 

(Section 
2.3.1  plus 
Appendix 

A. 2) 

not 

performed 

Simultaneous 
Choice  ot 

Pass i ve 

Actuator 
Design  and 

not 

not 

Act  l v< 

Feedback 

Gains  * 

performed 

performed 

in  Chapter  Two  and 


Minimum 
Quadratic  Cost 

Section 
2.1.3  (C) 


Section  2.3 
(Appendix  A. 2) 


Section  2.3 

(E) 

( Section  2.1.3 
plus  Appendix 
A. 2) 


*  Feedback  gains  derived  using  LQR  formulation 


Table  4.1;  Modal  Characteristics  of  a  Free-Free  Beam. 


E  = 

15.9xl06lbs/in2 

1.090x10 

^2Nt/m  Young’s 

Modulus 

I  - 

0.0009766  in4 

4.065x10 

-10  4 

m  Area  Moment  of  Inertia 

01  = 

0.007329  slugs/in 

4.120  kg/m 

Mass  per 

Unit  Length 

1  = 

288.0  in 

7.32  m 

Length 

b  = 

4.0  in 

0.1016  m 

Width 

t  = 

0.125  in 

0.003175 

m 

Thickness 

M  = 

2.10  slugs 

30.7  kg 

Total  Mass 

Mode  Modal  Frequency 

Modal 

Mass 

Group  Velocity 

w/4  actuators _ 

wo/actuaiors 

wo/actuators 

# 

Exper imantal 

Frii-  Unit- 

0\<5r>1 

Analytic 

( Hz  ) 

At  The  Tip 

At  The  1/4 

(m/sec) 

(Kg; 

Point  (kg) 

1 

0.33 

7.61 

694.81 

8.9 

•) 

0.63 

7.61 

21.02 

17.8 

3 

1.06 

7.61 

18.99 

26.7 

4 

1.88 

7.61 

112.96 

35.6 

5 

2.77 

7.45 

96.91 

44 . 5 

6 

3.76 

7.49 

17.39 

53.4 

7 

5 . 23 

7.51 

17.54 

62.3 

8 

6.70 

”.53 

102.33 

71.2 

9 

8.20 

7  55 

103.31 

80.1 

10 

10.06 

7.56 

17.74 

89.0 

]  1 

12.44 

7.57 

17.72 

97.9 

12 

14.27 

7.58 

103.51 

106.8 

13 

16.36 

7.59 

103.61 

115.7 

14 

19.38 

7.59 

17.77 

124.6 

lf» 

21.98 

7.60 

17.78 

133.5 

16 

22.37 

7.60 

103.83 

142.4 

Table  4.2;  Force  Actuator  Specifications  and  Components. 


Maximum  Actuator  Force 
Actuator  Mass  wo/Proof-Mass 
Proof-Mass  Mass  f 1 ) 

Arm  Length 
Motor  Resistance 
Motor  Inductance 
Motor’s  Peak  Amperage 
Motor’s  Peak  Voltage 
Torque  Constant 
Peak  Stall  Torque 
Field  Magnet 
DC  Motor 

Flexural  Spring 
Damper  Adjustable- 

Displacement  Transducer 
Velocity  Transducer 
Current  Source 


3. U  Nts. 

0.550  Kg. 

0.165  Kg. 

0.127  m 
5.000  Ohms 
1.760  mH 
4 . 920  Amps . 

19.10  Volts 
35.70  mN  m/Amp 

173.0  mN  m 

Samarium-Cobalt 
Pittman  Corp. 

Model  #  7214 
Adjustable 
Airpot,  Inc. 

Model  #  S95«!.5Y2 
See  Table  4.4 
Se-  Table  4.4 
EGkU  Torque  System 
110  Watt 
Model  #  PA-223 


Table  4.3;  Sensor  Specifications. 


Company: 

Accel eromet ers 
Endevco 

Tachometers 
Inland  Motor 

Potentiometers 

Maurey 

Quantity: 

4 

4 

4 

Type: 

Piezoresistive 

A 

25 , 000  Ohms 

Mode]  #: 

2262-25 

TG-0702 

SA  200-M608 

Range: 

+/~  25  g 

+/-  363  rad/sec 

+/-  178  Degs 

Sw? '  *  ivity: 

nn  p  ,,v  •  - 

0  •; /rad/sec 

0  TJ  V/Deg 

Resonant 
frequency : 

2.5- 2.8  KHz 

Table  4.4: 

Signal  Processing  Equipment 

Specif icat ions . 

Anti-Aliasing 

F l Iters 

Smoothing 

Filters 

Voltage 

Amplifiers 

Type: 

Sallen-key 

Sal ien-keye 

Instrumentation 

Output 

Voltage- 

Range: 

+/-  10  V 

+/-  10  V 

+/-  10  V 

Selectable 

Gains. 

- 

x  10,  50,  100, 
200,  Variable 

Corner 

Frequency: 

700  Hz 

100  Hz 

- 

Damping 

Ratio: 

0.7071 

0.7071 

[S 


Table  5. 

Mode 

.  1 ;  Open  and 
Modes . 

Freq  (Hz . ) 
O.L.  C.L. 

Closed- Loop  Characteristics  of  Controlled 

Damping  Ratio  Fractional  Change 

O.L.  C.L.  In  Damping  Ratio 

C.L./O. L. 

7 

5.2  5.1 

.0059  .0185 

3.14 

8 

6.7  6.7 

.0024  .0200 

8.33 

9 

8.3  8.2 

.0033  .0216 

6.55 

10 

10.3  10.3 

.0046  .0165 

3.59 

11 

12.5  12.3 

.0025  .0215 

8.60 

12 

14.5  14.5 

.0015  .0284 

18.90 

13 

16.5  16.5 

.0081  ? 

? 

14 

19.8  19.5 

.0069  .0113 

1.64 

Ave: 

.0044  .0197 

7.25 

.0197/0. 

0044-4.48 

O.L.  =  Open-Loop 

C.L.  =  Closed-Loop 

Table 

5.2;  Analytic  Results  for  Optimally  Tuned 
Absorber . 

Vibration 

Mode 

Absorber 

Effective 

Absorber/ 

Calculated 

Calculated 

Position 

Modal 

Mass 

(Kg.) 

Modal 

Mass  Ratio 

Opt imal 
Frequency 
Rat  10 

Modal 

Damping 

Ratio 

. . ■ 

- . —  —  ■ 

(Eq.  2.13) 

(Eg. .  5.12.) 

/ 

L  4 

18.089 

0.009 

0.991 

0.0  Jq 

H 

TIP 

8.081 

0.020 

0.980 

0.050 

9 

TIP 

8.096 

0.020 

0.980 

0.050 

Table  5.2; 

Experimental  Vibration  Absorber 
Lsing  a  Tip  Absorber  at  Various 

Results  for  8th  Mode 
Frequency  Ratios. 

Absorber/ 

Beam 

Frequency 
Rat  10 

Damping 

Ratio 

Fractional 
Increase 
in  Damping 
Ratio 

Figure 

( Baseline! 

0.0024 

_ 

5.5a 

0.387 

0.0031 

1.29 

5.13a 

0.493 

0.0082 

3.42 

5.13b 

0.671 

0.0117 

4.88 

5.13c 

0.846 

0.0300 

12.50 

5.13d 

0.884 

0.0120 

17.50 

5. 13e 

151 


(O  CO 


Table  5.4;  Summary  of  Experimental  Vibration  Absorber  Results  for 
Modes  Tested. 

Mode  Absorber  Frequency  Untuned  Tuned  Fractional  Figure 


Position 

Rat  io 

Damping 

Damping 

Increase  in 

Ratio 

Ratio 

Damping 

Ratio 

L/4 

0.859 

0.0059 

0.0120 

2.03 

5.  14 

Tip 

0.884 

0.0024 

0.0420 

17.50 

5.  I3e 

Tip 

0.864 

0.0033 

0.0170 

5.  15 

5.  15 

0.0039 

0.0237 

8.23 

Ratio  of  Experimental 
Damping  Ratio  to 
Analytic- 
Damping  Ratio 


Table  A.l;  Parameter  Optimization:  Initial  Conditions. 


Opt.  System  Initial  Condition  (z  ) 


_  9 

Zl 

_ _ 

z2. 

_ 2 _ 

o 

jz  3 

.  z  3  ... 

1 

lp 

0 

1.000 

0 

0.000 

- 

- 

2 

la 

0 

1.000 

0 

0.000 

- 

- 

3 

20p 

0 

1.000 

0 

1.618 

0 

0.000 

4 

20a 

0 

1.000 

0 

1.618 

0 

0.000 

5 

20p 

0 

-1.618 

0 

1.000 

0 

0.000 

6 

20a 

0 

-1.618 

0 

1.000 

0 

0.000 

7 

20p 

0 

-0.618 

0 

2.618 

0 

0.000 

8 

20a 

0 

-0.618 

0 

2.618 

0 

0.000 

9 

2Ip 

0 

-0.618 

0 

2.618 

0 

0.000 

10 

21a 

0 

-0.618 

0 

2.618 

0 

0.000 

Table  A. 2;  Quadratic  Cost  Versus  Frequency  Ratio  (6)  and 
Nondimensional  Damper  (u) .  System  lp  with 

z  =  ■'  0.0  0.0  1.0  0.0  jT 


6  n  values 


values : 

0. 

07000 

0. 

08607 

0. 

10580: 

0 

13010 

0 

16000 

0 

.  750  : 

31. 

650 

26. 

940 

23. 

380 

20 

820 

19 

150 

0 

.852 

16. 

730 

14. 

800 

13. 

510 

12 

790 

12 

620 

0 

.968 

9. 

654 

9. 

045 

8. 

-*• 

cc 

8 

984 

9 

528 

1 

.  100 

20. 

820 

18. 

130 

16. 

210 

14 

990 

14 

410 

1 

.250  : 

69. 

070 

57. 

400 

48. 

140 

40 

940 

35 

520 

6 

U  values 

values : 

0. 

10200 

0. 

10470 

0. 

10740: 

0 

11010 

0 

11300 

0 

.  960 

8. 

867 

8. 

865 

8. 

860 

8 

863 

8 

870 

0 

.965 

8. 

847 

8. 

837 

8. 

833  : 

8 

836 

8 

844 

0 

.970 

8. 

837 

8. 

827 

8. 

824* 

8 

827 

8 

835 

0 

.  975 

8. 

847 

8. 

837 

8. 

833 

8. 

836 

8 

844 

0 

.  980  , 

8. 

878 

8. 

867 

8. 

862 

8 

865 

8 

872 

*  minimum 


Table  A. 3;  Optimal  Passive  Parameters. 

Vibration  Damper  Type 


System 

Initial 

Passive 

Absorber 

Acti ve 

Act  uab 

Cond i t ion 

& 

.  .  M/& 

6 

M/6 

1 

0 . 970 

0.1110 

0 . 920 

0.0802 

20 

!. 

0.605 

0.0950 

0.540 

0 . 0528 

20 

H 

1.570 

0.0573 

1.540 

0.0377 

20 

0 . 690 

0.5130 

0.576 

0 . 085 1 

21 

0.630 

0.2510 

0.602 

0 . 0525 

Table  A. 4; 


Cost  Map,  System  Ip 


COARSE"  Quadratic  Cost  Table. 
.Passive  Vibration  Absorber) 

Initial  State  Vector  -  [0.  0. 


H  VALS 

6  VALS 

0.07000 

0.08607 

0.7500 

31.650 

26.940 

0.8522 

16.730 

14.800 

0 . 9682 

9.654 

9.045 

1.1000 

20.820 

18. 130 

1.2500 

69.070 

57.400 

0 

10580 

0 

13010: 

0 

16000 

23 

380 

20 

820  : 

19 

150 

13. 

510  ; 

12 

790 

12 

620 

8 

825 

8 

984  : 

9 

528 

16 

210  ; 

14 

990  : 

14 

410 

48 

140  : 

40 

940 

35 

520 

Table  A. 5;  "COARSE"  Quadratic  Cost  Table. 

Cost  Map,  System  la  Control  Actuator) 

Initial  State  Vector  -  0.  0.  1.  O.j^ 


U  VALS 

6  VALS 

0.04000 

0.05264. 

0.06928: 

0.09118: 

0.  12000 

0.7500 

7.095 

7 . 056 

7.035 

7.055 

7.  151 

0.8522 

6.283 

6.243 

6.221  . 

6.244 

6.350 

0 . 9682 

6.224 

6.177  ; 

6.145 

6.  154 

6.243 

1.1000 

7 . 393 

7.340 

7.300 

7.293 

7 . 355 

1 . 2500 

10.370 

10.310 

10 . 260 

10.230 

10.260 

Table  A. 6;  "COARSE”  Quadratic  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber; 

Initial  State  Vector  •  [0.  0.  0.  1.  1.1318  0.  T 


(j  VALS 

6  VALS 

_ 0.033C _ 0.0435 

0,0575 

0.0758 

0.1000 

0.5000 

217.80  :  174.80 

145.50 

127.40 

119.20 

0.6300 

81.00  71.23 

66.96 

67.87 

74.03 

0.7937 

828.90  :  638.70 

497.00 

393.70  : 

321.30 

1.0000 

4463.00  : 3402. 00 

2588.00 

1974.00 

1517.00 

1.2600 

0. 170E+05 : 0. 123E+05 

9468.00 

7205.00 

5479.00 

1.5870 

0.577E*-05;0. 4 14E+05 

0 . 297E  <-05 

0 . 227E+05 : 0 . 172E*05 

2.0000 

0. 124E+06 : 0. 104E+06 

0.800Ef05 

0. 687E+05  0.51 1E+05 

154 


Table  A. 7;  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator) 

Initial  State  Vector  -  (0.  0.  0.  1.  1.618  0.-^ 


n  VALS 

0.01000 

0.01732 

0.03000 

0.05196 

0.90000 

6  VALS 
0.5000 

32.83 

32.70 

32.65 

33.01 

34.81 

0.6300 

35.  14 

34.99 

34.86 

35.02 

36.39 

0.7937 

54.00 

53.84 

53.67 

53.72 

54 . 70 

1.0000 

103.40 

103.20 

102.90 

102.70 

102.90 

1.2600 

194.00 

193.80 

193.40 

192.90 

192.40 

1.5870 

341.30 

341.10 

340.60 

339.90 

339.00 

2.0000 

575.00 

574.70 

574.20 

573.50 

572.60 

Table  A. 8:  "COARSE"  Quadratic:  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber)  Initial  Condition  Vector 

;  (0.  0.  0.  -1.618  1.  o.;r 


!J  VALS 

6  VALS 

....  0.0300 

0.0520 

0.0900 

0. 1559; 

0.2700 

0.5000 

3562.00 

2060.00 

1201.00 

715.60 

152.20 

0.6300 

3114.00 

1802.00 

1052.00 

629.40 

402.80 

0.7937 

2460.00 

1427.00 

836.00 

505.20 

331.10 

1.0000 

1595.00 

926.40 

547.90 

339.00  , 

235.  10 

1.2600 

615.70 

363.00 

222.70 

151.30 

126.70 

1.5870 

68.43 

47.08 

40.32 

46.02  : 

65 . 96 

2.0000 

1727.00 

1003.00 

592.10 

364.40 

219.70 

Table  A. 9;  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator)  Initial  State  Vector 

=  t0.  0.  0.  -1.618  1.  0./ 


fi  VALS 

6  VALS 

0.03000 

0.05196 

0.09000 

' '  6  .15590 

0.27000 

0.5000 

87.05 

86.44 

85.60 

84.79 

85.19 

0.6300 

81.47 

80.87 

80.07 

79.37 

80.08 

0.7937 

72.57 

72.00 

71.29 

70.85 

72.  14 

1.0000 

58.64 

58.  16 

57.68 

57.77 

60.22 

1.2600 

38.62 

38.33 

38.35 

39.63 

44 . 38 

1.5870 

26.84 

26.59 

26.79 

28.71 

35.05 

2.0000 

61.41 

61.01 

60.72 

61. 15 

63 . 93 

Table  A. 10;  "COAfiSE"  Quadratic  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber)  Initial  State  Vector 

=  [0.  0.  0.  -0.618  2.618  0. ;T 


VALS 

0.  1000 

0.  1732 

0.3000 

0.5196 

0.9000, 

6  VALS 
0.5000 

1208.0 

782.8 

600.4 

603.9 

794.8 

0.6300 

1033.0 

681.8 

542.0 

570.2 

774.8 

0.7937 

1094.0 

717.0 

562.2 

582. 1 

781.4 

1.0000 

2045.0 

1268.0 

879.7 

765.4 

886.9  : 

1.2600 

5748.0 

3389.0 

2106.0 

1473.0 

1297.0 

1.5870 

0.  174E+05 

9952.0 

5883.0 

3637.0 

2547.0 

2.0000 

0. 520E+05 

0 . 286E+05 

0. 166K*05 

9674.0 

5948.0 

Table  A.  1.1;  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator)  Initial  State  Vector 

-  i 0.  0.  0.  -0.618  2.618  0.jr 


p  _  VALS 

6  VALS 

0.0150 

0. 0260 

0.0450 

0.0779 

0. 1350 

0.5000 

121.5 

121.1 

120.7 

121.2 

124.6 

0.6300 

120.6 

120.  1 

119.6 

119.7 

122.3 

0.7937 

132.2 

131.7 

131.2 

131.1 

133.2 

1.0000 

168.4 

167.8 

167.2 

166.9 

168.2 

1 . 2600 

241.2 

240.6 

239.9 

239.3 

239.9 

1.5870 

377.3 

376.4 

375.7 

374.7 

374.6 

2.0000 

636.6 

636.0 

635. 1 

634.2 

633.3 

Table  A. 12;  "COARSE"  Quadratic  Cost  Table.  Cost  Map,  System  2ip 
(Passive  Vibration  Absorber)  Initial  State*  Vector 

=  [0.  0.  0.  -0.618  2.618  0. ; 1 

p  VALS  0.0500  0.0866  0.1500.  0 . 2598  0.4500 

6  VALS  : 

0.5000  : 1218.0  :  786.3  .  597.6  ;  594.3  :  774.6 

0.6300  :  821.8  557.5  :  465.6  :  517.9  730.8 

0.7937  1800.0  1122.0  791.5  705.9  .  839.2 

1.0000  : 7477.0  4380.0  ; 2670.0  1787.0  : 1463.0 

1.2600  . 0. 270E+05 : 0. 157E+05 : 9071 . 0  :5490.0  :3588.0 

1.5870  :  0. 816E*05  :  0.486E-t-05  :  0. 275E-*-05  :  0. 164E+05  :  9847.0 

2.0000  .  0 . 259E+06 ;  0 . 155E )06  :  0 . 804 E 4 05  0 . 471E  +  05  :  0 . 281E+05  . 


Table  A.  13,  "COARSE"  Quadratic  Cost  Table.  Coit  Map,  System  21a 
.Control  Actuator)  Initial  State  Vector 

=  lO.  0.  0.  -0.618  2.618  0. )T 


p  VALS 

0.0100 

0.0173 

0.0300 

0.0520 

0.0900 

6  VALS 

0.5000 

99.88 

99.62 

99.49 

100.10 

103.30 

0.6300 

93.58 

93.28 

93 . 04 

93.44 

96.28 

0.7937 

117.50 

117.10 

116.80 

116.90 

118.70 

1.0000 

192.30 

191.90 

191.30 

190.70 

190.70 

1 . 2600 

325.80 

325.40 

324.70 

323.70 

322.90 

1.5870 

554.70 

554.30 

553.60 

552.30 

551.20 

2.0000 

952.50 

952.00 

951.10 

949.70 

947.70 

Table  A. 14;  "FINE"  Quadratic  Cost  Table.  Cost  Map,  System  Ip 
(Passive  Vibration  Absorber)  Initial  State  Vector 

=  (0.  0.  1.  0.]T 


p  VALS 

6  VALS 

0.1020 

:  0.1046 

0.1074 

0.1101 

0.1130 

0.960 

8.876 

8.865 

8.860 

8.863 

8.870 

0.965 

8.847 

8.837 

8.833 

8.836 

8.844 

0.970 

8.837 

8.827 

8.824 

8.827 

8.835 

0 . 975 

8.847 

8.837 

8.833 

8.836 

8.844 

0.980 

8.878 

8.867 

8.862 

8.865 

8.872 

Table  A. 15;  "FINE”  Quadratic  Cost  Table.  Cost  Map,  System  la 
(Control  Actuator)  Initial  State  Vector 

-  i o .  o .  i .  o.;T 


p  VALS 

6  VALS 

0.06800 

0 . 07082 

0.07376, 

0,07681 

0,08000 

0.900 

6.071 

6.070 

6.069 

6.070 

6.072 

0.910 

6.060 

6.059 

6.058 

6.059 

6.061 

0.920 

6.057 

6.055 

6.055  : 

6.055 

6.057 

0.930 

6.061 

6.059 

6.058 

6.058 

6.060 

0.940 

6.073 

6.071 

6.070 

6.070 

6.071 

Tab Le  A. 16; 


V  VALS 
6  VALS 
0.6025 
0.6037 
0.6050 
0.6062 
0.6075 


Table  A. 


iu  VALS 

6  VALS 

0.5200 

0.5297 

0.5396 

0.5497 

0.5600_ 


Table  A. 


H  VALS 
6  VALS 
1 . 565 
1.569 
1.572 
1.576 
1 . 580 


"FINE"  Quadratic  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber)  Initial  State  Vector 


=  (0. 

0. 

0. 

1.0 

1.618 

0. 

T 

:  ’  0 

05600 

d 

05674’ 

0 

05748 

0 

05824 

0 

05900; 

:  61 

11 

61 

10 

61 

09 

61 

10 

61 

11 

;  61 

06 

61 

05 

61 

04 

61 

05 

61 

06  : 

:  61 

05 

61 

04 

61 

03 

61 

04 

61 

05  : 

:  61 

06 

61 

05 

61 

04 

61 

05 

61 

06  : 

:  61 

10 

61 

09 

61 

08 

61 

09 

61 

10 

17;  "FINE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
Control  Actuator)  Initial  State*  Vector 

=  [0.  0.  0.  1.  1.618  0. ;T 


0.02200; 

0.02505: 

0.02853: 

0.03249: 

0.03700 

32.24 

32.22 

32.22 

32.23 

32.27 

32.14 

32.12 

32.11 

32.12 

32.  16 

32.10 

32.08 

32.07 

32.08 

32.11 

32.14 

32.12 

32.11 

32.11 

32.  14 

32.25 

32.22 

32.21 

nr) 

32.24 

18;  "FINE"  Quadratic  Cost  Table.  Cost.  Map,  System  20p 
(Passive  Vibration  Absorber)  Initial  State  Vector 

:  ; 0 .  0.  0.  1.618  1.  0. ! r 


0.08500 

0.08740: 

0.08986: 

0.09239: 

0.09500 

39.85 

39.79  : 

39.77 

39.78 

39.80 

39.77 

39.71 

39.69 

39.71 

39.73 

39.76 

39.70 

39.68 

39.70 

39.72 

39.82 

39.76  : 

39.74 

39.75 

39.77 

39.94 

39. 88  _  : 

39.86 

39.87 

39. 89 

t 

• 

/ 

4 

I 

1 


V-V- 


Table  A. 19; 


"FINE"  Quadratic  Cost  Table.  Cost  Map,  System  20a 
Control  Actuator/  Initial  State  Vector 

;  [0.  0.  0.  -1.618  1.  0. iT 


p  VALS 

"  6.04500 

0.05113 

~  "0.05809 

0.06601 

0.07500 

6  VALS 
1.530 

26.31 

26.27 

26.25 

26.26 

26.32  : 

1.542 

26.28 

26.24 

26.22 

26.23 

26.28 

1.555 

26.31 

26.27 

26.25 

26.26 

26.31  : 

1.567 

26.39 

26.35 

26.33 

26.34 

26.38 

1.580 

26.53 

26.49 

26.46 

26.47 

26.52  : 

Table  A. 20; 


"FINE"  Quadratic  Cost  Table.  Cost  Map,  System  20p 
(Passive  Vibration  Absorber)  Initial  State  Vector 


io.  0.  0.  -0.618  2.618  0.;T 


p  VALS 

6  VALS 

0.3300 

0.3418 

.0.3541 

:  0.3668 

.  0.3800  , 

0.6750 

527.9 

526.8 

526.  3 

526.5 

527.6 

0.6824 

527.5 

:  526.4 

: 525.9 

; 526.  1 

1527.3 

0 . 6898 

527.3 

526.2 

525.8 

525 . 9 

527.1 

0.6974 

527.3 

526.3 

: 525.8 

1526.0 

: 527 . 2 

0.7050 

527.7 

526.6 

526.  1 

526.7 

527.5 

Table  A. 21;  "FINE”  Quadratic  Cost  Table.  Cost  Map,  System  20a 
(Control  Actuator)  Initial  State  Vector 

-  0.  0.  0.  -0.618  2.618  0.'f 


p  VALS 

__0 .03000 

0.03834 

0.04899 

0 . 06260 

0.08000 

6  VALS 
0.5350 

;  119.6 

119.4 

119.3 

119.4 

119.8 

0 . 555 1 

:  119.2 

119.0 

118.9 

119.0 

1 1 9 . 3 

0.5759 

:  119.  1 

118.9 

118.7 

118.8 

1 19.  1 

0.5976 

:  119.3 

119.0 

118.9 

118.8 

1  19.  1 

0 . 6200 

119.7 

119.4 

119.2 

119.2 

119. 4 

Table  A.  22 ;  "FINE  Quadratic  Cost  Table.  Cost  Map,  System  2] 
(Passive  Vibration  Absorber)  Initial  State  Vector 

-  [0.  0.  0.  -0.618  2.618  0.jT 


p  VALS 
6  VALS 

0.1460 

0.  1520 

0.1593 

0. 1670 

0.1 

0.6200 

468.3 

466.6 

.465.5 

.465.7 

467.2 

0.6249 

467.5 

'  465 . 9 

1464.8 

:  465.0 

:466.6 

0.6299 

.467.2 

465.6 

464.5 

464.8 

466.3 

0.6349 

467 . 4 

;  465 . 7 

464.7 

464.9 

l  466 . 4 

0 . 6400 

468.0 

:  466.4 

.465.3 

465.5 

467.0 

Table  A. 28;  "FINK"  Quadratic  Cost  Table.  Cost  Map,  System 
(Control  Actuator’  Initial  State  Vector 

=  |0.  0.  0.  -0.618  2.618  0.;f 
u  VAI.S  0.02700  0.02921  0.02)61.  0.02-120  0.02700 


0.5900 

.  92.68 

.  92.67 

.  92.67 

92 . 69 

92.72 

0.5962 

.  92.61 

:  92.59 

92.59 

92 . 60 

.  92.62 

0.6024 

92.58 

.  92.56 

92 . 56 

92.57 

92.60 

0.6087 

:  92.60 

i  92.59 

l  92.58 

.  92.59 

-  92.62 

0.6150 

92.68 

:  92.66 

92.66 

92.67 

92.69 

ritWlMT 


Table  B.l;  fourth  Derivative  Finite  Difference  Operators  and 
Relative  Magnitudes  of  Terms  in  the  Operator. 


Order 


4 


-3 


-•> 


*0 


*1 


flo 


Operator 


MAGNITUDE  OF  RESPONSE  (X./A  )  VERSUS  NOND I  MF.NSIONAL 
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FIG.  2.5:  NONDIMENS IONAL  NATURAL  FREQUENCIES  (y>  OF  THE  TWO  SYSTEM 
MODES  AS  A  FUNCTION  OF  ABSORBER  FREQUENCY  RATIO  (6). 
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A:  RESPONSE  OF  ALL  SYSTEM  MODES 

MAGNITUDE  OF  RESPONSE  (X/A  )  VERSUS  NONDIMENS  1 ONA1,  FREQUENCY  (co/pQ)  FOR  AN  ABSORBER 
TUNED  TO  THE  HIGHER  STRUCTURAL  MODE  (CURVE  1)  AND  SOMEWHAT  TUNED  TO  BOTH  MODES 
(CURVE  2). 
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NDNDIMENSIQNAL  FREQUENCY 

DETAIL  IN  THE  VICINITY  OF  THE  HIGHER  MODAL  FREQUENCY 
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2.12:  RELATIVE  MODAL  MOTIONS  OF  MASS  ELEMENTS  FOR 

DIFFERENT  STRUCTURE? ABSORBER  FREQUENCY  RATIOS. 
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FIG.  3.3: 


A:  CONTROL  OPTION  PROVIDING  EXTENDED  BANDWIDTH  IN  THE 
LOW  AUTHORITY  CONTROLLER. 


B:  OVERLAPPING  LOW  AND  HIGH  AUTHORITY  CONTROL. 

OPTIONS  FOR  THE  SELECTION  OF  BANDWIDTHS  AND  INTERACTIONS 
OF  THE  VARIOUS  LEVELS  OF  CONTROL. 
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FIG.  3.4:  SINGLE  BUS  HARDWARE  ARRANGEMENT  FOR  MULTIPROCESSOR, 
HIERARCHIC  ARCHITECTURE. 


ELEMENT  1  ELEMENT  2 

FIG.  3.5:  MULTIBUS  HARDWARE  ARRANGEMENT  FOR  MULTIPROCESSOR, 
HIERARCHIC  ARCHITECTURE. 


PROCESSING  HARDWARE. 


RETURN  TO 
HUMAN 
INTERFACE 


FIG.  3.9:  FIRST  GENERATION  CONTROL  SOFTWARE  FLOW  CHART . 
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HUMAN  INTERFACE 


A:  TORQUE  CONFIGURATION 


B:  FORCE  CONFIGURATION 


C  :  SIDE  VIEW 

ELECTROMAGNETIC,  INERTIAL-REACTION  ACTUATOR 
CONFIGURATIONS . 


MODEL  OF  A  SINGLE  DOF  TORQUE  ACTUATOR  B: 
CONTROLLING  A  SINGLE  ROTATIONAL  MODE 
OF  A  STRUCTURE. 


MODEL  OF  A  — — - ►  0 

SINGLE  DOF  FORCE  ACTUATOR 
CONTROLLING  A  SINGLE  TRANS¬ 
LATIONAL  MODE  OF  A  STRUCTURE 


FIG.  4.4: 


MODELS  OF  ACTUATOR  CONFIGURATIONS  COUPLED  TO  A  SINGLE  DOF 
STRUCTURAL  REPRESENTATION . 


AD-A169  928  DEVELOPMENT  OF  FINITE  ACTIVE  CONTROL  ELEMENTS  FOR  LHR8C 
FLEXIBLE  SPACE  ST. .  <U)  MASSACHUSETTS  INST  OF  TECH 
CAMBRIDGE  HA  SPACE  SVSTEMS  LAB  D  H  MILLER  ET  AL. 


UNCLASSIFIED  82  JUN  85  MIT-SSL-6-85  AF0SR-TR-8S-84BB 


ivyaaasog  ^wflflaas  »&&&»  \&&&&&  \P£8& 


1.1 


*  ACTUATOR  LOCATIONS 


MODES  7-11 


— MODES  1,2 


NONOIMENSONAL  LENGTH 


THE  SUM  OF  THE  SQUARES  OF  THE  CONTROLLED  MODE  SHAPES 


nondiuensiokal  length 


FIG.  5 


THE  PRODUCT  OF  THE  SQUARES  OF  THE  CONTROLLED  MODE 
SHAPES. 

1:  ACTUATOR  EFFECTIVENESS  AS  A  FUNCTION  OF  POSITION  ON  THE 
BEAM  AS  INDICATED  BY  (A)  AND  (B) . 
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FIG.  5.3:  EXPERIMENTAL  INSTABILITY  II 
GAIN  BEAM  RATE  FEEDBACK  (S: 
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FIG.  5.4: 


B :  CLOSED-LOOP 

EXPERIMENTAL  TRANSIENT  DECAY  OF  SEVENTH  MODE. 
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B :  CLOSED-LOOP 

.  5.10:  EXPERIMENTAL  TRANSIENT  DECAY  OF  THIRTEENTH  MODE. 
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FIG.  5.12:  EXPERIMENTAL  FORCED  RESPONSE  VERSUS  SINUSOIDAL 
EXCITATION  FOR  OPEN  AND  CLOSED-LOOP  CONTROLLER. 
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TIME  (SEC) 

E:  ACTUATOR  FREQUENCY  =  5.93  HZ. 

EXPERIMENTAL  PASSIVE  TUNING  OF  VIBRATION  ABSORBER 
TO  THE  EIGHTH  MODE. 
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FIG.  5.14:  EXPERIMENTAL  PASSIVE  TUNING  TO  THE  SEVENTH  MODE. 


FIG.  5.15:  EXPERIMENTAL  PASSIVE  TUNING  TO  THE  NINTH  MODE. 
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10  11  12  13 


MINIMUM  ACHIFVKD  DAMP HR  SETT  INI 
FKI.Ut'ENCY  RATIOS. 


STRUCTURAL  REPRESENT AT I ONS  COUPLED  TO  ONE  DOF 
ABSORBER. 


CLOSED-LOOP  POLES  FOR  SYSTEM  20  OPTIMIZED 


CLOSED-LOOP  POLES  FOR  SYSTEM  20  OPTIMIZED  F 


21  OPTIMIZED  FOR  Z. 


MAG/(MAC  OF  DIAGONAL  ENTRY)  MAC/ (MAG  OF  DIAGONAL  TERM) 


TRUNCATION  ORDER 

TRIANGLE:  5 ^ 
DIAMOND:  4^ 
CROSS:  3rd 
SOUARE :  2nd 


SQUARES  :  STIFFNESS  MATRIX 
COLUMN  ENTRIES  DERIVED 
USING  A  SIX  ELEMENT 
FLEXIBILITY  INFLUENCE 
COEFF I C I ENT  FO RMULAT I  ON 


CROSSES:  FOURTH  ORDER 
FINITE  DIFFERENCE 
OPERATOR  CORRESPONDING  TO 
FOURTH  SPACIAL  DERIVATIVE 


OFF  DIAGONAL  POSITION 

COMPARISON  OF  BANDEDNESS  IN  FINITE  DIFFERENCE  OPERATORS 
AND  FLEXIBILITY  INFLUENCE  COEFFICIENTS. 

BANDEDNESS  IN  DISCRETE  REPRESENTATIONS  OF  ELASTIC 
STRUCTURES . 


